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We present new full-sky temperature maps in five frequency bands from 23 
to 94 GHz, based on data from the first three years of the WMAP sky survey. 
The new maps are consistent with the first-year maps and are more sensitive. 
The three-year maps incorporate several improvements in data processing made 
possible by the additional years of data and by a more complete analysis of the 
polarization signal. These include several new consistency tests as well as refi ne- 
ments in the gain calibration and beam response models (IJarosik et al.ll2006l ). 

We employ two forms of multi-frequency analysis to separate astrophysical 
foreground signals from the CMB, each of which improves on our first-year analy- 
ses. First, we form an improved "Internal Linear Combination" (LLC) map, based 
solely on WMAP data, by adding a bias correction step and by quantifying resid- 
ual uncertainties in the resulting map. Second, we fit and subtract new spatial 
templates that trace Galactic emission; in particular, we now use low-frequency 
WMAP data to trace synchrotron emission instead of the 408 MHz sky survey. 
The WMAP point source catalog is updated to include 115 new sources whose 
detection is made possible by the improved sky map sensitivity. 

We derive the angular power spectrum of the temperature anisotropy using 
a hybrid approach that combines a maximum likelihood estimate at low I (large 
angular scales) with a quadratic cross-power estimate for / > 30. The resulting 
multi-frequency spectra are analyzed for residual point source contamination. 
At 94 GHz the unmasked sources contribute 128 ± 27 /iK 2 to 1(1 + l)Ci/2n 
at I = 1000. After subtracting this contribution, our best estimate of the CMB 
power spectrum is derived by averaging cross-power spectra from 153 statistically 
independent channel pairs. The combined spectrum is cosmic variance limited to 
/ = 400, and the signal-to-noise ratio per /-mode exceeds unity up to I = 850. For 
bins of width Al/l = 3%, the signal-to-noise ratio exceeds unity up to I = 1000. 
The first two acoustic peaks are seen at I = 220.8 ± 0.7 and / = 530.9 ± 3.8, 
respectively, while the first two troughs are seen at I = 412 .4 ±1.9 and I = 
675.2 ± 11.1, respectively. The rise to the third peak is unambiguous; when 
the WMAP data are combined with higher resolution CMB measurements, the 
exi stence of a thi r d aco ustic peak is well established. 



Spergel et ajj (120061 ) use the three-year temperature and polarization data 



to constrain cosmological model parameters. A simple six parameter ACDM 
model continues to fit CMB data and other measures of lar ge scale structure 
remarkably well. The new polarization data (IPage et al.ll2006l ) produce a better 
measurement of the optical depth to re-ionization, r = 0.089±0.03. This new and 
tighter constraint on r helps break a degeneracy with the scalar spectral index 
which is now found to be n s = 0.958 ±0.016. If additional cosmological data sets 
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are included in the analysis, the spectral index is found to be n s = 0.947 ±0.015. 

Subject headings: cosmic microwave background, cosmology: observations, early 
universe, dark matter, space vehicles, space vehicles: instruments, instrumenta- 
tion: detectors, telescopes 



INTRODUCTION 



The Wilkinson Microwave Anisotropy Probe (WMAP) is a Medium-class Explorer 
(MIDEX) mission designed to elucidate cosmology by producing full-sky maps of the cosmic 
microwave background (CMB) anisotropy. Results from the first year of WMAP observations 
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I. The data were made available to the research community via the Legacy Archive for 



Microwave Background Data Analysis (LAMBDA), NASA's CMB T hematic Data Cent er 



and were described in detail in the WMAP Explanatory Supplement (ILimon et al.l 120031 ) . 



Papers based on the first-year WMAP results cover a wide range of topics, including: 
constraints on inflation, the nature of the dark energy, the dark matter density, implications 
for supersymmetry, the CMB and WMAP as the premier baryometer, intriguing features 
in the large-scale data, the topology of the universe, deviations from Gaussian statistics, 
time-variable cosmic parameters, the Galactic interstellar medium, microwave point sources, 
the Sunyaev-Zeldovich effect, and the ionization history of the universe. The WMAP data 
has also been used to establish the calibration of other CMB data sets. 

Our analysis of the first three years of WMAP data is now complete and the results are 



prese nted here and in companion papers (IJarosik et al.ll2006l : iPage et al.ll2006l : ISpergel et al. 



20061 ). The three- year WMAP results improve upon the first-year set in many ways, the 
most important of which are the following. (1) A thorough analysis of the polarization data 
has produced full-sky polarization maps and power spectra, and an improved understanding 
of many aspects of the data. (2) Additional data reduces the instrument noise, producing 
power spectra that are 3 times more sensitive in the noise limited regime. (3) Independent 
years of data enable cross-checks that were not previously possible. (4) The instrument 
calibration and beam response have been better characterized. 



This paper presents the analysis of the three-year temperature data, focusing on fore- 
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ground modeling and removal, evaluation of the angular power spectrum, and selected top- 
ics beyond the power spectrum. Com panion papers pr esent the new polarization maps and 
polarization-specific scientific results (IPage et al.l 120061). a n d discuss the cosmol ogical impli- 
cations of the three- year WMAP data fjSpergel et al. 12006 ). Jarosik et al. ( 2006 ) present our 
new data processing methods and place systematic error limits on the maps. 

In £J2] we summarize the major changes we have made to the data processing since the 
first-year analysis, and §3] presents a synopsis of the three-year temperature maps. In §0]we 
discuss Galactic foreground emission and our attempts to separate the emission components 
using a Maximum Entropy Method (MEM) analysis. $5] illustrates two methods we employ 
to remove Galactic foreground emission from the maps in preparation for CMB analysis. $6] 
updates the WMAP point source catalog and presents a search for the Sunyaev-Zeldovich 
effect in the three- year maps. §7] evaluates the angular power spectrum and compares it to 
the previous WMAP spectrum and to other contemporary CMB results. In £JH] we survey 
the claims that have been made regarding odd features in the WMAP first-year sky maps, 
and we offer conclusions in 39l 



2. CHANGES IN THE THREE- YEAR DATA ANALYSIS 



The first-year data analysis was described in detail in the suite of first-year WMAP 
papers listed above. In large part, the three- year analysis employs the same methods, with 
the following exceptions. 

In the first-year analysis we subtracted the COBE dipole from the time-ordered data to 
minimize the effect of signal aliasing that arises from pixelizing a signal with a steep gradient. 
Since the WMAP gain calibration procedure uses the Doppler effect induced by WMAP's 
velocity with respect to the Sun to establish the absolute calibration scale, WMAP data 
may be used independent ly to determine the C MB dipole. Consequently, we subtract the 
WMAP first-year dipole (IBennett et al.ll2003bl ) from the time-ordered data in the present 
analysis. 

A small temperature dependent pointing error (~ 1 arcmin) was found during the course 
of the first-year analysis. The effect is caused by thermal stresses on the spacecraft structure 
that induce slight movement of the star tracker with respect to the instrument. While the 
error was small enough to ignore in the fir st-year data, it is no w corrected with a temperature 
dependent model of the relative motion (I Jarosik et al.l 120061 ). 



The radiometer gain model described by iJarosik et al.l (l2003bl ) has been updated to 
include a dependence on the temperature of the warm-stage (RXB) amplifiers. While this 
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term was not required by the first-year data, it is required now for the model to fit the full 
three-year d ata with a single par ameterization. The new model, and its residual errors, are 
discussed by IJarosik et al.l (120061 ) . 



The WMAP beam response has now been measured with six independent "seasons" 
of Jupiter observations. In addition, we have now produced a physical model of one side 
of our symmetric optical sys tem, the A-side, b ased on simultaneous fits to all 10 A-side 
beam pattern measurements (IJarosik et al.l 120061 ) . We use this model to augment the beam 
response data at very low signal-to-noise ratio, which in turn allows us to determine better 
the total solid angle and window function of each beam. 

The far sidelobe response of the beam was determined from a combination of gr ound 
measurements and in-flight lunar data taken early in the mission (IBarnes et al.l 120031 ) . In 
the first-year processing we applied a small far-sidelobe correction to the K-band sky map. 
For the current analysis, we have implemented a n ew far sidelobe cor rection and gain re- 
calibration that operates on the time-ordered data (IJarosik et al.ll2006l ). These corrections 
have now been applied to data from all 10 differencing assemblies. 

When producing polarization maps, we account for differences in the freq uency pass- 



band between the two linear polarization channels in a differencing assembly (IPage et al. 



20061 ). If this difference is not accounted for, Galactic foreground signals would alias into 
linear polarization signals. 

Due to a combination of 1/f noise and observing strategy, the noise in the WMAP 
sky maps is correlated from pixel to pixel. This results in certain low-/ modes on the sky 
being less well measured than others. This effect can be completely ignored for temperature 
analysis since the low-/ signal-to-noise ratio is so high, and the effect is not important at 
high-/ ( §7.1.2p . However, it is very important for polarization analysis because the signal-to- 
noise ratio is so much lower. In order to handle this complexity, the map-making procedure 
has been overhauled to produce genuine maximum likelihood solutions that employ optimal 
filtering of the time-or dered data and a c onjugate- gradient algorithm to solve the linear 
map-making equations ( Jarosik et al. 20061 ). In conjunction with this we have written code 



to evaluate the full pixel-to-pixel weight (inverse covariance) matrix at low pixel resolution. 
(The H EALPix convention is to denote pixel resolution by the parameter A^de, with A" P i x = 
12A r s 2 ide (IGorski et al.ll2005l ). We define a resolution parameter r such that A^de = 2 r . The 
weight matrices have been evaluated at resolution r4, A^e = 16, Ap ix = 3072 .) The full 



noise covariance information is propagated through the power spectrum analysis (IPage et al. 



20061 ). 



When performing template-based Galactic foreground subtraction, we now use tern- 
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plates based on WM AP K- and Ka-band data in place of the 408 MHz synchrotron map 
(IHaslam et al.lll98ll ). As discussed in §5.3[ this substitution reduces errors caused by spec- 
tral index variations that change the spatial morphology of the synchrotron emission as a 
function of frequency. A simi lar model is used for subtracting polarized synchrotron emission 
from the polarization maps (IPage et al.ll2006l ). 



We have performed an error analysis of the internal linear combination (ILC) map and 
have now implemented a bias correction as part of the algorithm. We believe the map is now 
suitable for use in low-/ CMB signal characterization, though we have not performed a full 
battery of non-Gaussian tests on this map, so we must still advise users to exercise caution. 
Accordingly, we present full-sky multipole moments for I — 2, 3, derived from the three-year 
ILC map. 

We have improved the final temperature power spectrum (Cf T ) by using a maximum 
likelihood estimate for low-/ and a pseudo-C/ estimate for I > 30 (see ^7j). The pseudo-C/ 
estimate is simplified by using only V- and W-band data, and by reducing the number of 
pixel weighting schemes to two, "uniform" and "iV obs " ( §7.5p . With three individual years of 
data and six V- and W-band differencing assemblies (DAs) to choose from, we can now form 
individual cross-power spectra from 15 DA pairs within a year and from 36 DA pairs across 
3 year pairs, for a total of 153 independent cross-power spectra. In the first-year spectrum 
we included Q-band data, which gave us 8 DAs and 28 independent cross-power spectra. 
The arguments for dropping Q-band from the three-year spectrum are given in §7.21 

We have developed methods for estimating the polarization power spectra {Cf x for 
XX = TE, TB, EE, EB, BB) from temperature and polarization maps. The main technical 
hurdle we had to overcome in the pro cess was the prop er handling of low signal-to- noise ratio 
data with complex noise properties (IPage et al.ll2006l ). This step, in conjunction with the 
development of the new map-making process, was by far the most time consuming aspect of 
the three-year analysis. 

We have improved the form of the likelihoo d function used to infer cosmological pa- 
rameters from the Monte Carlo Markov Chains (jSpergel et al.ll2006l ). In addition to using 
an exact maximum likelihood form for the low-/ TT data, we have developed a method to 
self-consistently evaluate the joint likelihood of temperatur e and polarization data given a 
theoretical model (described in Appendix D of IPage et al.l (120061 )). We also now account 
for Sunyaev-Zeldovich (SZ) fluctuations when estimating parameters. Within the WMAP 
frequency range, it is difficult to distin guish between a pr i mord ial CMB spectrum and a 
thermal SZ spectrum, so we adopt the iKomatsu fc Seljaki (120021 ) model for the SZ power 
spectrum and marginalize over the amplitude as a nuisance parameter. 
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We now use the CAMB code (ILewis et al.ll2000l ) to compute angular power spectra from 
cosmological parameters. CAMB is derived from CMBFAST (ISeljak fc Zaldarriagalll996l ). 
but it runs faster on our Silicon Graphics (SGI) computers. 



OBSERVATIONS AND MAPS 



The three-year WMAP data encompass the period from 00:00:00 UT, 10 August 2001 
(day number 222) to 00:00:00 UT, 9 August 2004 (day number 222). The observing efficiency 
during this time is roughly 99%; Table [2] lists the fraction of data that was lost or flagged as 
suspect. The Table also gives the fraction of data that is flagged due to potential contami- 
nation by thermal emission from Mars, Jupiter, Saturn, Uranus, and N eptune. These dat a 
are not used in map- making, but are useful for in-flight beam mapping (ILimon et al.l 120061 ) . 



Sky map s are created from the time-ordered data using the procedure described by 
Jarosik et al.l (120061 ). For several reasons, we produce single- year maps for each year of 
the three-year observing period (after performing an end-to-end analysis of the instrument 
calibration). We produce three-year maps by averaging the annual maps. Figure [1] shows 
the three-year maps at each of the five WMAP observing frequencies: 23, 33, 41, 61, and 
94 GHz. The number of independent observations per pixel, A Q b s , is displayed in Figure [2J 
The noise per pixel, p, is given by a{p) = <7oN ohs ' (p), where a is the noise per observation, 
given in Table HI To a very good approximation, the noise per pixel in the three-year maps 
is a factor of v3 times lower t han in the one-year maps. The noise properties of the data 
are discussed in more detail in Ijarosik et ali ( 2006 ). 



The three- year maps are compared to the previously released maps in Figure [3j Both set 
of maps have been smoothed to 1° resolution to minimize the noise difference between them. 
When viewed side by side they look indistinguishable. The right column of Figure [3] shows 
the difference of the maps at each frequency on a scale of ±30 fiK. Aside from the noise 
reduction and a few bright variable quasars, such as 3C279, the main difference between 
the maps is in the large-scale (low-/) emission. This is largely due to improvements in our 
model of the instrument gain as a function of time, which is ma de possible by having a 
longer time span with which to fit the model (I Jarosik et al.ll2006l ). In the specific case of 
K-band, the improved far-sidelobe pickup correction produced an effective change in the 
absolute calibration scale by ~1%. This, in turn, is respons ible for the difference seen in the 
bright Galactic plane signal in K-band (I Jarosik et al.ll2006l ). We discuss the low-/ emission 
in detail in §7.41 and §3 but we stress here that the changes shown in Figure [3] are small, 
even compared to the low quadrupole moment seen in the first-year maps. Table [3] gives 
the amplitude of the dipole, quadrupole, and octupole moments in these difference maps. 
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For comparison, we estimate the CMB power at I = 2, 3 to be AT 2 = 236 and 1053 yuK 2 , 
respectively ( §7.4p . 

As discussed in §SJ several authors have noted unusual features in the large-scale signal 
recorded in the first-year maps. We have not attempted to reproduce the analyses presented 
in those papers, but based on the small fractional difference in the large-scale signal, we 
anticipate that most of the previously reported results will persist when the three-year maps 
are analyzed. 



4. GALACTIC FOREGROUND ANALYSIS 

The CMB signal in the WMAP sky maps is contaminated by microwave emission from 
the Milky Way Galaxy and from extragalactic sources. In order to use the maps reliably for 
cosmological studies, the foreground signals must be understood and removed from the maps. 
In this section we present an overview of the mechanisms that produce significant diffuse 
microwave emission in the Milky Way and we assess what can be learned about them using 
a Maximum Entropy Method (MEM) analysis of the WMAP data. We discuss foreground 
removal in £j5l 



4.1. Free-Free Emission 



Free-free emission arises from electron-ion scattering which produces microwaves with 
a brightness spectrum Ta ~ (EM/1 cm~ 6 pc) u~ 2AA for frequencies v > 10 GHz, where EM 
is the emission measure, / njdl, and we assume an electron gas temperatu re T e ~ 8000 K. 



1998; Haffneretal 



2003; 



As discussed in bennett et al l (l2003ch. high-resolution maps of Ha em ission ( IDennison et al.l 



Reynolds et al.l 120021 ; iGaustad et al.l l200ll ) can serve as approxi 



mate tracers of free-free emission. The intensity of Ha emission is given by 



I(R) =0.44£(r d ) (EM/1 cnr 6 pc) (T e /8000K) 



-0.5 



0.34 ln(T e /8000 K)] , 



where I is in Rayleighs (1 R = 2.42 x 10 -7 ergs cm _2 s _1 sr _1 at the Ha wavelength of 0.6563 
(Urn), the helium contribution is assumed to be small, and is an extinction factor that 
depends on the dust optical depth, r d , at the wavelength of Ha. If the emitting gas is co- 
extensive with dust, then C( T d) = [1 — exp(— T d )]/r d . Ha is in R-band, where th e extinction 



i s 0.75 times visible, A R = 0.75 Ay; thus, A R = 2.35E B _ V , and r d = 2.2E B _ V . iFinkbeiner 
(120031 ) assembled a full-sky Ha map using data from several surveys: the Wisconsin H-Alpha 
Mapper (WHAM), the Virginia Tech Spectral-Line Survey (VTSS), and the Southern H- 



Alpha Sky Survey Atlas (SHASSA). We use this map, together with the lSchlegel et al.l (119981 ) 
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(SFD) extinction map, to predict a map of free- free emission in regions where r r j<l, under the 



assump tion that the dust and ionized gas are co-extensive. As discussed in iBennett et al. 



( 12003d ). this template has known sources of uncertainty and error. We use it as a prior 
estimate in the MEM analysis ( §4.51) . and as a free- free estimate in the template-based 
foreground removal ( §5.31) . 



4.2. Synchrotron Emission 



Synchrotron emission arises from the acceleration of cosmic ray electrons in magnetic 



fields. In our Galaxy, discrete supernova remnants contribute only ~10% of t 



chrotron emission at 1.5 GHz flLisenfeld fc Vo lk 2000; lBiermannlll976l ; Ul vest ad 



re tot al syn- 
1982 f ) while 



~90% of the observed emission arises from a diffuse component. iHummel et al.l (119911 ) 
present maps of synchrotron emission at 610 MHz and 1.49 GHz from the edge-on spiral 
galaxy NGC891. They find that the synchrotron spectral index varies from (3 S ~ —2.6 in 
most of the galactic plane to (3 S ~ —3.1 in the halo. Similar spectral index variations are 
seen in the Milky Way at ~1 GHz, where the synchrotron signal is complex. Variations 
of the synchrotron spectral index are both expected and observed. Moreover, the emission 
is dominated at low frequencies by components with steep spectra, whereas at higher fre- 
quencies it is dominated by components with flatter spectra, usually with a different spatial 
distribution. As a result, great care m ust be taken when using low frequency maps, like the 
408 MHz map of lHaslam et al.l (Il98ll ). as tracers of the synchrotron emission at microwave 
frequencies. 

Synchrotron emission can be highly polarized. Theoretically, the linear polarization 
fraction ca n be a s high as ~75%, though values < 30% are more typically observed. See 
Page et al.l (120061 ) for a discussion of the new full-sky observations of polarized synchrotron 
emission in the three-year WMAP data. 



4.3. Thermal Dust Emission 



Thermal dust emission has been mapped over th e full sky in several infrared bands by 
the IRAS and COBE missions. ISchlegel et al.l (119981 ) combined data from both missions to 



produ ce an absolutely calibrated full-sky map of the thermal dust emission. iFinkbeiner et al. 



( 119991 ) (FDS) extended this work to far-infrared and microwave frequencies using the COBE- 
FIRAS and COBE-DMR data to constrain the low-frequency dust spectrum. They fit the 
data to a particular two-component model that gives power-law emissivity indices a\ = 
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1.67 and oli = 2.70, and temperatures of T\ = 9.4 K and T 2 = 16.2 K. The fraction of 
power emitted by each component is fx = 0.0363 and / 2 = 0.9637, and the relative ratio 
of IR thermal emission to optical opacity of the two components is gi/ffe = 13.0. The cold 
component is potentially identified as emission from amorphous silicate grains while the 
warm component is plausibly carbon based. Independent of the physical interpretation of 
the model, FDS found tha t it fit t he data moderately well, with xl — 1-85 for 118 degrees 
of freedom, bennett et all J2003ah noted that this model, call "Model 8" by FDS, did well 
predicting the first-year WMAP dust emission. 

It is reasonable to assume that the Milky Way is like other spiral galaxies and that the 
microwave properties of external galaxies should help to inform our understanding of the 
global properties of the Milky Way. It has long been known that a remarkably tight correla- 



ternal galaxies. This relation has been extensively studied and modeled ( 
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All theories attempting to explain this tight correlation are tied 



to the level of the star formation activity. During this cycle, stars form, heat, and destroy 
dust grains; create magnetic fields and relativistic electrons; and create the O- and B-stars 
that ionize the surrounding gas. However, it is not clear what these models predict on a 
"microscopic" (cloud by cloud) level within a galaxy. 



Bennett et al.l ((2003d) showed that the synchrotron and dust emission in our own Galaxy 
are spatially correlated at WMAP frequencies. Many authors have argued that this correla- 
tion is actually due to radio emission from dust grains themselves, rather than from a tight 
dust-synchrotron correlation. We review the evidence for this more fully in the next section. 



4.4. "Anomalous" Microwave Emission from Dust? 



With the advent of high-quality diffuse microwave emission maps in the early 1990s, it 
became possible to study the high-frequency ta il of the sync h rotron spectrum and the low- 
frequency tail of the interstellar dust spectrum. Kogut et al. ( 1996al lbl) analyzed foreground 
emission in the CQBE^DM R maps and report ed a signal that was significantly correlated 
with 240 11m dust em ission (lArendt et al.lll998l ) but not with 408 MHz synchrotron emission 
( IHaslam et al.lll98ll ). The correlated signal was notably brighter at 31 GHz than at 53 GHz 
{(3 ~ —2.2), hence they concluded it was consistent with free-free emission that was spatially 
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correlated with dust. The same conclusion was reached by Ide Oliveira-Costa et al.l (119971 ). 
who found the Saskatoon 40 GHz data to be correlated with infrared dust, but not with 
radio synchrotron emission. 



Leitch et al.l (119971 ); iLeitchl (119981 ). and iLeitch et al.l (120001 ) analyzed data from the 
"RING5m" experiment. A likelihood fit to their 14.5 GHz and 31.7 GHz data, assum- 
ing CMB anisotropy and a single foreground component, produced a foreground spectral 
index of (3 — —2.581^42. The data would have preferred a steeper value had it not been 
for an assumed prior limit of (3 > —3. This signal was fully consistent with synchrotron 
emission. However, a puzzle arose i n comparing the RIN G5m data with a the Wester- 
bork Northern Sky Survey (WENSS) faengelink et alill997h at 325 MHz: the WENSS data 
showed no detectable signal in the vicinity of the RING5m field. The 325 MHz limit implies 
that (3 > —2.1 and rules out conventional synchrotron emission as the dominant foreground. 
(Since the WENSS is an interferometric survey primarily designed to study discrete sources, 
the data are insensitive to zero-point flux from extended emission. It is not clear how much 
this affects the above conclusion.) As with the DMR and Saskatoon data, the 14.5 GHz 
foreground emission was correlated with dust, but it was difficult to attribute it to spatially 
correlated free-free emission because there was negligible Ha emission in the vicinity. To 
reconcile this, a gas temperature in excess of a million degrees would be needed to sup- 
press the Ha. Flat spectrum synchrotron was also suggested as a possible source; it had 
been previously observed in other sky regions and it would obviate the need for such a high 
temperature and pressure. 



Draine fc Lazarianl (119981 ) dismissed the hot ionized gas explanation on energetic grounds 
and instead suggested that the emission (which they described as "anomalous" ) be attributed 
to electri c dipole rotation al emission from very small dust grains - a mechanism first pro- 
posed by lEricksonl (119571 ) in a different astrophysical context. One of the hallmarks of this 
mechanism is that it produces a frequency spectrum that peaks in the 10-60 GHz range and 
falls off fairly steeply on either side. 



de Oliveira-Costa et al.l (119981 ) analyzed the nearly full-sky 19 GHz sky map (IBoughn et al 



19921 ) and found some correlation with the 408 MHz synchrotron emission, but found a 
stronger correlation with the COBE^DIRBE 240 fim dust emission. They concluded the 19 
GHz data were consistent with either free-free or spinning dust emission. 



Leitchl (119981 ) commented that the preferred model of iDraine fc Lazarianl (119981 ) could 



produce the RING5m foreground component at 31.7 GHz, but that it only accounted for 
at most 30% of the 1 4.5 GHz emission, even when adopting unlikely values of the grain 
dipole moment. Since ILeitch et al.l (120001 ) were primarily interested in studying the CMB 
anisotropy, they considered using the IRAS 100 fim map as a foreground template to remove 
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the "anomalous" emission, regardless of its physical origin. They found, however, that fitting 
only a CMB component and a dust-correlated component produced an unacceptably high 
xl = 10 per degree of freedom. Thus, while the radio foreground morphology correlates with 
dust, the correlation is not perfect. 



Finkbeiner et al.l (120021 ) used the Green Bank 140 foot telescope to search for spinning 



dust emission in a set of dusty sources selected to be promising for detection. Ten infrared- 
selected dust clouds were observed at 5, 8, and 10 GHz. Eight of the ten sources yielded 
negative results, one was marginal, and one (the only diffuse H II region of the ten, LPH 
201.6) was claimed as a tentative detection based on its spectral index of > —2. Recogniz- 



i ng th at this spectrum does not necessarily imply spinning dust emission, iFinkbeiner et al. 



( 120021 ) offer three additional requirements to convincingly demonstrate the detection of spin- 
ning dust, and concluded that none of the three requirements was met by the existing data. 
The absence of a rising spectrum in most of the sources may be taken as evidence that 
spinning dust emission is not typically dominant in this spectral region, at least for this 
type of infrared-selected cloud. The tentative detection in LPH 201.6 has met with three 
criticisms: (1) lack of eviden ce for the premise t hat its radio emission is proportional to 



its far-infrared dust emission (ICasassus et al.ll20Q4J) , (2) the p utative spinning dust emission 
is stronger than theory predicts (jMcCullough fc Chenll2002l ). and (3) the positive spectral 



index may be accounted for by unresolved optically thick emission (jMcCullough &: Chen 
20021 ). On the latter point, follow-up observations failed to identify a compact HII region 
candidate (McCullough, private communication). 



Bennett et al.l (j2003d ) fit the first-year WMAP foreground data to within ~1% using a 



Maximum Entropy Method (MEM) analysis (see also §4.51) . As with the above-cited results, 
WMAP found that the 22 GHz to 33 GHz foregrounds are domin ated by a component 
with a synchrotron-like spectrum, but a dust-like spatial morphology. iBennett et al.l (12003a ) 
suggested that this may be due to spatially varying synchrotron spectral indices acting over 
a large frequency range, significantly altering the synchrotron morphology with frequency. 
A spinning dust component with a thermal dust morphology and the Draine and Lazarian 
spectrum could not account for more than ~ 5% of the emission at 33 GHz. Of course, the 
WMAP fit did not rule out spinning dust as a sub-dominant emission source (as it surely 
must be at some level), nor did it rule out spinning dust models with other spectra or spatial 
morphologies. 



Casassus et al.l (12004 ) report evidence for anomalous microwave emission in the Helix 



planetary nebula at 31 GHz, where at least 20% of the emission is correlated with 100 /zm 
dust emission. They rejected several explanations. The observed features are not seen in 
H/?, ruling out free-free emission as the source. Cold grains are also ruled out as the source 
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by the absence of 250 GHz continuum emission. Very small grains are not expected to 
survive in planetary nebulae, a nd none have b e en de tected in the Helix, but Fe is strongly 
depleted in the gas. Instead, ICasassus et al.l (120041 ) favor the notion of magnetic dipole 
emissio n (produced by variation s in grain magnetization) from hot ferromagnetic classical 
grains (IDraine fe Lazarianl 119991 ). Although the derived emissivity per n ucleon in the Helix 



i s a fa ctor of ~ 5 larger than the highest end of the range predicted by 



Draine fe Lazarian 



( 1999 ). this excess could be explained by a high dust temperature, since IDraine fe Lazarian 
( 119991 ) assume an ISM temperature of 18 K instead of a typical planetary nebula dust 
temperature of ~ 100 K. The fraction of 31 GHz Helix emission attributable to free-free 
is estimated to be in the range 36 — 80%. This low level of free- free emission implies an 
electron temperature of T e = 4600 ± 1200 K, which is much lower than the value Td ~ 9000 
K based on collisionally exci ted lines. This discrep ancy may be du e to strong temperature 
variations within the nebula. ICasassus et al.l (120041 ) suggest that the lFinkbeiner et al.l (120021 ) 
measurement of LPH 201.6 may also be produced by magnetic dipole emission from classical 
dust grains. 



de Oliveira-Costa et al.l (120041 ) correlate the Tenerife 10 and 15 GHz data (iGutierrez et al 
20001 ) with the WMAP non-thermal ( "synchrotron" ) map that was produced as part of the 



first-year Maximum Entropy Method (MEM) analysis of the WMAP foreground signal. They 
detect a low frequency roll-off in the correlated emission, as shown in Figure [6j We consider 
this result after discussing the three-year MEM analysis in the next section. 



Fernandez-Cerezo et al.l (120061 ) report new measurements with the COSMOSOMAS ex- 
periment covering 9000 square degrees with ~1° resolution at frequencies of 12.7, 14.7, and 
16.3 GHz. In addition to CMB signal and what is interpreted to be a population of un- 
resolved radio sources, they find evidence for diffuse emission that is correlated with the 
DIRBE 100 /im and 240 /xm bands. As with many of the above-mentioned results, they 
find that the correlated signal amplitude rises from 22 GHz (using WMAP data) to 16.3 
GHz, and that it shows signs of flattening below 16.3 GHz "compatible with predictions for 
anomalous microwave emission related to spinning dust." 

The topic of anomalous dust emission remains unsettled, and is likely to remain so 
until high-quality diffuse measurements are available over a modest fraction of sky in the 
5-15 GHz frequency range. We offer some further comments after presenting the three-year 
MEM results in the following section. 
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4.5. Maximum Entropy Method (MEM) Foreground Analysis 



Bennett et al.l (12003d ) described a MEM-based approach to modeling the multifrequency 
WMAP sky maps on a pixel-by-pixel basis. Since the method is non-linear, it produces maps 
with complicated noise properties that are difficult to propagate in cosmological analyses. 
As a result, it is not a promising method for foreground removal. However, the method is 
quite useful in helping to separate Galactic foreground components by emission mechanism, 
which in turn informs our understanding of the foregrounds. 

We model the temperature map at each frequency, u, and pixel, p, as 

T m (v,p) = T cmb (p) + S s (v,p) T s (p) + S s (v,p)T s (p) + Sd(u,p)T d (jp), (2) 

where the subscripts cmb, s, ff, and d denote the CMB, synchrotron (including any anomalous 
dust component), free-free, and thermal dust components, respectively. T c (p) is the spatial 
distribution of emission component c in pixel p, and S c (v,p) is the spectrum of emission c, 
which is not assumed to be uniform across the sky. We normalize the spectra (S = 1) at K- 
band for the synchrotron and free-free components, and at W-band for the dust component. 



The model is fit in each pixel by minimizing the functional H = A + XB f lPress et al. 



19921 ). where A = ~Yl v \T{i>,p) — T m (u,p)] 2 /al, is the standard x 2 °f the model fit, and B = 
'Yuf.Tdjp) ln[T c (p) / P c (p)] is the MEM functional (see below). The parameter A controls the 
relative weight of A (the data) and B (the prior information) in the fit. In the functional B, 
the sum over c is restricted to Galactic emission components, and P c (p) is a prior estimate of 
T c (p). The form of B ensures the positivity of the solution T c (p) for the Galactic components, 
which greatly alleviates degeneracy between the foreground components. 

Throughout the MEM analysis, we smooth all maps to a uniform 1° angular resolution. 
To improve our ability to constrain and understand the foreground components, we first 
subtract a prior estimate of the CMB signal from the data rather than fit for it. We use 
the ILC map described in §5.21 for this purpose and subtract it from all 5 frequency band 
maps. Since WMAP employs differential receivers, the zero level of each temperature map is 
unspecified. For the MEM analysis we adopt the following convention. In the limit that the 
Galactic emission is described by a plane-parallel slab, we have T(|6|) = T p esc where T p 
is the temperature at the Galactic pole. For each of the 5 frequency band maps, we assign the 
zero level such that a fit of the form T(|6|) = T p esc \b\ + c, over the range —90° < b < —15°, 
yields c = 0. 



We construct a prior estimate for dust emission, Pd{p), using Model 8 of lFinkbeiner et al. 



(119991 ) . evaluated at 94 GHz. The dust spectrum is modeled as a straight power law, Sd(v) 



(z//Vyi/) +2 ' - For free-free emission, we estimate the prior, P$(p), using the extinction- 
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corrected Ha map (lFinkbeinerll2003l ) . This is converted to a free-free signal using a conversion 
factor of 11.4 /iK R -1 (units of antenna temperature at K-band). We model the spectrum 
as a straight power law, Ss(u) = {y /vk)~ 2 ' 1a ■ As noted in §4.1| the main source of uncer- 
tainty in this free-free estimate is the level of extinction correction (in addition to any Ha 
photometry errors) . We reduce A in regions of high dust optical depth to minimize the effect 
of errors in the prior. 

For the synchrotron emission, we construct a prior estimate, Pjp), by subtractin g an 



extragalactic brightness of 5.9 K from the Haslam 408 MHz map (lLawson et al.l 119871 ) and 
scaling the result to K-band assuming (3 S = —2.9. Since the synchrotron spectrum varies 
with position on the sky, this prior estimate is expected to be imperfect. We account for this 
in choosing A, as described below. We construct an initial spectral model for the synchrotron, 
S B (u,p), using the spectral index map @ s (p) = /5(408 MHz, 23 GHz). Specifically, we form 
S B (u,p) = ( z/ / z/ ^)A-o.25.[/3 s +3.5] for Ka-band, and S a (v,p) = S B (u Ka ,p)-(u/u Ka )' 3s ' - 7 ^ +3 ^ for 
Q-, V-, and W-band (S s = 1 for K-band). This allows for a /3-dependent steepening of the 
synchrotron spectrum at microwave frequencies. In our first-year results, use of this initial 
spectral model produced solutions with zero synchrotron signal, T s (p), in a few low latitude 
pixels for which the K-Ka spectrum is flatter than free-free emission. For the three-year 
analysis, this problem is handled by setting the initial model, S s , to be flatter than free-free 
emission in these pixels. 

For all three emission components, the priors P c and the spectra S c are fixed during each 
minimization of H. As described below, we iteratively improve the synchrotron spectrum 
model based on the residuals of the fit. 

The parameter A controls the degree to which the solution follows the prior. In regions 
where the signal is strong, the data alone should constrain the model without the need for 
prior information, so A can be small (though we maintain A > to naturally impose positivity 
on T c and reduce degeneracy among the emission components). As in the first-year analysis, 
we base \(p) on the foreground signal strength: X(p) = 0.6 [Tk(p)/1 mK]~ 15 , where Tk(p) is 
the K-band map (with the ILC map subtracted) in units of mK, antenna temperature. This 
gives 0.2<A<3. 

After each minimization of H, we update S s for K- through V-band according to 

Sr(u,p) T s (p) = S s (u,p) T B (p) + 9 ■ R(u,p) (3) 

where R(u,p) = T(u,p) — T m (u,p) is the solution residual, and g is a gain factor that we set 
to 0.5. For W-band we update S s by power-law extrapolation using the inferred Q,V-band 
spectral index. For some low signal-to- noise ratio pixels, we occasionally find (3 s {vq, vy) > 0. 
In a change from our first-year analysis, we now extrapolate the synchrotron spectrum from 
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V- to W-band using j3 s = —3.3. This change affects only 7% of pixels. 

We iterate the minimization of H and the update of S s eleven times. At the end of this 
cycle, the residual, R(v,p), is less than 1% of the total signal T(v,p). However, there are 
still several sources of potential error in the component decomposition, including: 

(1) Zero-level uncertainty. As noted above, we use a plane-parallel slab model to assign 
the sky map zero point at each frequency. We estimate the uncertainty in this convention 
by fitting the model separately in the northern and southern Galactic hemispheres. These 
separate fits change the zero levels by as much as 4 fiK. When these differences are propa- 
gated into T(u,p), the output maps, T c , change by <15%, <5%, and <2% for the free-free, 
synchrotron, and dust components, respectively, at low Galactic latitudes. 

(2) Dust spectrum uncertainty. Changing (3& by 0.2 changes the component maps by 
<10%, <3%, and <2% for the free-free, synchrotron, and dust components, respectively, at 
low Galactic latitudes. 

(3) Dipole subtraction uncertainty. A 0.5% dipole error would result in a systematic 
0.5% gain error in all bands. The dominant effect would be to rescale each component map 
by 0.5%. 

(4) CMB signal subtraction uncertainty. Errors in the ILC estimate of the CMB signal 
will produce errors in the corrected Galactic data, T(v,p) — T ILC (p), used in the MEM 
analysis. We quantify ILC errors in §5.2} but note here that they are small compared 
to the total Galactic signal at low latitudes. However, they plausibly dominate the total 
uncertainty at high Galactic latitudes, where the WMAP data are (fortunately) dominated 
by CMB signal. Fractional uncertainties of ~100% are not unlikely in the foreground model 
at high latitudes. 

(5) While the total model residuals are small, there is still potentially significant un- 
certainty in the individual foreground components. The program outlined above produces 
three component maps, T s , Tg, and Tj, and the synchrotron spectrum model S s (v,p). To 
illustrate the degeneracy between these outputs, we consider a simplified single-pixel model 
of the form 

T m (u) = (S s (u,p))T s + S s (u)T s + S d (v)T d , (4) 

where the angle brackets indicate a full-sky average, and we explicitly evaluate the average 
MEM functional, 

H = J2 1(T(M) ; Tm0/)]2 + (A(P)> £T C ln[T c /<P c (p)>], (5) 

V u c 

for selected pairs of parameters, while marginalizing over the rest. Contour plots of H are 
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shown in Figure HI For the panels that explore the shape of the synchrotron spectrum, S s , 
we parameterize it as a power law with a steepening parameter, (3 B (v) = @ + d(3 B /d\ogv ■ 
[logy — log^R-] and evaluate if as a function of df3 s /d log u, while marginalizing over j3 . 
For the rest, we iteratively update S B as per equation fl3j). For the most part, the output 
parameters are only weakly correlated; the most notable degeneracy is between the free-free 
amplitude and the synchrotron amplitude. WMAP data tightly constrain the sum of the two, 
but their difference is determined by the relative amplitude of the prior estimates for these 
two components, and by the initial synchrotron spectrum model. The synchrotron spectrum 
is found with modest significance to be steepening with increasing frequency. However, 
the dust index, (3& (not shown), is poorly constrained; thus our conclusion in the first-year 
analysis that (j3d) = 2.2 was not well founded upon further investigation. 

Figure shows the three input prior maps, P c {p), and the corresponding output compo- 
nent maps, T c (p), obtained from the three-year data. These maps are available on LAMBDA 
as part of the three-year data release. The maps are displayed using a logarithmic color 
stretch to highlight a range of intensity levels. The morphology and amplitude of the ther- 
mal dust emission are well predicted by the prior (FDS) dust map (see also §5.31) . The 
free-free emission is generally over-predicted by the prior discussed above, especially in re- 
gions of high extinction. But even in regions where the extinction is low, we find the mean 
free-free to Ha ratio at K-band to be closer to ~8 /zK R -1 than the value of 11.4 /iK R _1 
assumed in generating the prior. Moreover, we find considerable variation in this ratio (a 
factor of ~2) depending on location. 

The most notable discrepancy between prior and output maps is seen in the synchrotron 
emission. Specifically, the K-band signal has a much more extended Galactic longitude distri- 
bution than does the 408 MHz emission, and it is remarkably well correlated with the thermal 
dust emission. Is this K-band non-thermal component due to anomalous dust emission or to 
mostly flat-spectrum synchrotron emission that dominates at microwave frequencies and is 
well correlated with dusty active star-forming regions? We cannot answer this question with 
WMAP data alone because the frequency range of WMAP does not extend low enough to 
see the predicted rollover in the low-frequency anomalous dust spectrum. However, we note 
the following points. 

In Figure [6] we show the mean spectra of the three Galactic emission components ob- 
served by WMAP, in addition to the sum of the three. For comparison, we also show the sig- 
nal observed at 408 MHz and we infer the signals at 10 and 15 GH z ba sed on correlation analy- 



ses of the Tenerife CMB data by lde Oliveira-Costa et al.l (j2004d ) and lde Oliveira-Costa et al. 



( 119991 ). The curves show the mean signal in the range 20° < |6| < 50°, as computed from 



the output MEM component maps: blue is dust, green is free-free, red is the non-thermal 
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signal, and black is the sum of the three. The total intensity of the 408 MHz emission is 
remarkably well matched to a simple power-law extrapolation of the total WMAP signal 
measured from K-band to V-band. The spectral index of the dashed black curve is —2.65 
between 408 MHz and K-band, and —2.69 between 408 MHz and Q-band. If we interpret 
the non-thermal emission as synchrotron, the implied spectral index between 408 MHz and 
the K-band non-thermal component is —2.73. 



As noted in §4.4[ Ide Oliveira-Costa et al.l (120041 ) correlate the Tenerife 10 and 15 GHz 
data with the first-year WMAP non-thermal MEM map. We infer the non-thermal signal at 
10 and 15 GHz (shown in red) by scaling the three-year non-thermal map using the reported 



correl ation coefficients. We find a frequency rollover consistent with Ide Oliveira-Costa et al. 



( 12004] ) . Using the same scaling method, we also infer the 10 and 15 GHz signals derived from 
correlations with Ha emission (green arrows) and the Haslam 408 MHz (synchrotron) map 
(grey points). The derived free-free emission is lower than the extrapolation of the free- free 
emission inferred from WMAP - which is not physically tenable. Similarly, the Haslam- 
correlated emission at 15 GHz is substantially lower than it is at either 10 or 23 GHz. Thus, 
the sum of all correlated components in the 10 and 15 GHz data requires a substantial dip 
in the spectrum of the total emission (solid black points), which is not hinted at by the 
Haslam or WMAP data (dashed black points). This may be a sign of substantial anomalous 
emission, but one must be cautious interpreting the spectra of correlation coefficients. 



Page et al.l (120061 ) have used the three-year WMAP polarization data to construct a 
novel decomposition of the intensity foreground signal at each WMAP frequency band. In 
brief, they predict a synchrotron polarization fraction from a model of the Galaxy's mag- 
netic field strength and electron density. This fraction is de-rated by an empirical factor 
to account for missing structure in the model, then a "high-polarization-fraction" compo- 
nent of the intensity signal is formed as T high (z/, p) = P(u,p)/ f(p), where P(u,p) is the 
polarization intensity at frequency v, in mK, and f(p) is the model polarization fraction. 
After subtracting an estimate of the CMB and free-free signal from the intensity maps, 
the remaining non-thermal signal is attributed to a "low-polarization fraction" component, 
Ti ow (z/, p) = T'{y,p) — T high (z/, p), where T' is the temperature map corrected for CMB and 
free-free emission. T he high- frac t ion co mponent has a morphology similar to the 408 MHz 
emission (Figure 9 in Page et al. J 20061)'). w hile low- fraction component has a very dust-like 
morphology (Figure 13 in lPage et al.l (120061 ) ) . While the accuracy of this decomposition de- 
pends on the model fraction, f(p), the basic picture should inform future studies of anomalous 



emission. 
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5. GALACTIC FOREGROUND REMOVAL 

The primary goal of foreground removal is to provide a clean map of the CMB for 
cosmological analysis; an improved understanding of foreground astrophysics is a secondary 
goal. Removal techniques typically rely on the fact that the foreground signals have quite 
different spatial and spectral distributions than the CMB. In this section we describe two- 
and-a-half approaches to foreground removal that use complementary information. The first 
is an update of the Internal Linea r Combination (ILC) method we employed in the first- 



year analysis (IBennett et al.ll2003d ). The second is an updated approach to fitting Galactic 
emission templates to each WMAP frequency band map. The remaining strategy is to mask 
regions of the sky that are too contaminated to be reliably cleaned. Extragalactic sources are 
treated in £j6j For our primary CMB results, we analyze the masked, template-subtracted 
maps, but for some low-/ applications we also analyze the ILC map as a consistency check. 



5.1. Temperature Masks 



Many regions of the sky are so strongly contaminated by foreground signals that reliable 
cleaning cannot be assured. These regions are masked for cosmological analy sis, though the 



extent of the masking required depends on the type of analysis being done. IBennett et al. 



((2003d) defined a set of pixel masks based on the first-year K-band temperature map. Since 
these masks were based on high signal-to-noise ratio Galactic signal contours in the K-band 
data, we have not modified the diffuse emission masks for the three-year analysis. 

In addition to diffuse Galactic emission, point sources also contaminate the WMAP data. 
A po int source ma s k was constructed for the first-year analysis that included all of the source s 



from lStickel et all (11994 ); s ources with 22 Gl 



flat spectrum o bjects from 
Perlman et all (119981 ) and 



Terasranta et al. 



Landt et all (I200lh . 



z fluxe s > 0.5 Jy from lHirabayashi et all (120001 ); 
J200ll ); and sources from the blazar survey of 



The mask contained nearly 700 objects, 
including all 208 of the sources directly detected by WMAP in the first-year data. Each 
source was masked to a radius of 0.6°. For the three-year analysis, we have supplemented 
the source mask with objects from the three-year WMAP source catalog discussed in §6.11 Of 
the newly detected sources, 81 were not included in the previous mask and have been added 
to the three-year mask. Weaker, undetected sources still contribute to the high-/ angular 
power spectrum. As discussed in §7.2[ we fit and subtract a residual source contribution to 



the multi-frequency power spectrum data. 

Figure [7] gives an overview of the microwave sky and indicates the extent of the various 
foreground masks. The yellow, salmon, and red shaded bands indicate the diffuse masks 
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defined in iBennett et al.l (12003d). Th e violet shading shows the "P06" polarization analysis 
mask described in iPage et al.l (120061 ). The small blue dots indicate point sources detected 
by WMAP (to alleviate crowding, the full source mask described above is not shown). In 
addition, some well-known sources and regions are specifically called out. 



5.2. The Internal Linear Combination (ILC) Method 

Linear combinations of the multi-frequency WMAP sky maps can be formed using 
coefficients that approximately cancel Galactic signals while preserving the CMB signal. This 
approach exploits the fact that the frequency spectrum of foreground emission is different 
from that of the CMB. The method is "internal" in that it relies only on WMAP data, so 
the calibration and systematic errors of other experiments do not enter. There are a number 
of ways the coefficients can be determined, some of which require only minimal assumptions 
about the nature of the foreground signals. In the first year WMAP papers we introduced a 
method in which the coefficients were determined by minimizing the variance of the resulting 
map subject to the constraint that the coefficients sum to unity, in order to preserve the 
CMB signal. We called the resulting map the "ILC" map. In this section we elaborate on 
the strengths and limitations of the ILC method and quantify the uncertainties in the ILC 
map. 



Eriksen et al.l (120041 ) have also analyzed the method as an approach to foreground re- 
moval. They devised an approach to variance minimization that employed a Lagrange mul- 
tiplier to linearize the problem and dubbed the resulting map the "LILC" map, where the 
first L denotes Lagrange. They found their LILC map differed somewhat from the ILC map 
in certain regions of the sky. We have since verified that the two minimization methods 
produce identical results for a given set of inputs and that the differences were due to an 
ambiguity in the way the regions were defined in the original ILC description. Because the 
linearized algorithm is considerably faster than our original nonlinear minimization, we have 
adopted it in the present work. 



5.2.1. Uniform Foreground Spectra 

In order to better understand how errors arise in the ILC map, we first consider a 
simple scenario in which instrument noise is negligible and the spectrum of the foreground 
emission is uniform across the sky, or within a defined region of the sky. In this case, 
a frequency map, Tj(p) = T(z/j,p), may be written as a superposition of a CMB term, 
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T c (p), and a foreground term, SiTf(p), where Si = S'(i'i) describes the composite frequency 
spectrum of the foreground emission, and T f (jo) describes the spatial distribution, so that 
Tj(p) = T c (p) + SiTf(p). A linear combination map has the form 

Ti L c(p) = ^(P) = E Ci ( T c(p) + £T f (p)] = T c (p) + r T f (p), (6) 

j i 

where we have imposed the constraint — 1? an d have defined T = X^Ci'S'i- 

Suppose we choose to determine the coefficients Q by minimizing the variance of Tilc- 
Then, 

= <^ C (P)> " (TilM) 2 (7) 
= <T C 2 > - (T c > 2 + 2r [(T c T f ) - (T c > (T f )] + L 2 [<T f 2 > - (T f ) 2 ] (8) 
= a 2 + 2IV cf + T V f 2 (9) 

where the angle brackets indicate an average over pixels, and we have defined the variance 
and covariance in terms of these averages. Note that this expression would still hold if we 
added an arbitrary constant to each frequency map, Tj — > Tj + T 0ji . The ILC variance will 
be minimized when 

- -w = V rf + 2r ^' (10) 

Thus the coefficients C« that minimize a 2 LC give T = — a ci /af, and in the absence of noise, 
the corresponding ILC solution is 

Tilc(p) = T c (p) - a* /of T t (p), (11) 

with 

<&c = °l ~ (12) 
In this ideal case, the frequency maps combine in such a way as to maximize the cancellation 
between CMB signal and foreground signal, producing a biased CMB map with <r 2 LC < o\. 
We have tested this result with ideal simulations in which we generate 5 frequency maps, 
Tj, which include a Galaxy signal with a constant spectrum, Si, and random realizations of 
CMB signal and instrument noise. We then generate ILC maps from each realization and 
compare the residual map, Tilc — T c , to the bias prediction, — a c f/afT{. The results confirm 
that the above description is correct, and that instrument noise is not a significant concern 
in this situation. The level of the bias is typically ~10 /iK in the Galactic plane. 



5.2.2. Non-uniform Foreground Spectra 



To minimize the anti-correlation bias we should choose regions that minimize the co- 
variance between the CMB and the foreground, (T c T f ). However, in the previous analysis 
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we assumed that the spectra of the foreground signals were constant over the sky. In reality 
these will vary as the ratio of synchrotron, free-free, and dust emission varies across the 
sky (and as the intrinsic synchrotron and dust spectra vary). In this case, the bias analysis 
becomes more complex. Specifically, the foreground component at each frequency may be 
written as Si(p)Tf(p), and the ILC map takes the form 

T 1LC (p) = T c (p) + T(p)T i (p), (13) 

where T(p) = YliCiS>i(p)- The ILC variance then generalizes to 

^ilc = (T?) - (T c f + 2 [<T c rT f > - (T c > (TT f )] + [<r 2 T f 2 > - (rT f ) 2 ] . (14) 

Using the same reasoning that led to equation (flOj) . we obtain the following result for the 
minimum variance solution 

(rr f ■ SiT f ) = - (T c • SiT f ) . (15) 

This has the same interpretation as equation (flOl) . in the sense that it relates the foreground 
variance to the CMB-foreground covariance. We can solve this equation for T(p) by noting 
that r(p) = J2iCiSi(p), so that 

&T { ■ SjTt) Q = - (T c ■ S^) . (16) 

3 

Now define Fy = (SiTf ■ SjTf) and Cj = (T c ■ SjTf), whereby 

T = J2QS l = -J^S l -(F-%-C J , (17) 

i ij 

which is the multi-frequency analog of equation (fill . Once again though, the bias in the 
ILC solution is proportional to (minus) the CMB-foreground covariance. 

We have tested this expression with simulations like the ones described above, except 
this time we employ a three-component Galaxy model with variable spectra, S(vi,p), based 
on the first-year MEM model. As we discuss in more detail below, the simulations verify that 
the output ILC map is biased, T IL c(p) — T c (p) = r(p)T f (p), with T as given in equation (fT7|) . 
Unfortunately, we do not know If and T c a priori, and it has proven difficult to relate 
this bias expression to the frequency band maps, Tj, in a way the can be used to minimize 
the bias. As a result, we have primarily resorted to correcting the bias with Monte Carlo 
simulations, as we describe below. 

Given the results above, and our previous experience with the ILC method, it is clear 
that one should subdivide the sky into regions selected by foreground spectra, in order to 
reduce bias prior to correcting it. We have carried out such a program for the three-year 
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analysis and have found it very difficu lt to improve on the region selection made in the 



first-year analysis (IBennett et al.ll2003d ). Nonetheless, we have adopted a few changes: (1' 



we eliminat e the Taurus A reg ion, as it is too small to ensure a reliable CMB-foreground 



separation (lEriksen et al.ll2004j ); (2) we add a new region to minimize dust residuals in the 
Galactic plane. This region is based on a Ty — Tw color selection and encompasses the outer 
Galactic plane within the Kp2 cut. The new ILC region map is shown in Figure [HJ The 
region designated 1, shown in red, replaces the old Taurus A region; the remaining regions 
are unchanged. This map is available on LAMBDA as part of the three-year data release. 

For each region n, we determine a set of band weights, £n,ij by minimizing the variance of 
the linear combination map T n (p) = Y^i=i Cn,%T%{jp) in that region, subject to the constraint 
'Ylii Cn,i — 1- There are two exceptions to note. The coefficients for region were derived from 
a subset of the data in that region, specifically pixels inside the Kp2 cut with \l\ > 60°. The 
coefficients for region 1 were derived from a slightly larger region of data, specifically pixels 
inside the Kp2 cut with |/| > 50°, that pass the Ty — Tw color cut. To ensure uniformity, 
the band maps have been smoothed to a common resolution of 1° FWHM. The coefficients 
£ n ,i are given in Table 

We form a full-sky map by combining the A" region maps, T n ; but to minimize edge 
effects, we blend the region maps as follows. We create a set of A" full-sky weight maps, w n 
such that w n {p) = 1 for p £ R n and w n {p) = otherwise. We smooth these maps (which 
contain only ones and zeros) with a 1?5 smoothing kernel, to get smoothed weight maps, 
w n (p). The final full-sky map is then given by 

Tip) = S "jf" (P) f r " (p) , (18) 

Ln W n{P) 

where the sum is over the twelve sky regions. 

In order to obtain the final bias correction, we generate multi-region Monte Carlo sim- 
ulations, using the variable spectrum, MEM-based Galaxy model as input. We evaluate the 
error, Tilc — T c , for each realization and compute a bias from the mean error averaged over 
100 realizations. The result, shown in Figure [HI is roughly 20-30 fiK in the Galactic plane, 
but substantially less off the plane. This map was used to correct the three-year ILC map, 
which is shown in the middle panel of Figure [9j This Figure also shows the first-year ILC 
map (top panel) for comparison. The difference between the two (bottom panel) is primarily 
due to the new bias correction, but a small quadrupole difference, due to the changes noted 
in Figure [3], is also visible. 



Based on the Monte Carlo simulations carried out for this ILC study, we estimate that 
residual Galactic removal errors in the three-year ILC map are less than 5 fiK on angular 
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scales greater than ~ 10°. But we caution that on smaller scales, there is significant structure 
in the bias correction map that is still uncertain. On larger scales, we believe the three-year 
ILC map provides a reliable estimate of the CMB signal, with negligible instrument noise, 
over the full sky. We analyze the low-/ multipole moments of this map in §7.31 



5.3. Foreground Template Subtraction 

The ILC method discussed above produces a CMB map with complicated noise prop- 
erties, while the MEM method discussed in §4.51 is primarily used to identify and separate 
foreground components from each other. For most cosmological analyses one must retain 
the well-defined noise properties of the WMAP frequency band maps. To achieve this we 
form low-noise model templates of each foreground emission component and fit them to the 
WMAP sky maps at each frequency. After subtracting the best-fit model we mask regions 
that cannot be reliably cleaned because of limitations in the template models. In this section 
we describe the model templates we use for synchrotron, free-free, and dust emission and we 
estimate the residual foreground uncertainties that remain after these templates have been 
fit and subtracted. The WMAP band maps are calibrated in thermodynamic temperature 
units; where appropriate, we convert Galactic signals to units of antenna temperature using 
the factors g v given in Table [TJ 

In our first-year model we used the Haslam 408 MHz map as a template for synchrotron 
emission. We now use the WMAP K- and Ka-band data to provide a synchrotron template, 
as described below. This is preferable because: (1) the intrinsic systematic measurement 
errors are smaller in the WMAP data than in the Haslam data, and (2) the non-uniform 
synchrotro n spectrum produces morphological changes in the brightness as a function of 



frequency (IBennett et al.ll2003d ). so that the low frequency Haslam map is less reliable at 



tracing microwave synchrotron emission than the WMAP data. 

There are two potential pitfalls associated with using the K- and Ka-band data for clean- 
ing: (1) the data are somewhat noisy, and since the template subtraction will be common 
to all cleaned channels, there can be a noise bias introduced in the inferred angular power 
spectrum. (But note that we use separate templates for each year of data, so the correlation 
only acts across frequency bands within a single year.) (2) Since the K- and Ka- band data 
are contaminated with point sources, this signal could interfere with the primary goal of 
cleaning the diffuse emission. Using the fitting coefficients obtained below and the known 
noise properties of the K- and Ka-band data, we estimate the noise bias in the final power 
spectrum to be <5 /zK 2 near the 1st acoustic peak (<0.1% of the CMB signal), and even 
smaller at lower and higher multipoles. Further, assuming the point source model given in 
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equation (|43p . and the fact that the template has been smoothed to an effective resolution of 
1?0 FWHM, we estimate that sources contribute <1 /1K 2 to the power spectrum at I = 400 
in the T K — T Ka template, and thus may be safely ignored. In the end, these pitfalls are not 
a source of concern for the three-year analysis. 

The difference map Tk — Tx a , in thermodynamic units, cancels CMB signal while it 
contains a specific linear combination of synchrotron and free-free emission (and a minimal 
level of thermal dust emission). We use this map as the first te mplate in the mo del. For 
the second template we use t he full-sky Ha map c ompiled by iFinkbeinerl (120031 ) with a 
correction for dust extinction (IBennett et al.l 12003a ). This template independently traces 
free-free emission, allowing the model to produce an arbitrary ratio of synchrotron to free- 
free emission at a given frequency (the limitations of Ha as a proxy for free-fr e e are discussed 
below). For dust emission, we adopt "Model 8" from the iFinkbeiner et al.l (119991 ) analysis 
of IRAS and COBE data, evaluated at 94 GHz (see §4.31) . The full model has the form 



M(v,p) = &iO)(T K - T Ka ) + b 2 (v)I Ha + b 3 (u)M d , 



(19) 



where bi(u) are the fit coefficients for each template at frequency z/, and M d is the dust 
map. As discussed below, this model is simultaneously fit to the Q-, V-, and W-band maps, 
and we constrain the coefficients 62 and 63 to follow specified frequency spectra to minimize 
component degeneracy. 

To clarify the physical interpretation of b\ and 62 we first note that Tk — Tx a may be 
rewritten in terms of synchrotron and free-free emission as 



Tk — Tka — Rs T s + Rg Tff , 



(20) 



where T s and Tg are the synchrotron and free-free maps in antenna temperature at ri- 
band, R c = gKSdvKiV) — 5'Ka>S'c( z/ Ka,p) is the surviving fraction of emission component c 
(synchrotron or free-free) in Tk — TK a , and S c is the spectrum of component c, in antenna 
temperature, relative to K-band. To a very good approximation, the spectrum of free-free 
emission is S s = (z//z/ K )" 214 (SZD, so that R s = °- 552 - For synchrotron emission, variations 
in the spectrum as a function of position will produce variations in R s . For spectral indices 
in the range (3 S = —2.9 ± 0.2 we have R s = 0.667 ± 0.026. In the remainder of this section 
we assume R s = 0.67 and neglect the ~2.5% error introduced by spectral index variations 
between K and Ka bands. Note that this value of R s is only used to estimate the level of 
synchrotron emission in the template Tk — TK a ; we do not constrain the fit coefficients b\ to 
follow a specified frequency spectrum. 



By adding the Ha map to the model, we allow the synchrotron to free-free ratio to vary 
as a function of frequency, but we must be cognizant of potential errors introduced by the use 
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of Ha as a proxy for the free-free emission. Nominally we have = hsln a , where /hq is the 
Ha intensity in Rayleighs and is th e free-free to Ha ratio. At K-band kg is predicted to be 



~11.4 /iK R -1 OBennett et al.ll2003d ) but the actual ratio is both uncertain and dependent 
locally on extinction and reflection. These effects make the Ha proxy unacceptable in the 
Galactic plane, and force one to mask these regions for CMB analysis. Outside the masked 
regions the variations in kg are primarily due to residual extinction and to variations in the 
temperature of the emitting gas. Here higher fractional errors can be tolerated because the 
total free-free signal is fainter. 

Due to the uncertainties in the free-free to Ha ratio, and the fact that Tk — Xk a contains 
a mixture of synchrotron and free-free emission, care must be taken to interpret the model 
correctly. Let the combined synchrotron and free-free emission in the data at frequency v be 

T{v,p) = g{v) [S s (v)T s (p) + S s (v)T s (p)\ , (21) 

where the terms are as defined above. The synchrotron and free-free terms in the model may 
be written as 

& 1 0)(T K -T Ka ) + b 2 {v)I Ha = 6 1 (z/)[i? s T s + J R ff T ff ] + [6 2 (z/)//i ff ]T ff (22) 

= [h(v)R*\ T s + [b^Rs + b 2 (v)/h s ] T s . (23) 

Comparing the synchrotron terms in equations (1231 and (j2ip we can infer the mean syn- 
chrotron spectral index returned by the fit 

\og[S s (u)} hglb^Rjgju)] 

Ps{.VK, V) - — ; r - -—, r • (24) 

log(u/v K ) log(v/u K ) 

Comparing the free-free terms, and assuming Sg is known, we can solve for the free-free to 
Ha ratio at K-band 

ns gMSiM-hMR*' 1 ] 

We fit the template model, equation ( |T9|) . simultaneously to each of the eight Q- through 
W-band differencing assembly (DA) maps (the three-year maps smoothed to 1° resolution) 
by minimizing x 2 

^2 _ ^ P^'P) ~ b M( T K - T Ka) - b 2 {v)lHa ~ h(v)M d } 2 ^ 
i.p * 

where T(vi,p) is the WMAP sky map from DA i (in thermodynamic units), of is the mean 
noise variance per pixel for DA i, and the second sum is over pixels outside the Kp2 sky 
cut. To regularize the model, we impose the following constraints on the fit coefficients: 
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(1) all coefficients must be positive-definite, (2) the dust coefficients must follow a spectrum 
b${vi) = 63 ■ g(z/j)(z/j/z/\y 1 ) +2 '°, and (3) the free-free coefficients for each DA must follow a 
free-free spectrum, which leads to the following form 

b 2 (vi) = b 2 ■ [g(vi)(vi/vK)- 2 - U - hMRs] ■ (27) 

The synchrotron coefficients are fit separately for each differencing assembly Given the 10 
coefficients from the three- year fit, we subtract the model from each single-year DA map to 
produce a set of cleaned maps. In doing so, we form separate single-year maps of Tk — Tk a 
to maintain rigorously independent noise between separate years of data. 

The fit coefficients 6j are given in Table H] along with derived values for f3 s and hg. 
To facilitate model subtraction, we tabulate values for b 2 and 63 for each DA using the 
above constraints. Note that the FDS dust model, which predates WMAP by a few years, 
predicts the 94 GHz dust signal remarkably well. The synchrotron e mission shows a stead y 



steepening with increasing frequency, as seen in the first-year data (jBennett et al.ll2003d j. 
Also, the free-free to Ha ratio is seen to be ~6.5 /zK R _1 , which is roughly half of the 11.4 
fiK R _1 prediction. Taken together, this fit finds a remarkably low total Galactic foreground 
amplitude at V-band. 

Figure [TU] shows the three- year band maps before and after subtracting the above model. 
For comparison, the Figure also shows the same th ree-year maps after subtracting the first- 



year template-based model (jBennett et al.l 12003d ). In all panels an estimate of the CMB 



signal (the ILC map) has been subtracted to better show residual foreground errors. The 
main visible difference between the first-year and three-year residual maps is the synchrotron 
subtraction error in the first-year model due to the use of the Haslam 408 MHz map. This is 
especially visible in the region of the North Polar Spur and around the inner Galaxy. Note 
also the significant model errors visible inside the Kp2 sky cut. This is presumably caused by 
a combination of synchrotron spectral index variations and errors in the extinction correction 
applied to the Ha template. Indeed, errors of up to 30 fiK are also clearly visible in isolated 
regions outside the cut, especially in the vicinity of the Gum Nebula and the Ophiuchus 
complex. In §7.21 we compare the foreground signal to the CMB and assess the degree to 
which these residual errors contaminate the CMB power spectrum. 
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EXTRAGALACTIC FOREGROUNDS 



6.1. Point Sources 



Extragalactic point sources contaminate the WMAP anisotropy data and a few hundred 
of them are strong enough that they should be masked and discarded prior to undertaking 
any CMB analysis. In this section we describe a new direct search for sources in the three- 
year WMAP band maps. Based on this search, we update the source mask that was used in 
the first-year analys is. In jj7.2l we des cribe our approach to fitting and subtracting residual 
sources in the data. iPage et al.l (120061 ) discuss the treatment of polarized sources. 



For the first-year analysis, we constru cted a catalog of sou rces surveyed at 4.85 GHz us- 
ing the northe rn hemisphere GB6 catalo g ( Gregory et allll996l) an d the southern hemisphere 



PMN catalog ([Griffith et al.lll994 Il995l : IWright et al 



1994 



19961 ). The GB6 catalog covers 



the declination range 0°< 5 < + 75° to a flux limit of 18 mJy, while the PMN catalog covers 
— 87°<5< + 10° to a flux limit between 20 and 72 mJy. Combined, these catalogs contain 
119,619 sources, with 93,799 in the region |6|>10°. We have examined the three-year WMAP 
sky maps for evidence of these sources as follows: we bin the catalog by source brightness 
and, for each bin, we cull the corresponding sky map pixels that contain those sources. The 
data show a clear correlation between source strength and mean sky map temperature that 
disappears if the sky map pixels are randomized. The multi-frequency WMAP data suggest 
that the detected sources are primarily flat-spectrum, with a ~ 0. 



In the first-year analysis, we produced a catalog of bright point sources in the WMAP sky 

maps, independent of their presence in external surveys. This process has been repeated with 

1 si 

the three-year maps as follows. We filter the weighted maps, N / T (N b R is the number of ob- 
servations per pixel) in harm onic space by bi/(bfCi mb +C^ olse ), (ITegmark fc de Oliveira-Costa 



998 



Refregier et all 120001). where bi is the transfer function of the WMAP beam response 
(IPage et al.ll2003al : Ijarosik et alJl2006h . Cf mh is the CMB angular oower spectrum, and C" olse 



is the noise power. Peaks that are >5<r in the filtered maps are fit in the unfiltered maps to 
a Gaussian profile plus a planar baseline. The Gaussian ampl itude is converted t o a source 
flux density using the conversion factors given in Table 5 of IPage et al.l (j2003al ). When a 
source is identified with >5cr confidence in any band, the flux densities for other bands are 
given if they are >2cr and the fit source width is within a factor of 2 of the true bea m width. 
We cross-correlate detected sources with the GB6, PMN, and iKuhr et al.l (Il98ll ) catalogs 
to identify 5 GHz counterparts. If a 5 GHz source is within 11' of the WMAP source posi- 
tion (the WMAP source position uncertainty is 4') we tag the WMAP source. When more 
than one source lies within the cutoff radius the brightest one is assumed to be the WMAP 
counterpart. 
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The catalog of 323 sources obtained from the three-year maps is listed in Table |9j In 
the first-year catalog, source ID numbers were assigned on the basis of position (sorted by 
galactic longitude). Now, rather than assigning new numbers to the newly detected sources, 
we recommend that WMAP sources be referred to by their coordinates, e.g., WMAP J0006- 
0622. For reference, we give the first-year source ID in column 3 of Table [91 The 5 GHz 
IDs are given in the last column. The source count distribution, dN/ dS, obtained from the 
Q-band data is shown in Figure [TTJ 

The first-year catalog contained 208 sources. Given the increased sensitivity in the 
three-year maps, the number of new sources detected is consistent with expectations based 
on differential source count models. By the same token, 6 sources from the first-year catalog 
were dropped from the three-year list (numbers 15, 31, 61, 96, 156, and 168). Simulations 
of the first-year catalog suggested that it contained 5 ± 4 false detections, so the number 
of dropped sources is consistent with expectations. Five of the 208 sources in the first-year 
catalog could not be identified with 5 GHz counterparts; now 6 out of 323 can not be. Of 
the 6 sources that dropped off the first-year catalog, source numbers 15 and 61 did not have 
5 GHz identifications. The remaining 4 may have dropped off because of variable (declining) 
flux density. 



Trushkinl (120031 ) has compiled multifrequency radio spectra and high resolution radio 



maps of the sources in the first-year WMAP catalog. Reliable identifi cations are claim ed for 



205 of the 208 sources. Of the 203 sources with optical identifications, iTrushkinl (120031 ) finds 
141 quasars, 29 galaxies, 19 active galactic nuclei, 19 BL Lac-type objects and one planetary 
nebula, IC418. Forty percent of the sources are identified as having flat and inverted radio 
spectra, 13% have GHz-peaked spectra, 8% are classical power-law sources, and 7% have a 
classical low frequency p ower-law combined with a flat or inverted spectrum component (like 



3C84). ITrushkinl (120031 ) suggests that WMAP source number 116 is likely to be spurious 



and, for source 61 no radio component was found. 



6.2. Sunyaev-Zeldovich (SZ) Effect 

Hot gas in galaxy clusters produces secondary anisotropy in the CMB via the Sunyaev- 
Zeldovich effect: a systematic frequency shift of CMB photons produced by Compton scat- 
tering off hot electrons in the cluster gas. At frequencies less than 217 GHz, this pro- 
duces a temperature decrement in the CMB along sight lines that pass through clusters. 
The effect is relatively small for WMAP due to its moderate angular resolution; neverthe- 
less, nearby clusters are large enough to produce a measurable effect in the sky maps. In 
our first year analysis we reported detections of an SZ signal from the Coma cluster and 
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from an aggregate sample of X-ray selected Abell clusters, the XBAC catalog. Since then, 



data ( 


Fosalba & Gaztanaeall2004; 


Fosalba et al.l 


2003; 


Mvers et al.ll2004 


;lAfshordi et al. 


2004 




Hernandez- Monteaeudo e 


t al. 


2004; 


Hernandez-Monteaeudo & Rubino-Martmll2004 


Afshordi et al. 


20051 


Hansen et al. 


20051 


Atrio-Barandela & Muecke 


J 2006h. In this section, we update our 



results for the Coma cluster and the XBAC catalog; in subsequent cosmological studies, we 
mask these clusters from the data. 

The brightest SZ source is the Coma cluster. For SZ analysis we model its temperature 
profile with an isothermal /9-model of the form 

AT SZ (6) = ATsz(0) [1 + (6/6 c ) 2 } (1 ' 3/3)/2 , (28) 

where 9 = r jDA z) is the angular distance from the core, located at (l,b) = (56?75, 88?05); 
Briel et al.l Jl992h give 6 C = 10.5' ± 0.6' and (3 = 0.75 ± 0.03. Usin g these values, we 



convolve the profil e with the WMAP beam response at each frequency (jPage et al.l 12003a 
Jarosik et al1l20od ) and fit for ATsz(0) by minimizing 



X 



[Ti - T - AT sz (0i)] (M-% [Tj - T - AT SZ (^-)] 



(29) 



where T is the temperature in pixel i, T is a temperature baseline, 0i is the angular distance 
of pixel % from the cluster core, and is the pixel-pixel covariance matrix of the sky map. 



We model Mj,- as C(n 



+ crfSij, where C(9) is the 2-point correlation function of the 



CMB and of is the variance due to instrument noise in pixel i. 

For W-band we find a Coma decrement of ATsz(0) = — 0.46±0.16 mK, with a reduced y 



of 1.05 for 196 degrees of freedom. In V-band the decrement is — 0.31±0.16 mK. lHerbig et al. 



(119951 ) measured a Rayleigh- Jeans decrement of —0.51 ± 0.09 mK at 32 GHz. Given the SZ 
frequency spectrum, their result predicts — 0.47± 0.08 mK at V-band and —0.42 ± 0.07 mK 
at W-band, consistent with what we observe. 



The XBAC catalog of X-ray clusters produced by lEbeling et al.l ( 119961 ) is an essentially 
complete flux-limited sample of 242 Abell clusters selected from the ROSAT All-Sky Survey 
data. We cross-correlate this catalog with the WMAP 94 GHz map. Treating the XBAC 
clusters as poin t sourc es, we estimate their flux density at 94 GHz using the expression from 
Refreeier et all (hoooh 

kT e \ 5/2 
lkeV/ 



s< 



91 



11.44 



300 Mpc 



D z) 



/gas 

or! 



[mJy], 



(30) 



where D(z) is the angular diameter distance to the cluster, / gas = M gas /M is the gas mass 
fraction, and T e is the electron temperature. The overall normalization of this relation is 



31 



uncertain due to our ignorance of the correct gas fraction and cluster virialization state. We 
fix /gas = 0.11. Extended clusters (more than one-third the extent of Coma) are omitted; the 
remaining fluxes are all < 1 Jy. A templat e map is construc ted by convolving the clusters 



with the WMAP W-band beam response (IPage et al.ll2003al ) and fitting the result to the 



WMAP 94 GHz map. We find a template normalization of —0.32 ± 0.14 for the three-year 
data (compared to —0.36 ± 0.14 in the first year data). Since the fluxes used to construct 
the template are positive, the negative scaling is consistent with a 2.5a SZ decrement. 

CMB photons that travel to us through the plane of the Milky Way undergo an SZ 
distortion of y ~ kn e T e aTL / m e c 2 , where <jt is the Thomson scattering cross-section and L is 
the electron-pressure-weighted path length through our Galaxy. Taking n e T e = 10 3 K cm -3 
and L = 50 kpc we find i/«2x 10 -8 ; thus our Galaxy does not significantly affect CMB 
photons. The SZ effect can safely be ignored as a diffuse contaminating foreground signal. 



7. THE ANGULAR POWER SPECTRUM 

Full-sky maps provide the most compact record of CMB anisotropy without loss of 
information. They permit a wide variety of statistics to be computed from the data, the 
most fundamental of which is the CMB angular power spectrum. Indeed, if the tempera- 
ture fluctuations are Gaussian distributed, with random phase, the angular power spectrum 
provides a complete description of the statistical properties of the CMB. While there have 
been a number of papers based on the first-year data that claim evidence of non-Gaussianity 
and/or non-random phases in the fluctuations (see below) it is clear that the fluctuations 
are not strongly deviant from Gaussian, random phase. Thus, the measured power spectrum 
does provide the primary point of contact between data and cosmological parameters. This 
section presents the angular power spectrum obtained from the first three years of WMAP 
observations. 

A sky map T(n) defined over the full sky may be decomposed in spherical harmonics as 

oo I 

T(n) = a ^™( n ) (31) 

1=0 m=-l 

with 

a lm = J dnT(n)t(n), (32) 

where n is a unit direction vector. If the CMB anisotropy is Gaussian distributed with 
random phases, then each ai m is an independent Gaussian deviate with (tty m ) = 0, and 



\ a lm a l'm') ~ drum' Cj, 



(33) 
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where C\ is the angular power spectrum and S is the Kronecker symbol. C\ is the mean 
variance per I that would be observed by a hypothetical ensemble of observers distributed 
throughout the universe. The actual power spectrum observed in our sky, by a supposedly 
typical member of this ensemble, would be 



£<sky 



2Z + 1 ^ 



\ a lr 



(34) 



If we had noiseless CMB data over the full sky, equation (1321) could be evaluated exactly 
and equation (13"1"]) would give an unbiased estimate of the true power spectrum, in the sense 
that (C* ky ) = Q when averaged over the ensemble. However, we do not have ideal data (see 
below), and even if we did, since we only measure 21 + 1 modes per I (per sky), the above 
estimate of the variance has an intrinsic uncertainty (or "cosmic variance" ) of 



21 + 1 



(35) 



Beyond cosmic variance, there are two effects that preclude using equations (1321) and (l34l) to 
estimate the power spectrum; (1) real CMB data contains noise and other sources of error 
that cause the quadratic expression in equation (1341 to be biased; and (2) data near the 
Galactic plane must be masked because Galactic emission cannot be reliably cleaned there. 
Masking precludes the proper evaluation of the integral in equation (132|) . so other methods 
must be found to estimate ai m and C/. 

In the first-year analysis, we addressed problem (1) by adopting a "cross-power" es- 
timator in which we replace | a/ m .[ 2 by (a ^a^J, where i and j denote data channe l s wit h 



uncorrelated noise (IHinshaw et al.ll2003bl ). see also iPolenta et al.l (120051 ); iPatanchonl (120031 ); 
Tristram et al.l (12004 ) . This form removes the noise bia s from the estimate of |a/ m | 2 . We ad - 



dress ed problem (2) by using a "pseudo-CV' estimator (IPeebles fc Hauserlll974l ; iHivon et al. 



20021 ) that statistically corrects for the aliasing introduced by non-uniform pixel weights, of 
which the sky cut is an extreme case. 

An alternative approach is to invoke a maximum likelihood estimator that employs 
Bayes' theorem to relate the likelihood of the power spectrum to the likelihood of the data, 



L{d | d) oc L(d | Ci)P{Cy 



(36) 



Here L[C\ | d) is the likelihood of the model (the underlying power spectrum, C{) given the 
data d (usually a sky map, or its transform, the a; m coefficients), L(d | C{) is the likelihood 
of the data given the model, and P{C{) is the prior probability of the power spectrum. It is 
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standard to assume that the data follow a multivariate Gaussian distribution (as predicted 
by most inflationary models, for example) in which case the likelihood takes the form 



L(Q | d) oc 



exp(-|d T C- 1 d) 
V / detC 



■p(a: 



(37) 



where C is the covariance matrix of the data, including contributions from both signal and 
noise, C(Ci) = S(C/) + N. Since the covariance matrix includes both terms, the confidence 
regions on C\ deduced from L(Ci | d) automatically incorporate the uncertainty from both 
cosmic variance and instrument noise. In addition, if the data are restricted to lie outside 
a Galaxy cut, the confidence regions on Ci will incorporate the uncertainty due to aliasing. 
Finally, it is also standard to assume a uniform prior distribution for the model, P(Cj) = 1, 
though given sufficiently robust data, this choice will be unimportant. There are many 
technical challenges one must overcome when evaluating equation (13D) and most of these 



have been well covered in the literature (Peebles 



1997afi iTegmarkl Il997l : bond et all Il998 



Borrill 



973: Peebles fc Hauscr 1974 



19991 : lOh et all Il999h . See IPage et al. 



Hamilton 



20061 ) for an extension of this methodology to a simultaneous analysis of the temperature 
and polarization data. 

Since the WMAP first-year data release, several authors have re-derived the angular 
power spectrum from the sky maps using a variety of other estimators. Some of these anal- 
yses cover the full I range to which WMAP is sensitive, while others restrict attention to the 
low / regime where certain technical issues matter most, and where hints of unusual beh avior 
have been suggested. Of the analyses covering the full / range. iFosalba fc Szapudil (120041 ) em- 
ploy a pixel-space estimator called SpICE that produces a spectrum that agrees remarkably 
well with the released WMAP spectrum. These authors do find small discrepancies at the 
very highest I values (where the WMAP data is noise-dominated ) but these differences do 



not sig nificantly affect parameter determination. lO'Dwyer et al.l (120041 ) and lEriksen et al. 



( I2004bl ) use a novel method based on Gibbs Sampling to undertake a full Bayesian analysis 
of the power spectrum. This method also produces results that are in good agreement with 
the WMAP spectrum. In addition, the method is able to generate a full conditional likeli- 
hood for the spectrum that allows one to rigorously evaluate cosmological models without 
resorting to likelihood approximations. 



Efstathioul (j2004bl ) has surveyed a variety of methods for estimating the power spectrum 



of large data sets such as WMAP and concluded that the best approach is a hybrid one 
which employs a maximum likelihood estimate at low I and a pseudo-Q based estimate at 
high I. We agree with this approach and adopt it for the three-year analysis, making the 
transition at I = 30. This gives results that agree well with the first-year spectrum for I > 30, 
and reasonably well for low-/, as detailed below. In the original version of this paper, we 
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chose to make the t ransition at / = 12, but a detailed study of the three-year spectrum by 
Eriksen et al.l (120061 ) suggested that the transition needed to be at / ~ 30 to avoid biasing 



cosmo logic al parameter estima tes, particularly the scalar spectral index, n s . Eriksen et al. 
(bood ) and lSpergel et al.l (bood ) (Appendix A) discuss the effect of this choice on the spectral 
index. 

In the following sections we discuss the instrumental properties that are important to 
know for accurate power spectrum estimation, and assess the ability of the new Galactic 
foreground models to clean the power spectrum data. We then analyze the low and high-/ 
power spectrum in detail and compare the result with the previous WMAP spectrum and 
with other contemporary CMB data. 



7.1. Instrumental properties 



The temperature measured on the sky is modified by the properties of the instrument. 
The most important properties that affect the angular power spectrum are finite resolution 
and instrument noise. Let C l ^ %3 denote the auto or cross-power spectrum evaluated from 
two sky maps, ij and i'j', where ij is a composite DA/year index (i denotes the DA, j the 
year). Further, define the shorthand i = (ij,i'j') to denote a pair of composite indices, e.g., 
(Ql/yr-l,V2/yr-l), (Ql/yr-l,Ql/yr-2), etc. This spectrum will have the form 



Cl = wlC^ + Ni 



(3* 



where w\ = b\ bj' pf is the window function that describes the combined smoothing effects of 
the beam an d the finite sky ma p pixel size. Her e b) is the be a m tra nsfer function for DA 
i, defined by iPage et al.l (j2003al ) and updated by IJarosik et al.l (120061 ) [note that we reserve 
the term "beam window funct ion" for (6') 2 1, an d pi is the pixel transfer function supplied 



with the HEALPix package (IGorski et al.ll2005l ). is the noise spectrum realized in this 
particular measurement. On average, the observed spectrum estimates the underlying power 
spectrum, Ci, 

{C\)=w\C l + {Nt)8 ii/ , (39) 

where (N{) is the noise bias for differencing assembly i, and the Kronecker symbol indicates 
that the noise is uncorrelated between differencing assemblies and between years of data. 
To estimate the underlying power spectrum on the sky, Ci, the effects of the noise bias and 
beam convolution must be removed. The ability to determine these terms accurately is a 
critical element of any CMB experiment design. 



In §7.1.1l we summarize the results of IJarosik et al.l (120061 ) on the WMAP window func- 
tions and their uncertainties. We propagate these uncertainties through to the final Fisher 
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matrix for the angular power spectrum. In §7.1.21 we present a model of the WMAP noise 
properties appropriate to power spectrum evaluation. For cross power spectra (i 7^ %' above), 
the noise bias term drops out of the signal estimate if the noise between the two DAs (or 
between years of data from a single DA) is uncorrelated. However, the noise bias term still 
enters into the error estimate, even for cross spectra. Therefore, the noise model is used 
primarily in error propagation. 



7.1.1. Window Functions 



The instrument beam response was map ped in flight using observations of the planet 
Jupiter (IJarosik et al.ll2006l ; iPage et al.ll2003al ). The beam widths, measured in flight, range 
from 0?82 (FWHM) at K band down to 0?20 in some of the W band channels. The signal- 
to-noise ratio is such that the response, relative to the peak of the beam, is measured to 
approximately —35 dB in W band. As part of the three-year analysis, we have produced a 
physical model of the A-side optics based on simultaneous fits to all 10 A-side beam pattern 
measurements ( IJarosik et al.ll2006l ). We use this model to augment the beam response data 
at very low signal-to-noise ratio (—30 to —38 dB, depending on frequency band) which 
in turn allows us to better determine the total solid angle and window function of each 
A-side beam. For the B-side response we scale the A-side model by fitting it to the B- 
side Jupiter measurements in the high signal-to-noise regime. We then form similar hybrid 
response maps by augmenting the B-side data with the scaled model in the low signal-to-noise 
regime. (We plan to update the A-side model and produce an independent B-side model in 
the near future.) These hybrid beam response maps are available on LAMBDA as part of 
the WMAP three-year data release. The radial beam profiles obtained from these maps have 
been fit to a model consisting of a sum of Hermite polynomials that accurately characterize 
the main Gaussian lobe and small deviations from it. The model profil es are then Legendr e 
transformed to obtain the beam transfer functions h\ for each DA i (IJarosik et al.ll2006l ). 
These transfer functions are also provided with the three-year data release. 

The constraints imposed by the physical optics model have allowed us to extend the 
beam analysis out to much larger radii than was possible with the first-year analysis. As a 
result, we have determined that the beam solid angles were being systematically underes- 
timated by ~1%. Since we normalize the transfer functions to 1 at I — 1 this meant that 
the normalized functions were being systematically over-estimated by a comparable amount 
for Z>50. Consequently, as we discuss in §7.61 the final three-year power spectrum is ~l-2% 
higher than the first-year spectrum for Z>50. We believe this new determination of the beam 
response is more accurate than that given in the first-year analysis, but until we can complete 
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the B -side model analys is, we have held the window function uncertainties at the first-year 



level ([Page et al.l l2003al ). The first-year uncertainties are large enough to encompass the 
differences between the previous and current estimates of 6/. The propagation of random 
beam errors through to the power spectrum Fisher matrix is discussed in Appendix |A] 

An additional source of systemati c error in our treatment of the beam response arises 



from non-circularity of the main beam (jWu et al.ll200ll ; ISouradeep fc Ratrall2001 



Mitra et al. 



20041 ) . The effects of non-circularity are mitigated by WMAP's scan strategy, which causes 
most sky pixels to be observed over a wide range of azimuth angles. However, residual 
asymmetry in the effective beam response on the sky will, in general, bias estimates of the 
power spectrum at high I. In Appendix |B] we compute the bias induced in a cross-power 
spectrum due to residual beam asymmetry. In principle the formalism is exact, but in order 
to make the calculation tractable, we make two approximations: (1) that WMAP's scan 
coverage is independent of ecliptic longitude, and (2) that the azimuthal structure in the 
WMAP beams is adequately described by modes up to azimuthal quantum number m = 16. 
We define the bias as the ratio of the true measured power spectrum to that which would 
be inferred assuming a perfectly symmetric beam, 



Oil 



where Cf d is a fiducial power spectrum, h\ is the symmetrized beam transfer function noted 
above, and G\ v is the coupling matrix, defined in Appendix [B], that accounts for the full 
beam structure. 

Plots of ol\ for selected channel combinations are shown in Appendix [B] These results 
assume three-year sky coverage and beam multipole moments fit to the hybrid beam maps 
with Z max = 1500 and m max = 16. As shown in the Figure, the Q-band cross-power spectra 
(e.g., QlxQ2) require the greatest correction (6.3% at / = 600). This arises because the Q- 
band beams are farthest from the optic axis of the telescope and are thus the most elliptical 
of the three high frequency bands; for example, the major and minor axes of the Q2A beam 
are 0?60 and 0?42 (FWHM) respectively. The corrections to the V and W-band spectra 
are all <1% for / < 1000. Because of the relatively large bias in Q-band, and because of 
the potential to confuse it with the frequency- dependent point source correction, Q-band 
data was excluded from the final three- year power spectrum. The remaining V and W-band 
cross-power spectra were not corrected because the bias is already substantially smaller than 
the random instrument noise in the pertinent / range. 
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7.1.2. Instrument Noise Properties 

The noise bias term in equation fl39l) is the noise per ai m coefficient on the sky. If auto- 
power spectra are used in the final power spectrum estimate, the noise bias term must be 
known very accurately because it exponentially dominates the convolved power spectrum at 
high I. If only cross-power spectra are used, the noise bias is only required for estimating 
errors. Our final combined spectrum is based only on cross-power spectra for I > 30; while 
for I < 30, the noise bias is negligible compared to the sky signal, so that it need not be 
known to high precision for the maximum likelihood estimation (see §7.4p . 

In the limit that the time-domain instrument noise is white, the noise bias will be a 
constant, independent of I. If the noise has a 1/f component, the bias term will rise at low I. 
While the WMAP radiometer noise is nearly white by design, with 9 out of 10 differencing 



assemblies having 1/f knee frequencies of <10 mHz (I Jarosik et al.l 12003 bl ) . deviations of 



(Nf) from a cons tant must be accounted for, especially for the low-/ polarization analysis 



(IPage et al.ll2006l ). In the first-year analysis, we used end-to-end Monte Carlo simulations to 
generate realizations of the noise bias, then fit the mean of the realizations to a parameterized 
model. In the present analysis we use two sets of independent single-year sky maps from 
each DA to form the following estimate of the noise bias 

2 



w) = I E [ c i m - G i ii,ij,) ] ■ ( 41 ) 



In words, we take the difference between the auto-power spectra (years j = f) and the 
cross-power spectra (years j ^ j') for a given DA to estimate the noise bias. We then fit the 
estimate to a model of the form 

(Nf) =N + N 1 (^j , (42) 

where Nq and Ni are fit coefficients (we enforce Ni > 0), and the fit is performed over the 
range 33 < / < 1024. This model is used to compute the weights that enter into the final 
combined spectrum ( §7. 5[) and to estimate its noise properties. 



7.1.3. Systematic Errors 



Jarosik et al.l (120061 ) present limits on systematic errors in the three-year sky maps. 
They consider the effects of absolute and relative calibration errors, artifacts induced by 
environmental disturbances (thermal and electrical), errors from the map-making process, 
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pointing errors, and other miscellaneous effects. The combined uncertainty due to relative 
calibration errors, environmental effects, and map-making errors are limited to <29 /zK 2 
(2a) in the quadrupole moment C2 in any of the 8 high-frequency DAs. We estimate the 
absolute calibration uncertainty in the three-year WMAP data to be 0.5%. 

The noise in the WMAP sky maps is weakly correlated from pixel-to-pixel due to 1/f 
noise in the radiometers, and to the map-making process that infers pixel temperatures from 
filtered differential data. Neglecting these correlations is a form of systematic error that 
must be quantified. As part of the three-year analysis, we have developed code to compute 
the pixel-pixel inverse covariance matrix at pixel resolution r4. The resulting information 
is propagated through the com putation of the Fisher matrix to estimate power spectrum 



uncertainties (IPage et al.ll2006l ). Sample C\ uncertainty results for single- year cross-power 
spectra are shown in Figure [12j The red curves show the result when pixel-pixel noise 
correlations are ignored. The smooth rise at low I reflects our approximate representation 
of 1/f noise in the noise bias model, equation (T4"2]) . The black curves account for the full 
structure of the pixel-pixel inverse covariance matrix, including 1/f noise and map-making 
covariance. For the TT spectrum (left pair of curves in each panel), the noise is negligible 
compared to the signal (~1000 /iK 2 ), so this structure can be safely ignored. 

Random pointing errors are accounted for in the beam mapping procedure; the beam 



transfer functions presented in lJarosik et al.l (120061 ) incorporate random pointing errors auto 



matically. A systematic pointing error of ~1' at the spin period is suspected in the quaternion 
solution that defines the spacecraft pointing. This is much smaller than the smallest beam 
width (~12' at W band), and we estimate that it would produce <1% error i n the angular 



power spectrum at / = 1000; thus we do not attempt to correct for this effect. iBarnes et al. 



( 120031 ) place limits on spurious contributions due to stray light pickup through the far side- 
lobes of the instrument. As shown in Figure 4 of that paper, they place limits of <10 /iK 2 
on spurious contributions to 1(1 + l)C//27r, at Q through W band, due to far sidelobe pickup. 



Our template-based approach to Galactic foreground subtraction is detailed in §5.3j 
the effect of this frequency-by-frequency cleaning on the power spectra is shown in §7.21 
Diffuse foreground emission is a modest source of contamination at large angular scales ( > 2°). 
Systematic errors on these angular scales are negligible compared to the (modest) level of 
foreground emission. On intermediate angular scale (~2°), the 1-2% uncertainty in the beam 
transfer functions is the largest source of uncertainty, while for multipole moments greater 
than ~400, random white noise from the instrument is the largest source of uncertainty. 
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7.2. Galactic and Extragalactic Foregrounds 



In this subsection we determine the level of foreground contamination in the angular 
power spectrum. On large angular scales (Z ^ 50) the primary source of contamination is 
diffuse emission from the Milky Way, while on small scales (Z J; 400) the primary culprit is 
extragalactic point source emission. To simplify the discussion, we present foreground results 
using the pseudo-C/ estimate of the power spectra over the entire /-range under study, despite 
the fact that we use a maximum likelihood estimate for the final low-/ power spectrum. Thus, 
while the low-/ estimates presented in the section may not be optimal, they are consistently 
applied as a function of frequency, so that the relative foreground contamination is reliably 
determined. 

Figure [13] shows the cross-power spectra obtained from the high-frequency band maps 
(Q-W) prior to any foreground subtraction. All DA/yr combinations within a band pair 
were averaged to form these spectra, and the results are color coded by effective frequency, 
\j v i v % \ where v % is the frequency of differencing assembly i. As indicated in the plot legend, 
the low frequency data are shown in red, and so forth. To clarify the result, we plot the ratio 
of each frequency band spectrum to the final combined spectrum (see §7.51) . The red trace 
in the top panel shows very clearly that the low frequency data are contaminated by diffuse 
Galactic emission at low / and by point sources at higher /. The contamination appears to 
be least at the highest frequencies, as expected in this frequency range. This is consistent 
with the dominance of synchrotron and free-free emission over thermal dust emission. 

Not surprisingly, the most contaminated multipoles are / = 2 Sz 4, which most closely 
trace the Galactic plane morphology. Specifically, the total quadrupolar emission at Q-band 
is nearly 10 times brighter (in power units) than the CMB signal, while at W-band it is 
nearly 5 times brighter. The foreground emission reaches a minimum near V-band where 
the total / = 2 emission is less than a factor of 2 brighter than the CMB, and even less for 
/ > 2. Note also the modest foreground features in the range / ~ 10 — 30. 

The contribution from extragalactic radio sources can be seen in Figure d3] as excess 
emission in Q-band at high /. As discussed in §6.11 we performed a direct search for sources 
in the WMAP data and found 323. Based on this search, we augmented the mask we use for 
CMB analysis; so the contribution shown in the Figure is mostly due to sources just below 
our detection threshold. In the limit that these sources are not clustered, their contribution 
to the cross-power spectra has the form 



where A is an overall amplitude, measured antenna temperature, the factors gt convert the 




(43) 
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result to thermodynamic temperature, vq = 40.7 GHz (note that this differs from our first- 
year definition of 45 GHz), and we assume a frequency spectrum (3 = —2.0. The amplitude 
A is determined by fitting to the Q -, V-, and W-b a nd cro ss-power spectra. Following the 
procedure outlined in Appendix B of iHinshaw et al.l (l2003bl ). we find A = 0.014 ±0.003 /iK 2 
sr, in antenna temperature. (As discussed below, this value has been revised down from 
0.017 ± 0.002 /iK 2 sr since the original version of this paper appeared.) In the first-year 
analysis we found A = 0.022 /iK 2 sr (0.015 /^K 2 sr, referred to 45 GHz). The three-year 
source amplitude is expected to be lower than the first-year value due to the enlargement of 
the three-year sour ce mask and the cons equent lowering of the effective flux cut applied to the 
source population. ISpergel et al.l (120061 ) evaluate the bispectrum of the WMAP data and are 
able to fit a non-Gaussian source component to a particular configuration of the bispectrum 
data. They find the source model, equation (T4"3"j) . fits that data as well. We thus adopt the 
model given above with A = 0.014 /iK 2 sr and use the procedure given in Appendix IA1 to 
marginalize the over the uncertainty in A. At Q band, the correction to 1(1 + 1)Ci/2tt is 608 
and 2431 /iK 2 at Z=500 and 1000, respectively. At W band, the correction is only 32 and 
128 /iK 2 at the same I values. For comparison, the CMB power in this I range is ~ 2000 
/iK 2 . 



Since the release of the three-year data, iHuffenberger et al.l (120061 ) have reanalyzed the 



multi-frequency spectra for residual sources using a similar methodology. They find a some- 
what smaller amplitude of 0.011 ±0.001 /iK 2 sr. Differences in the inferred source amplitude 
at this level affect the final power spectrum by about 1%, which is also about the level at 
which second-order effects due to beam asymmetry can alter the spectrum. We expect our 
understanding of these effects to improve significantly with additional years of data. In the 
meantime, our revised estimate of A = 0.014 ± 0.003 uK 2 sr encompasses both our original 



Spergel et all (I2006h . 



Huffenbereer et al. 


( 


2006 


) estimate. r . 


ers is discussed in 


Huffenbereer et al. 



The bottom panel of Figure [13] shows the band-averaged cross-power spectra after sub- 
traction of the template-based Galactic foreground model and the above source model. The 
Q-band spectra exhibit clear deviations of order 10% from the V and W-band spectra at 
low I, while the higher frequency combinations all agree with each other to better than 5%. 
Also, while subtracting the source model in equation (T43|) brings the high-Z Q-band spec- 
trum into good agreement with the higher frequency results up to / ~ 400, beyond this point 
the 0?48 angular resolution of Q-band limits the sensitivity of this band. This is also the / 
range where the effects of the Q-band beam asymmetry start to bias our pseudo-C^-based 
power spectrum estimates. Since the V and W-band data alone provide a cosmic variance 
limited measurement of the power spectrum up to I = 400 (see §7.5[) we have decided to 
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omit Q-band data entirely from our final power spectrum estimate. But note that Q band 
serves a valuable role in fixing the amplitude of the residual point source contribution, and 
in helping us to assess the quality of the Galactic foreground subtraction. 



7.3. Low-/ Moments From The ILC Map (/ = 1, 2, 3) 

Based on our analysis of the ILC method presented in §5.21 we conclude that the newly 
debiased three-year ILC map is suitable for analysis over the full sky up to I < 10, though we 
have not performed a full battery of n on- G aussian tests on this map, so we still advise users 
to exercise caution. iTegmark et al.l (120031 ) arrived at a similar conclusion based on their 
foreground analysis of the first-year data. In this section we use the ILC map to evaluate 
the low-/ ai m coefficients by direct full-sky integration of equation (1321) . In §7.41 we estimate 
the low-/ power spectrum using a maximum likelihood estimate based on the ILC map using 
only data outside the Kp2 sky cut. We compare the two results, as a cross-check, in the 
latter section. 



Sky maps of the modes from / = 2 — 8, derived from the ILC map, are shown in Figured^ 
There has been considerable comment on the non-random appearance of these modes. We 
discuss this topic in more detail in §|E1 but note here that some non-random appearances 
may be deceiving and that high-fidelity simulations and a critical assessment of posterior 
bias are very important in assessing significance. 



7.3.1. Dipole (1=1) 



Owing to its large amplitude and its role as a WMAP calibration source, the / = 1 dipole 
signal requires specia l handling in t h e WM AP dat a processing. The c alibration process is 
described in detail in iHinshaw et al.l (l2003bl ) and in iJarosik et al. r|2006h . After establishing; 
the calibration, but prior to map-making, we subtract a dipole term from the time-ordered 
data to minimize signal aliasing that would arise from binning a large differential signal into 
two finite size pixels. In the first-year analysis we removed the COB E-determined dipole from 
the data, then fit the final maps for a residual dipole. From this, a best-fit WMAP dipole 
was reported, which was limited by the 0.5% calibration uncertainty. For the three-year data 
analysis we subtract the WMAP first-year dipole from the time-ordered data and fit for a 
residual dipole in the final ILC map. The results are reported in Table [6j The errors reported 
for the ai m coefficients in the body of the Table are obtained from Monte Carlo simulations 
of the ILC procedure that attempt to include the effect of Galactic removal uncertainty. In 
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the Table end notes we report the dipole direction and magnitude. The direction uncertainty 
is dominated by Galactic removal errors, while the magnitude uncertainty is dominated by 
the 0.5% calibration uncertainty. 

7.3.2. Quadrupole (1=2) 

The quadrupole moment computed from the full-sky decomposition of the ILC map is 
given in Table [6j (The released ILC map, and the quadrupole moment reported here have 
not been corrected for the 1.2 /iK "kinematic" quadrupole, the second-order Doppler term.) 
As with the dipole components, the errors reported for the ai m coefficients are obtained from 
Monte Carlo simulations that attempt to include the effect of Galactic removal uncertainty. 
The magnitude of this quadrupole moment, computed using equation (1341) . is AT| = 249 
/iK 2 [AT 2 = 1(1 + 1)Ci/2tt]. Here, we do not attempt to correct for the bias associated with 
this estimate, and we postpone a more complete error analysis of AT| to §7.41 However, we 
note that this estimate is quite consistent with the maximum likelihood estimate presented 
in §E3J 

A corresponding analysis of the first-year ILC map gives AT| = 196 /iK 2 . The increase 
in the three-year amplitude relative to this value comes mostly from the ILC bias correction 
and, to a lesser degree, from the improved three-year gain model. When our new ILC 
algorithm is applied to the first-year sky maps, we get AT| = 237 /iK 2 . While this 41 /iK 2 
bias correction is small in absolute terms, it produces a relatively large fractional correction 
to the small quadrupole. Further discussion of the low quadrupole amplitude is deferred to 
§731 

7.3.3. Octupole (1=3) 

The octupole moment computed from the full-sky decomposition of the ILC map is given 
in Table [6j The reported errors attempt to include the uncertainty due to Galactic foreground 
removal errors. The magnitude of this octupole moment, computed using equation (1541) . is 
AT 3 2 = 1051 ^K 2 . As above, we do not attempt to correct for the bias associated with this 
estimate, but again we note that it is quite consistent with the maximum likelihood estimate 
presented in §7.41 



It has been noted b y several authors that the orient ation of the quadrupole and octopole 
are closely aligned, e.g., Ide Oliveira-Costa et al.l (j2004al ). and that the distribution of power 
amongst the ai m coefficients is possibly non-random. We discuss these results in more detail 
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in §H1 but we note here that the basic structure of the low I modes is largely unchanged from 
the first-year data. Thus we expect that most, if not all, of the "odd" features claimed to 
exist in the first-year maps will survive. 



7.4. Low-/ Spectrum From Maximum Likelihood Estimate (1 = 2 — 30) 



If the temperature fluctuations are Gaussian, random phase, and the a priori probability 
of a given set of cosmological parameters is uniform, the power spectrum may be optimally 
estimated by maximizing the multi-variate Gaussian likelihood function, equation ( 1371) . This 
approach to spectrum estimation gives the optimal (minimum variance) estimate for a given 
data set, and a means for assessing confidence intervals in a rigorous way. Unfortunately, 
the approach is computationally expensive for large data sets, and for />30, the results do 



not significantly improve o n the quadratic pseudo-C; based estimate used in §7.51 (jEfstathiou 
2004bl ; lEriksen et al.ll2006l ). In this section we present the results of a pixel-space approach 
to spectrum estimation for the multipole range / = 2 — 30 and compare the results with the 
pseudo-C; based estimate. In Appendix Owe discuss several aspects of maximum likelihood 
estimation in the limit that instrument noise may be ignored (a very good approximation 
for the low-/ WMAP data). In this case, several simple results can be derived analytically 
and used SIS db guide to a more complete analysis. 

We begin with equat ion (l3Tj) and work in pixel space, where the data, d = t is a sky map. 
As in lSlosar et al.l (120041 ). we account for residual Galactic uncertainty by marginalizing the 
likelihood over one or more Galactic template fit coefficient (s). That is, we compute the 
joint likelihood 

nr Itl exp[-l(t- a . tf fC-'(t- a -t,)] 

HQ, a 1 1) oc -f= , (44) 

where tf is one or more foreground emission templates, and a is the corresponding set of fit 
coefficients, and we marginalize over a 



L(d\t) 



da L(Ci, a | t) 



(45) 



This can be evaluated analytically, as in Appendix B of iHinshaw et al.l (J2003bJ). 



In practice, we start with the r9 ILC map, we further smooth it with a Gaussian kernel of 
width 9?183 (FWHM); degrade the map to pixel resolution r4; mask it with a degraded Kp2 
mask (accepting all r4 pixels that have more than 50% of its r9 pixels outside the original 
Kp2 cut); then add 1 fiK rms of white noise t o each pixel to aid n umerical regularization of 
the likelihood evaluation. (See Appendix A of lEriksen et al.l (120061 ) for a complete discussion 
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of the numerical aspects of likelihood evaluation.) We evaluate the marginalized likelihood 
function with the following expression for the covariance matrix 

C(m, n,) = 2^ C i w i P ^ ■ n j) + N - ( 46 ) 

I 

Here wf is the effective window function of the smoothed map evaluated at low pixel res- 
olution, Pi(cos6) is the Legendre polynomial of order /, and N is a diagonal matrix that 
describes the 1 fiK rms white noise that was added to the map. This expression is evaluated 
to the Nyquist sampling limit of an r4 map, namely Z max = 47. The foreground template 
we marginalize over is the map T V i — T IL c, where T V1 is the VI sky map and T ILC is the 
ILC map. As discussed below, we have tried usin g numerous other com binations of data 



and sky cut and obtain consistent estimates of Ci. lEriksen et al.l (120061 ) have also studied 
this question and obtain similar results using a variety of data combinations and likelihood 
estimation methodologies. 

The full likelihood curves, L(C{) for / = 2 — 10, are shown in Figure [15] and listed in 
Table H (Note that the likelihood code delivered with the three-year data release employs 
the abov e method up to I = 30, then uses the quadratic form employed in the first-year 



analysis (jVerde et al.l 120031 ) for I > 30.) The predicted C; values from the best-fit ACDM 
model (fit to WMAP data only) are shown as vertical red lines, while the values derived 
from the pseudo-C; estimate ( §7.5[) are shown as vertical blue lines. Maximum likelihood 
estimates from the ILC map with and without a sky cut are also shown, as indicated in the 
Figure. There are several comments to be made in connection with these results: 

(1) All of the maximum likelihood estimates (the black curve and the vertical dot-dashed 
lines) are all consistent with each other, in the sense that they all cluster in a range that 
is much smaller than the overall 68% confidence interval, which is largely set by cosmic 
variance. We conclude from this that foreground removal errors and the effects of masking 
are not limiting factors in cosmological model inference from the low-Z power spectrum. 
However, they may still play an important role in determining the significance of low-/ 
features beyond the power spectrum. 

(2) The Ci values based on the pseudo-C/ estimates (shown in blue) are generally con- 
sistent with the maximum likelihood estimates. However, the results for I = 2, 3, 7 are all 
nearly a factor of 2 lower than the ML values, and lie where the likelihood function is roughly 
half of its peak value. This discrepancy may be related to the different assumptions the two 
methods make regarding the distribution of power in the cii m when faced with cut-sky data. 
Both methods have been demonstrated to be unbiased, as long as the noise properties of the 
data are correctly specified, but the maximum likelihood estimate has a smaller variance. 
Additionally, we note that the two methods are identical in the limit of uniformly-weighted, 
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full-sky data; and the maximum likelihood estimate based on cut-sky data is consistent with 
the full-sky estimate. In light of this, we adopt the maximum likelihood estimate of C/ for 
I = 2 - 30. 

(3) The best-fit ACDM model lies well within the 95% confidence interval of L(C{) for 
all I < 10, including the quadrupole. Indeed, while the observed quadrupole amplitude 
is still low compared to the best-fit ACDM prediction (236 /iK 2 vs. 1252 /iK 2 , for AT 2 ) 
the probability that the ensemble-average value of AT| is as large or larger than than the 
model value is 16%. This relatively high probability is due to the long tail in the posterior 
distribution, which reflects the xl distribution (i> — 21 + 1) that Gaussian models predict 
the C\ will follow. Clearly the residual uncertainty associated with the exact location of the 
maximum likelihood peak in Figure [J5] will not fundamentally change this situation. 

There have been several other focused studies of the low-/ power spectrum in the first- 



nature as the largest observable structure in 


the universe 


(Efstathiou 


2003; 


Gaztanaea et al. 


2003 


> 


Efstathiou 


2003: 


Bielewicz et al. 


2004 




Efstathiou 


2004al: 


Slosar et al. 


2004 


). These 



authors agree that the quadrupole amplitude is indeed low, but not low enough to rule out 
ACDM. 



7.5. High-/ Spectrum From Combined Pseudo-Q Estimate (/ = 30 — 1000) 
The final spectrum for / > 30 is obtained by forming a weighted average of the individua l 



cross-power spectra, equation (I3"5j) . computed as per Appendix A in lHinshaw et al.l (l2003bl ). 
For the three-year analysis, we evaluate the constituent pseudo-a^ m data using uniform pixel 
weights for / < 500 and Af b s weights for / > 500. For the latter we use the high-resolution rlO 
maps to reduce pixelization smearing. Given the pseudo-a^ data, we evaluate cross-power 
spectra for all 153 independent combinations of V- and W-band data ( §7.11) . 

The weights used to average the spectra are obtained as follows. First, the diagonal 
elements of a Fishe r matrix, F™, a r e comp uted for each DA i, following the method outlined 



in Appendix D of iHinshaw et al.l (l2003bl ). Since the sky coverage is so similar from year 
to year (see Figure [2]) we compute only one set of F™ for all three years using the year-1 
A^ bg data. To account for 1/f noise, the low-/ elements of each Fisher matrix are decreased 
according to the noise bias model described in §7.1.21 F™ gives the relative noise per DA 
(as measured by the noise bias) and the relative beam response at each /. We weight the 
cross-power spectra using the product F^ = (Fjf'Ff 1 ) 1 / 2 . 

The noise bias model, equation (|42|) . is propagated through the averaging to produce 
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an effective noise bias 



N, 



eff 



E«W)-W) 



(47) 



where N{ is the noise bias model for DA i, and the sum is over all DA pairs used in the final 
spectrum. In the end, the noise error in the combined spectrum is expressed as 



AG 



^eff 



21 + l / sk y ' 



(4S 



Here, f f ky is the effective sky fraction observed, which we calibrate using Monte Carlo simu- 
lations (IVerde et al.ll2003l ). The factor a is 1 for TE, TB, and EB spectra, and 2 for TT, EE, 
and BB spectra. Note that the final error estimates are independent of the original Fisher 
elements, F™, which are only used as relative weights. Note also that N° s is separately 
calibrated for both regimes of pixel weighting, uniform and N i, s . 

The full covariance matrix for the final spectrum includes contributions from cosmic 
variance, instrument noise, mode coupling, and beam and point source uncertainties. We 
have revised our handling of beam and source errors in the three-year analysis, as discussed 
in Appendix IA1 However, the remaining contribut ions are treated as we did in the first-year 
analysis (IHinshaw et al.ll2003bl ; I Verde et al.l 120031 ) . A good estimate of the uncertainty per 
C\ is given by 

' " (Cf d + iVf ff ), (19) 



AC; 



/sky 



21 + 1 



where Cf d is a fiducial model spectrum. This estimate includes the effects of cosmic variance, 
instrument noise, and mode coupling, but not beam and point source uncertainties. However, 
these effects are accounted for in the likelihood code delivered with the three-year release. 



7.6. The Full Power Spectrum 

To construct the final power spectrum we combine the maximum likelihood results from 
Table [7] for I < 30 with the pseudo-C based cross-power spectra, discussed above, for I > 30. 
The results are shown in Figure [16] where the WMAP data are shown in black with noise- 
only errors, the best-fit ACDM model, fit to the three-year data, is shown in red, and the 
la error band due to cosmic variance shown in lighter red. The data have been averaged 
in I bands of increasing width, and the cosmic variance band has been binned accordingly. 
To see the effect that binning has on the model prediction, we average the model curve in 
the same I bins as the data and show the results as dark red diamonds. For the most part, 
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the binned model is indistinguishable from the unbinned model except in the vicinity of the 
second acoustic peak and trough. 

Based on the noise estimates presented in §7.1.21 we determine that the three-year 
spectrum is cosmic variance limited to I = 400. The signal-to-noise ratio per /-mode exceeds 
unity up to I = 850, and for bins of width A/// = 3%, the signal-to-noise ratio exceeds unity 
up to I = 1000. In the noise-dominated, high-/ portion of the spectrum, the three-year data 
are more than three times quieter than the first-year data due to (1) the additional years 
of data, and (2) the use of finer pixels in the V and W band sky maps, which reduces pixel 
smearing at high /. The xt °f the full p ower spectrum relat ive to the best-fit ACDM model 



is 1.068 for 988 D.O.F. (13 < / < 1000) flSpergel et all 120061 ). The distribution of x vs. / is 



shown in Figure [TO and is discussed further below. 

The first two acoustic peaks are now measured with high precision in the three-year 
spectrum. The 2nd trough and the subsequent rise to a 3rd peak are also well established. 
To quantify these results, we rep eat the mo d el-inde pendent peak and trough fits that were 



applied to the first-year data by iPage et al.l (]2003bl ). The results of this analysis are listed 



in Table [SJ We note here that the first two acoustic peaks are seen at / = 220.8 ± 0.7 
and / = 530.9 ± 3.8, respectively, while in the first-year spectrum, they were located at 
/ = 220.1 ± 0.8 and / = 546 ± 10. Table [8] also shows that the second trough is now well 
measured and that the rise to the third peak is unambiguous, but the position and amplitude 
of the 3rd peak are not yet well constrained by WMAP data alone. 

Figure [18] shows the three-year WMAP spectrum compared to a set of recent balloon 
and ground-based measurements that were selected to most complement the WMAP data in 
terms of frequency coverage and / range. The non- WMAP data points are plotted with errors 
that include both measurement uncertainty and cosmic variance, while the WMAP data in 
this / range are largely noise dominated, so the effective error is comparable. When the 
WMAP data are combined with these higher resolution CMB measurements, the existence 
of a third acoustic peak is well established, as is the onset of Silk damping beyond the 3rd 
peak. 

The three- year spectrum is compared to the first-year spectrum in Figure [TJB We show 
the new spectrum in black and the old one in red. The best-fit ACDM model, fit to the three- 
year data, is shown in grey. In the top panel, the as-published first-year spectrum is shown. 
The most noticeable differences between the two spectra are, (1) the change at low-/ due to 
the adoption of the maximum likelihood estimate for / < 30, (2) the smaller uncertainties in 
the noise-dominated high-/ regime, discussed further below, and (3) a small but systematic 
difference in the mid-/ range due to improvements in our determination of the beam window 
functions ( §7.1.11) . The middle panel shows the ratio of the new spectrum to the old. For 
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comparison, the red curve shows the (inverse) ratio of the three-year and first-year window 
functions, which differ by up to 2%. The spectrum ratio tracks the window function ratio 
well up to / ~ 500 at which point the sensitivity of the first-year spectrum starts to diminish. 
For I < 30 in this panel, we have substituted the pseudo-C; based spectrum from the three- 
year data into this ratio to show the stability of the underlying sky map data. When based 
on the same estimation method, the two spectra agree to within ~2%, despite changes in the 
gain model and in the details of the Galactic foreground subtraction, both of which affect 
the low-/ data. (The gain model changes are important for the polarization data, while the 
temperature-based foreground model has no direct bearing on polarization.) The bottom 
panel shows the three-year and first-year spectra again, but this time the first-year data 
have been deconvolved with the three-year window functions and we have substituted the 
maximum likelihood estimate into the first-year spectrum. The agreement between the two 
is now excellent. 

As noted above, the power spectrum is now measured with cosmic variance limited 
sensitivity to I = 400, which happens to coincide with the measured position of the first 
trough in the acoustic spectrum (Tabled]). Figure [2171 shows the measurement of the first 
acoustic peak with no binning in /. The black trace shows the three-year measurement, while 
the grey trace shows the first-year result for comparison. The background error band gives 
the la uncertainty per I for the three-year data, including the effects of Galaxy masking and 
a minimal contribution from instrument noise at the high-/ limit. In the first-year power 
spectrum there were several localized features, which have become known, technically, as 
"glitches". The most visible was a "bite" near the top of the first acoustic peak, from 
/ = 205 to 210. Figure [20] shows that the feature still exists in the three-year data, but it is 
not nearly as prominent, and it disappears almost entirely from the binned spectrum. In this 
/ range, the pixel weights used to evalu ate the pseudo-C'/ data changed from "transitional" 
in the first-year analysis [Appendix A of iHinshaw et al.l (J2003bj)] to uniform in the three-year 
analysis. This reduces the effective level of cosmic variance in this / range, while increasing 
the small noise contribution somewhat. We conclude that the feature was likely a noise 
fluctuation superposed on a moderate signal fluctuation. 

Several other glitches remain in the low-/ power spectrum, perhaps includi ng the low 
quadrupole. Several authors have commented on the signi ficance of these features (jEfstathiou 
20031 ; lLewisll2004l ; iBielewicz et al.ll2004l ; ISlosar et al.ll2004l ) , and several more have used the Cy 
data t o search for features in the unde r lying primordial spectrum, P(k ) (IShafieloo fc Souradeep 



2004]: Martin & Ringevall I2004L 120051 : iTocchini-Valentini et al.ll2006T ). see also ISpergel et al. 



(120061 ) . Figure [T7| shows that the binned \ 2 P er /-band is slightly elevated at low /, but 
not to a compelling degree. In the absence of an established theoretical framework in which 
to interpret these glitches (beyond the Gaussian, random phase paradigm), they will likely 
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remain curiosities. 



8. BEYOND THE ANGULAR POWER SPECTRUM - LARGE-SCALE 

FEATURES 

The low / CMB modes trace the largest structure observable in the universe today. In 
an inflationary scenario, these modes were the first to leave the horizon and the most recent 
to re-enter. Since they are largely unaffected by sub-horizon scale acoustic evolution, they 
also present us with the cleanest remnant of primordial physics available today. Unusual 
behavior in these modes would be of potentially great and unique importance. 

Indeed there do appear to be some intriguing features in the data, but their significance is 
difficult to ascertain and has been a topic of much debate. In this section we briefly review the 
claims that have been made to date and comment on them in light of the three-year WMAP 
data. These features include: low power, especially in the quadrupole moment; alignment 
of modes, particularly along an "axis of evil"; unequal fluctuation power in the northern 
and southern sky; a surprisingly low three-point correlation function in the northern sky; an 
unusually deep/large cold spot in the southern sky; and various "ringing" features, "glitches" , 
and/or "bites" in the power spectrum. Of course one expects to encounter low probability 
features with any a posteriori data analysis, and there are no clear cut rules for assigning 
the degree of posterior bias to any given observation, so some amount of judgment is called 
for. Ultimately, the most scientifically compelling development would be the introduction of 
a new model that explains a number of currently disparate phenomena in cosmology (such 
as the behavior of the low / modes and the nature of the dark energy) while also making 
testable predictions of new phenomena. 



8.1. Summary of First-Year Results 



The 2-point correlation function computed from the first-year WMAP data showed a 
notable lack of signal on large ang ular scales (ppergel et al.l 120031 ) . Similar behavior was 
also seen in the COBE-DMR data (IHinshaw et al.lll996l ). making unidentified experimental 
systematic effects an unlikely cause. Because the 2-point function is the Legendre transform 
of the angular power spectrum, the large-scale behavior of C{6) is dominated by the lowest 
I modes, especially the quadrupole. So this signature is at least partially a reflection of the 
low quadrupole amplitude. 



To study this further, ISpergel et al 



(120031 ) characterized the feature by an integral, 



-50- 



S = j^ /2 C 2 (9)d cos8, which measures the power in C(9) for 9 > 60°. In a Monte Carlo 
simulation, only 0.15% of the realizations produced a value of S as small as did the WMAP 
data. It remains to assess the degree to which Gal actic errors and poster ior bias (the selection 
of the 60° integration limit) affect this result. iNiarchou et al.l (120041 ) concluded that the 
evidence for measurement error was at least as likely as the evidence for a cutoff in the 
primordial power spectrum, but that both in both cases the evidence was weak. 



Tegmark et al.l (120031 ) produced a high resolu tion full-sky CMB map by usi ng a variant 
of the foreground cleaning approach developed by iTegmark fc Efstathiou Jl996h . They cau- 
tiously quote a quadrupole amplitude of 202 /zK 2 , evaluated over the full-sky. They further 
remark that the quadrupole and octupole phases are notably aligned with each other and 
that the octupole is unusually "planar" with most of its power aligned approximately with 
the Galactic plane. Ide Oliveira-Costa et al.l (l2004bl ) estimate the a priori chance of observ- 
ing both the low quadrupole, the I — 2, 3 alignment, and the / = 3 planarity to be ~1 in 
24000. Note that this estimate does not formally attempt to account for G alactic modeling 



uncert ainty, which will tend to reduce the significance of the noted features. iBielewicz et al. 
(120041 ) confirmed the above result, and they further conclude that the I = 3 properties are 



stable with respect to the applied mask and the level of foreground correction, but that the 
quadrupole is much less so. 



Schwarz et al.l (120041 ) claimed that the quadrupole and octopole are even more correlated 
(99.97% CL), with the quadrupole plane and the three octopole planes "remarkably aligned." 
Further, they claimed that three of these planes are orthogonal to the ecliptic and the 
normals to these planes are aligned with the direction of the cosmological dipole and with 
the equinoxes. Th is had led to speculation that the low-Z signal is not cosmological in origin. 
Copi et al.l (120031 ) use "multipole vectors" to characterize the geometry of the I modes. They 
conclude that the "oriented area of planes defined by these vectors ... is incon sistent with the 
i sotrop ic Gaussian hypothesis at the 99.4% level for the ILC map." See also iKatz fc Weeks 
d2004h . 



Eriksen et al.l (120041 ). in their study of the ILC method, confirm that the quadrupole and 
octopole are strongly aligned. They also note that the I = 5 mode is "spherically symmetric" 
at ~ 3<j, and the / = 6 mode is planar at ~ 2cr confidence (see Figure [14]). But they add 
that the fi rst-year ILC map i s prob ably not clean enough to use for cosmological analyses of 
this type. lLand fc Magueijol (120051 ) point out that the / = 3 and 5 modes are aligned in both 
direction and azimuth, thereby "rejecting statistical isotropy with a probability in excess of 
99.9%." 



Hansen et al.l (j2004al ) fit cosmological parameters to the WMAP data separately in the 
northern and southern hemispheres in three coordinate systems. They conclude that, "it may 
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be necessary to question the assumption of cosmological isotropy." lEriksen et al.l (j2004al ). 
evaluate the ratio of low-/ power between two hemispheres and conclude that only 0.3% of 
simulated skies have as low a ratio as obser ved, even when allowing the (simulated) data 
to define the direction. lHansen et al.l (l2004bl ) reach a similar conclusion and note that it is 
"hard to explain in terms of residual foregrounds and known systematic effects." 

Faced with so many apparently unlikely features in the data, our foremost priority 
is to re-evaluate potential sources of systematic error in the maps. Indeed, a by-product 
of the exhaustive three-year polarization analysis is enhanced confidence in the accuracy 
of the temperature signal. Figure [3] shows that, even after all of the processing changes 
outlined in §j2J the three-year maps are consistent with the first-year maps up to a small 



quadru pole difference (Tabled]) that is well within the first-year error budget (IHinshaw et al. 
2003al ). Furthermore, the improvements applied to the ILC processing to reduce Galactic 



foreground residuals, did not visibly alter the low-/ phase structure in the three-year ILC 
map (Figure fl4"|) . nor the north-south asymmetry that is plainly visible in Figure 1211 We 
have made no attempt to evaluate the above-noted statistics using the three-year maps, so 
the degree to which they persist remains to be seen. 



9. SUMMARY AND CONCLUSIONS 



For the three-year WMAP data release, we have made improvements in nearly every 
aspect of the data processing pipeline, many of which were driven by the need to successfully 
characterize the instrument noise to the level required to measure the polarization signal, 
~0.1 /xK. The impr o vements are spelle d out in detail here and in the companion papers by 
Jarosik et al.l (120061 ). iPage et al.l (120061 ). and ISpergel et al.l (120061 ) . The key points follow. 



1. Improved models of the instrument gain and beam response are now accurate to better 
than 1% and do not presently limit the scientific conclusions that can be drawn from the 
data. Polarization-specific effects such as spurious pickup from radiometer bandpass 
mismatch have also been accounted for in the data processing, and do not limit the 
data. 



The map-making procedure has been overhauled to produce genuine maximum likeli- 
hood maps of the temperature and polarization signal. In concert with this, we have 
produced code to evaluate the corresponding pixel-pixel weight matrix (inverse covari- 
ance matrix) at pixel resolution r4 . This was required to adequately characterize the 
noise for the polarization analysis (I Jarosik et al.ll2006l ). 
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3. The three-year temperature maps are consistent with the first-year maps ($3J) and, to 
a good approximation, have 3 times lower variance. A by-product of the exhaustive 
polarization analysis is enhanced confidence in the accuracy of the temperature signal. 

4. We have updated the MEM, ILC, and template-based temperature foreground emission 
models and assessed their uncertainties. We have developed new models of the polar- 
ized foreground emission that allow us to extract the cosmolo gical reionization signal 
and that pave the way for future studies of CMB polarization (IPage et al.l 120061 ) . 



5. Our analysis of the temperature power spectrum is improved at low I by employing 
a maximum likelihood estimate. At high I, the spectrum is > 3 times more sensitive 
due to the additional years of data and the use of rlO sky maps to reduce the effects 
of pixel smoothing. We have also placed new limits on systematic effects in the power 
spectrum due to beam ellipticity. 

6. We have updated the WMAP point source catalog. The contribution of unresolved 
sources to the power spectrum has been updated and subtracted. 

7. We have developed new methods to evaluate polarization power spectra from sky map 
data. Our approach accounts for the following key issues: (1) sky cuts, (2) low signal 
amplitude, (3) correlated noise from 1/f and scan-related effects, (4) other polarization- 
specific effects, such as baseline sensitivity, and bandpass mismatch. Thi s effort was 



devel oped in conjunction with the processing of polarization sky maps (IPage et al. 



2006|). 



8. Figure 1221 which updates Figure 12 fro mlBennett et al.l (I2003bf) . shows the three-year 
TT and TE power spectra ( §7.61 and IPage et al.l (120061 )). The TT measurement is 
cosmic variance limited up to I — 400, and has a signal-to-noise ratio greater than 
one up to I = 1000, in / bands of width 3%. The high-/ TE signal is consistent with 
the first-year result, while for / < 6 the signal is reduced. A det ailed compa r ison o f 
the three-year TE signal to the first-year is shown in Figure 24 of IPage et al.l (120061 ). 
The implications for the inferred optical depth, r, are shown in their Figure 26, which 
includes a detailed comparison of the three-year results with the first-year. 

9. We have refined the evaluation of the likelihood function used to infer cosmological 
parameters. The new approach uses pixel-space inputs for the low-/ temperature and 
polarization data and c orrectly account s for the joint probability of the temperature 
and polarization signal (IPage et al.ll2006l ). We use this to infer cosmological parameters 
from the three-year data. Many parameter degeneracies that existed in the first-year 
results are greatly reduced reduced by having new EE data and more sensitive high-/ 
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TT data. Figure 1 of ISpergel et al.l (120061 ) illustrates this for the parameter pairs (n s ,r) 
and (n s ,flbh 2 ). A cosmological model with only 6 parameters still provides an excellent 
fit to WMAP and other cosmological data. 



10. With new constraints on r from the EE data, we get correspondingly tighter con- 
straints on the scalar spectral index, n s . Using WMAP data only, we now find n s = 
0.958 ± 0.016, while combining WMAP data with other data gives n s = 0.947 ± 0.015 
( ISpergel et al.ll2006l). Joint limit s on n s and the tensor to scalar ratio, r, are shown in 
Figure 14 of ISpergel etaD tood ). 



11. The WMAP observatory continues to operate flawlessly, and many results that are 
currently measured at 2-3 o confidence, such as the apparent deviation from scale 
invariance, hints of a running spectral index, and the details of the reionization history, 
will be significantly clarified by additional years of data. 



10. DATA PRODUCTS 

All of the three-year WMAP data products are being made available through the Legacy 
Archive for Microwave Background Data Analysis (LAMBDA), NASA's CMB Thematic 
Data Center. The low-level products include the time-ordered data (calibrated and raw), 
the functions used to filter the time-ordered data, and the beam response data. The pro- 
cessed temperature and polarization maps in both single- year and three-year forms, with and 
without foreground subtraction, are supplied at pixel resolutions r4, r9, and rlO. Full pixel- 
pixel inverse covariance matrices are supplied for analyzing the r4 polarization maps. The 
processed foreground products include the MEM component maps, the ILC map, and the 
temperature and polarization templates used for foreground subtraction. The analyzed CMB 
products include the angular power spectra, the likelihood used in the three-year parameter 
analysis, and a complete web-based set of parameter results for every combination of data set 
and cosmological model that was run for the three- year analysis. For each of the model/data 
combinations, we also supply the best-fit model spectra and the full Markov chains. Th e 



products are described in detail in the WMAP Explanatory Supplement (ILimon et al.ll2006l ). 
The LAMBDA URL is http://lambda.gsfc.nasa.gov/. 
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A. IMPLEMENTATION OF BEAM AND POINT SOURCE ERRORS 

In this Appendix we discuss our simplified handling of beam deconvolution uncertain- 
ties and point source subtraction errors. Since both of these errors primarily affect the 
high-/ spectrum, we adopt the approximation that the likelihood is Gaussian and develop 
methodology for computing the change in the likelihood induced by these errors. 



A.l. Likelihood Decomposition 

At all but the lowest multipoles, I, the likelihood of the data, d = Ci, given a theo- 
retical model, m = C/, is well approximated by treating the power spectrum as Guassian 
distributed, 

C = -21nL(d|m) = J^(Ci - Ci)T^}(d v - C v ) + lndetS (Al) 



where S is the covariance matrix of the data, £^ = (ACjACV). As discussed in lHinshaw et al 



( I2003bl ). the full covariance matrix may be decomposed into its constituent contributions 

S S C v ~\~ S no j se ~\~ S ma sk ~\~ ^beam ~\~ S src . v^^) 

Since we are only interested in the final two terms in this Appendix, we introduce a simplified 
notation whereby 

S = S + S 1; (A3) 

where So consists of the first three terms: cosmic variance, instrument noise and mode cou- 
pling due to the sky mask; while Si consists of the last two terms, the beam decovolution and 
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point source subtraction uncertainties. With this decomposition of the covariance matrix, 
the likelihood becomes 

C = - C ^ + ^)u>(Cv ~ C v ) + lndet(S + S x ). (A4) 

IV 

The key point in the new treatment is that the beam and point source uncertainties have 
a very restricted form in I space. The source model has the form given in equation ( )43l) . in 
which we treat A as uncertain, while the beam deconvolution uncertainty may be described 
by only a handful of modes. As discussed in detail below, this allows us to decompose Si in 
the form 

El « UU T , (A5) 

where U is an Ni x M matrix with M <C A; (typically M ~ 10). This approximation allows 
us to efficiently compute S -1 using the Sherman- Morrison- Woodbury formula, 

(So + Si)- 1 « (Eo + UUT 1 (A6) 
= So 1 -S 1 U(l + U T S 1 U) _1 U T S 1 . (A7) 

With this result the likelihood decomposes into £ = £ + £\, where 

£o = £(Ci - C0(S o )^(4 - C v ) + lndet S , (A8) 

XV 

and 

£i = - ^{Ci - Ci) [s l U (I + U T S 1 U) _1 U T S 1 ~\ f {C v - C v ) + lndet (I + U T S l \J). 

(A9) 

Note that we have rewritten the In det term, 

lndet (I + Eo 1 !!!^) = In det (I + l^E^U), (A10) 

to reduce the dimensionality of the matrix in the determinant from x JVj to M x M. As 
a further simplification, we take S to be diagonal when evaluating C\. We have checked 
that this has a negligible effect on the accuracy of the expression. 



A. 2. Beam Deconvolution Uncertainty 



The detailed form of the covariance matrix due to beam uncertainties was discussed in 



Hinshaw et al.l (l2003bl ). but we review the salient features here. For a pair of cross-power 
spectra, C\ and Cj, the beam covariance has the form 



Cj(Sb)«/CV, 



(All) 
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where 

(Sb)J = + V) + flJCfy + (A12) 

Here z, i' denote the pair of DA's in the spectrum C\, and similarly for j,j'. B\ v is the 
fractional covariance matrix, B\ v = (u\u\,), where u\ = A&j/fe] is the fractional error in the 
beam transfer function b\. The complete beam covariance matrix for the final combined 
spectrum is 

(£ b V = Q(S b VCV = ^a^(Sb)>iC„ (A13) 

ij 

where the w] are the weights used to form the combined spectrum, Cj = J2i w i^i ( see §7.5p . 

The beam covariance matrix is dominated by a small number of modes. We take its 
singular value decomposition, S5 = PwQ T , where P and Q are orthogonal matrices (P = Q, 
since S& is symmetric), and w = diag(tOi, . . . , iujvJ is diagonal with singular values W\ > . . . > 
wn v Only 10 modes have Wi/w\ > 10 -3 , and only 18 have Wi/w\ > 10~ 6 . Thus we can 
approximate E& as 

M 

(? b ) u , Rt^CtPuWiPfi'Cv, (A14) 
where M ^ N\ is the number of modes used in the approximation. Then, defining an iVj x M 

1 /2 7" 1 

matrix t/j fc = CiPuw/ , we have S fe w UU . For the three- year analysis we use M = 9; 
tests with additional modes changed the TT likelihood by A% 2 < 0.01. 



A. 3. Point Source Subtraction Uncertainty 

The amplitude of the point source correction is uncertain at the 20% level. To include 
this uncertainty in the likelihood evaluation we simply add another mode to the matrix U 
defined above. Specifically, U^m+i) = oaC/j 10 , where Cf IC = Y^i wl fii' STC ^ s the point source 
model defined in equation fH3l) . and a a = 0.21 is the fractional uncertainty of the point 
source amplitude A. 



B. BEAM ASYMMETRIES 



In iJarosik et al.l (120061 ) we outline our approach to determining the beam transfer func- 
tions from flight observations of Jupiter. This method implicitly assumes that the effective 
beam response is equivalent to the azimuthal average of the response, which is appropriate 
in the limit that each sky map pixel is observed with equal weight over a full range of az- 
imuth. In practice this is a good approximation near the ecliptic poles, but is less so near the 
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ecliptic plane. In this Appendix we quantify how WMAP's beam asymmetry distorts our 
estimate of the angular power spectrum by taking into account: (1) the intrinsic shape of the 
beam response in spacecraft-fixed coordinates, (2) the map-making algorithm, including the 
spacecraft scan strategy, and (3) the form of the estimator used to measure the cross-power 
spectra. 



Wu et al.l (120011) tre ated the MAXIMA beams using an effective symmetric beam. 



Wandelt fc Gorskil (120011 ) developed a fast method of convolving two arbitrary band-limited 
functions on the sphere by transforming it into a 3-dimensional Fourier transform. 



B.l. Formalism 

The calibrated, differential time-ordered data, d, is related to the sky map, t, via 

d = Mt + n, (Bl) 



where M is the mapping funct ion and n is the noise which we assume to have zero mean 



(n) = (IHinshaw et al.ll2003al ). The maximum likelihood estimate for t is t = Wd, where 
W = (M T N~ 1 M)~ 1 M T N~ 1 , and N = (nn T ) is the noise covariance matrix of the time- 
ordered data. Using the fact that WM = I, where I is the identity matrix, it follows that 
t = Wd = t + Wn. The maximum likelihood map is unbiased because the sky map noise, 
Wn, has zero mean. 

For a differential instrument like WMAP that observes temperature differences between 
two sides A and B, the mapping function M can be decomposed into two pieces, — M B , 
which denote the A- and B-side contributions, respectively. For a given side, the matrix 
elements Mj P have the form gi(p)Q p , where gi(p) is the beam response in pixel p for the ith 
observation, and Q p is the pixel solid angle. The normalization is such that J dQgi(p) — > 

In practice, the solution t = Wd is too costly to evaluate exactly given the properties 
of M and N _1 . Rather, we adopt approximate forms for these quantities, M' and N' _1 , and 
estimate the sky map using t' = Wd, where W = (M /T N' _1 M / ) _1 M' T N /_1 . The resulting 
sky map is biased relative to the true sky according to 

t' = WMt = Jt, (B2) 

where J = WM. (The distorted sky map noise, Wn, still has zero mean.) The map-making 
algorithm used by WMAP assumes that the beams are circularly symmetric with infinite 
resolution. That is, each row (observation) of M' contains a +1 in the column (pixel) seen 
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by the A-side beam and a —1 in the column (pixel) seen by the B-side beam. In addition, 
since the WMAP instrument noise is nearly white (N w (TqI) the algorithm solves for t' using 
W = (M /T M') _1 M' T . The matrix M' T M' is diagonally dominant with diagonal elements 
N p = N A + Np equal to the number of times pixel p has been observed by either the A- or 
B-side beam. The action of W on d is thus approximately 

W;^^(<$-<©, (B3) 

where 5^ is one if the A-side beam is in pixel p and zero otherwise, and likewise for 5^. 
Combining equation ( 1B8I) and equation ( 1B3I) . we obtain the following expression for the 
distortion matrix J 

£ 9f( P ') + £ gfip') J = j£, + J£, (B4) 

i\Aep i\Bep J 

Here "i\A G p" denotes summation over all observations i where the A-side beam observes 
pixel p. This algorithm produces a sky map in which the effective beam response is an 
average of the A- and B-side responses convolved with the scan pattern at each pixel. 




At small angular scales where beam asymmetries may become important, a pseudo-C, 



method is used to estimate the WMAP angular power spectrum from these maps (IHauser fc Peebles 



19731 : IWandelt et al.ll200ll ; iHivon et al J 12002 : Iffinshaw et alJl2003bh . This method estimates 



the power spectrum C\ by first transforming a weighted map Qt', where Q is a weight matrix 
used to mask regions of strong foreground contamination and/or high instrument noise. This 
produces the pseudo-C/ spectrum, C\, which is then corrected for the statistical bias induced 
by Q under the assumption that the underlying signal is statistically isotropic. 

In the present context we wish to correct the power spectrum for the biases induced 
by both Q and J. We transform the distorted map into the spherical harmonic basis, a' = 
fipY^Qt' = fipY^QJt = fipY^QJYa, where a is the (full-sky) spherical harmonic transform 
of t@ The covariance matrix of a' is then related to the underlying signal covariance by 
(a'a't) = JSJ\ where J = fi p YtQJY and S is the undistorted CMB signal matrix S = 
(aat) = Ci5(i m )(ii m i). The distorted pseudo-C; spectrum is thus related to the true power 
spectrum, in ensemble- average, by 

W> = 2/TT £<KJ 2 > = £ G "'^'>' 



^^Here Y is the pixelized spherical harmonic transform matrix with elements ~Y p n m ) = Yim(p)> where p is 
the pixel index. It obeys YY^ = 5 pp >/Q p and Y^Y = Sn m \ni m '\/Qp, where fi p is the pixel solid angle. 
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where 



G IV = 21+1 \^(lrn)(l'm')\ 2 (B6) 



is the mode coupling matrix. For cross-power spectra (aiaJ,), |J7"| 2 is replaced by J\J£ where 
J\ is evaluated with the weights and beam properties appropriate to channel 1, and so on. 
The unbiased power spectrum estimate is then G~[^C' V . 

We now turn to the explicit evaluation of the distortion matrix J . First, we assume 
that the time dependence of the beam response on the sky is due only to changes in the 
spacecraft orientation. That is, gi(p') = g{Pi,&i',p f ) where Pi is the pixel observed by the 
beam centroid during the ith observation, and a.i is the azimuth angle of the observation at 
that time. If we expand the spacecraft-fixed beam response in spherical harmonics, 

g(z,0;p) = ^g lm Y lm (p), (B7) 

Im 

then the response on the sky may be expressed as 

9i (j/) = D Lm'(Pi, Oi)9l m 'Y lm (p'), (B8) 

Imm' 

where the D l mm , are the Wigner D-functions. 

Since J" is a linear function of J, it follows that J(3 A + J B ) = J A + J B . Furthermore, 
if Q is diagonal, we can write 3 s (with S=A,B) as 



s _ Q p 
V _ N p 

P i\S£p 



£ 9fip') = ^vj ^w s (p,a)g s (p,a;p>), (B9) 



where 

w 



p i\sep 



The Wigner D-functions form an orthogonal basis for functions defined over the the 3- 
dimensional rotation group, SO(3), thus we can expand w s (p, a) as 

w s (p,a) = £ ^mm'D l mmf (p,a), (Bll) 

Imm' 

where the sum is over all integral values of l,m,m f . Then, using the product rule 

«)^(P. «) = E C\7m^mP^m'PmrAV, «), (B12) 
Imm' 
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where the C 1 ^^ are the Clebsch- Gordon coefficients, along with the fact that 



/ ^ D l mm ,(p, a) = \f^yL{p) Sm'o, (B13) 



we find 



QppJppi — ^ ^m; W m^^/i^ii2m2^m' 1 «2m2y 21+1 ^™^^ imi ^ ' (^^) 

lm,!imim'j ,l 2 m 2 m' 2 

Here we have defined bi m = a/47t/ (21 + 1) gi m , and we note that when the be a ms are circular 
bim = bi6 m0 , where bi is the beam transfer function defined in iPage et al.l (j2003al ). Now, 
recalling that J = Q P Y^Q3Y, we have 



7 s - h s * v, s ' h * r lm r 10 j zh—- mi^ 

,l 2 m 2 m' 2 

Note that the sums over m,2 and m' 2 effectively drop out of equation (1B15I) . since the Clebsch- 
Gordon coefficients are zero unless mi + 1712 = m and m! x + m' 2 = 0. 

Calculating Gw is now straightforward; for cross-power spectra, (aia^), we can expand 

it as 

G lv = + Gf^ + Gf,} M + G^ B \ (B16) 

where A\ refers to side A of channel 1, and so forth. Then, using 

^ Cv^xhmflhmil'zmiz = ^W^rmm'^ + 1)/(2Z 2 + 1), (B17) 



rami 



we can write 



where 



(tSi\W nS2*\U' 

nSiS 2 \ y \ JhrniV 1 - Jl 1 m 1 /ni n\ 

Gw =p ^1 > ( B18 ) 



ffihmi = Yl b l'rn' W m[m[ C lLn im []j^r^Y- ( B19 ) 

In the case of uniform sky coverage, equation (IB18I) reduces to a familiar limit. Specif- 
ically, if all pixels are observed for the same length of time with a uniform distribution of 
azimuth angles, then w(p, a) = 1, which gives w l mm , = This result, along with 

the identity C^ o0 = 8u>5 mm >, gives Gw = bf Q 5w- Thus, even if the beams are not circularly 
symmetric, the full azimuthal coverage effectively symmetrizes the response. 
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If we retain full azimuthal coverage but introduce a mask, w(p, a) = q(p), we have 

21 + 1 



w. 



in" 



21 + 1 
Air 



dtt daq(p) D^ m ,(p,a 



Qlm &m'0i 



(B20) 

(B21) 



where the q\ m are the spherical harmonic coefficients of q(p). Then, expressing the Clebsch- 
Gordon coefficients in terms of the Wigner-3j symbols, 



Vm'V'm" 



■l) m+l ' +l "^2T+l 



V V' 



m m 



I 



-m 



(B22) 



we can write 



I l V'm" ~ bl'O qv'm" ( — 1) 



V+V' 



V I" I 




(2i> + i)(2i" + r 

47T 



(B23) 



which gives the standard result for the coupling matrix fjHivon et al.ll2002l ) 

2 



21' + 1 2 

— 2^ — 7Z — ^'0 \Ql"m"\ 



l"m" 



An 



V I" I 




(B24) 



Note that the beam window function is naturally incorporated in this expression, because 
of the way the mapping matrix was defined. This result also generalizes to the case of 
non-uniform noise per pixel by allowing the mask function, q(p), to become a weight map. 
Finally, if the beams are circularly symmetric, we can drop the requirement that pixels be 
observed with uniform azimuthal weight and derive the same result. 



B.2. Results 

Evaluating equation (IB18I) requires 0(1^^) operations and is thus impractical when 
I m£H ~ 10 3 . However, if we ignore sky cuts and apply uniform sky weight, Q = I, we 
can exploit an approximate symmetry of the WMAP scan pattern to significantly speed 
up the calculation. Specifically, the WMAP scan strategy is approximately independent of 
ecliptic longitude, implying w s (9, 0, a) ~ w s (9, 0, a) where 9 and are ecliptic latitude and 
longitude, respectively. This eliminates the sum over mi in equation (1B18I) . 

The weight coefficients are 

= / dQdaw s (9, ( j>,a)D l * ml (9, ( j ) ,a), (B25) 
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and since 



4-7T 

D l * ml (9,<P,a) = \j ^_ m ,(0, a ), (B26) 



the weight coefficients simplify to 



^ = (B27) 



where 



= J d(cos6)daw s (6,a)Y l * m (6,a) (B28) 
is a standard spherical harmonic transform. Then 

(jS\W _ r / (2/' + l)(2l" +1) >^ g 5 ( V I" l\ v+v , 

I 1 )l"m" ~ °m"0\l ^ 2^ b l'm' W l"m' ^ m / _ m , q J l" 1 ) > ( B29 ) 



and 



^ 2 = ^ E fe^fc^f L l 'i l)( L V L \ )• (B30) 



47r i'm' l"m' I'm" w /"m" 1 / , Q I \ n n g 

«"m'm" V 7 V 

If the sums over m',m" are truncated at m max , this requires £}(m m ax^max) operations. For 
elliptical beams, h\ m falls off rapidly with m, becoming negligible for m > 6. 

To quantify how beam asymmetry distorts the power spectrum, we compute the ratio 

a ^5^E G ^' d ' (B31) 

where C ; fid is a fiducial power spectrum, is the symmetrized beam transfer function noted 
above, and G\ v is the coupling matrix appropriate to cross-power spectrum i (this index 
notation is defined in §7.11) . The matrices Gw were computed usin g the three-year sk y 



coverage with beam transforms derived from the hybrid beam maps (IJarosik et al.l 120061 ). 
and with Z max = 1500, m max = 16. The asymmetry corrections for each of the Q-, V-, and 
W-band auto- and cross-power spectra are shown in Figure [231 The QQ spectrum show 
the largest effect (6.3% at I = 600), because the Q-band beams are relatively elliptical. For 
example, the FWHM along the major and minor axes of the Q2A beam are 0?60 and 0?42, 
respectively. The rise of a.\ at high-Z indicates that the symmetric window underestimates the 
power in the beam, and thus overestimates the power spectrum. Because of the magnitude 
of this effect in the Q-band data, and because of the large foreground correction needed at 
low-/, Q-band was excluded from the final combined power spectrum. The corrections to the 
V- and W-band spectra are all < 1% for I < 1000. 
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C. MAXIMUM LIKELIHOOD NOTES 

C.l. The General Case Without Noise 

The log of the likelihood function has the form 

-21nL(d|m) =x 2 + lndetC, (CI) 

where d is the data, which is assumed to be drawn from a multivariate Gaussian distribution, 
m is a set of model parameters, and C is the covariance matrix of the data. In the limit of 
no noise and full sky coverage, it is convenient to represent the data in terms of the spherical 
harmonic coefficients, d = a = ai m , in which case 

C = diag[C ,C 1 ,...,C 2 ,...] (C2) 

and 

Wlr. "" 



X 2 = a r .C- 1 -a = ^^ (C3) 

lm 

where Ci is the power spectrum that encodes the variance of the a\ m coefficients 

( a lmdl m ) = Ci8i V 8 mm i. (C4) 

It is straightforward to show that the maximum likelihood estimate of the power spectrum, 
C Z ML , is given by 



iML \^ \ a lm 



2 



m=— I 

To set up the case of incomplete sky coverage, consider a new representation of the data 
which is related to a by a simple linear transformation 

d = L • a, (C6) 

where L is an Nd x (/ max + l) 2 matrix, Nd is the number of points in d and / max is the 
maximum spherical harmonic order being analyzed. For example, if d = t is a sky map, the 
matrix L = Y is the matrix of spherical harmonics evaluated at each pixel %. (We need not 
assume that t covers the full sky.) 

The likelihood of d, given the model spectrum, C/, is 

-21nL(d|m) =d T -CT 1 -d + lndetC, (C7) 

where C is is the covariance matrix of d which has the form 

C = <dd T ) = L • <aa T ) • L T = L • C • L T , (C8) 
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with 

C _1 = (L T ) -1 • C _1 • L -1 . (C9) 

The determinant of C factors as 

IndetC = ln(det L det C det L T ) = In det C + ln(det L) 2 , (CIO) 
and since C is diagonal and L is independent of Ci, this reduces to 

IndetC = J^ln Q + const. (Cll) 

/ m 

Thus we can rewrite the likelihood as 

- 2 In L(d|m) = (L _1 d) T ■ C" 1 • (L _1 d) + In det C + const. (C12) 

This form is a trivial recasting of the original form of L in terms of L M. By the same 
argument that led to equation flC5j) we have the following expression for the maximum 
likelihood estimate of Ci 

c i - 2^ 2Z + 1 • {LU} 

m=—l 

So, in the limit of no noise, a complete maximum likelihood estimate of the low I spectrum 
requires only a single matrix inversion to evaluate L _1 . This enables very rapid computation 
of the likelihood and it isolates the source of systematic errors to the form of L and its 
inverse. 



C.2. The Pseudo-a; m Case 

One representation of the data are the pseudo-a; m coefficients, di m , which are obtained 
by transforming the data on the cut sky 

aim = ^2^i Y i* m (ni) ti, (C14) 

i 

where the sum is over all pixels i, Wi is the weight per pixel, which includes any mask applied, 
Y* m {fii) is the spherical harmonic in the direction of pixel i, and ti is the temperature in pixel 
i. If we expand the temperature in spherical harmonics, we can express the a\ m in terms of 
the true ai m coefficients as 

m'Yl'm'(p>i) — ^ ^ M \lm)(lm)' ai'm' ; 

(C15) 

i I'm' I'm' 
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where the coupling matrix is defined as 

Afym)(im)' = ^ Wj Y* m (nj)Yi /m/ (hi) , (C16) 

i 

which depends only on the weight array, w. (If one expands w in spherical harmonics, the 
coupling matrix can be expressed in terms of the Wigner 3-j symbols, but this form is not 
especially useful in this context.) We express the pseudo-a; m in matrix notation as a = M a, 
which has the same form as equation (IC6I) . It follows that the maximum likelihood solution 
for Ci is then 

cr = e l(M 2/ !f ) ; m|2 - (ci7) 

m=— I 

In practice, the coupling matrix depends on both Z max and l' max , the harmonic orders to 
which the true and pseudo-a/ m are analyzed, respectively. In general M need not be square. 
It is instructive to look at the structure of M using SVD 

a = M • a = uwv T ■ a. (C18) 

Since u and v are orthogonal, the values of w encode the level of deprojection that occurs 
when solving for a = M _1 • a. Specifically, since 

1VT 1 = vw-'u T , (C19) 

it follows that pseudo-c^ m modes corresponding to small values of w will be greatly amplified 
in the recovery of the true a\ m . Thus small errors in the a~i m can be also be greatly magnified. 
In the case of full sky, w = 1, while for the Kp2 cut sky and / max = l' max = 10, we find 
0.05 < w < 1. Note also that any error in the a\ m will bias the estimate of C\ obtained with 
this method, since the estimator for C ; ML is quadratic. 



C.3. The Pixel-Space Case 

As noted above, the data may be represented as a (cut) sky map, t, which can be 
expanded in spherical harmonics, t = Y • a. Here Y is an iVpi x x (Z max + l) 2 matrix of 
spherical harmonic functions evaluated at each pixel % that survives the sky cut. 

The likelihood of t, given the model spectrum, Ci, is 

-21nL(t|m) = t T • CT 1 • t + lndet C, (C20) 
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where C is is the covariance matrix of t which has the form 

C = <tt T ) = Y ■ <aa T ) ■ Y T = Y ■ C ■ Y T . (C21) 

This may be expressed in the more familiar form of the 2-point correlation function, C(%) 

2/ + 1 

Cij = C(%) = ^ / ——C l Pi(coBe ij ) = J2 c iYL(n i )Y lm (h j ). (C22) 

l lm 

where the last equality follows from the addition theorem for spherical harmonics. The 
inverse of C has the form 

C- 1 = (Y T y 1 -C- 1 -Y~\ (C23) 
and the determinant of C factors as 

In det C = ln(det Y det C det Y T ) = In det C + ln(det Y) 2 = In det C + const. (C24) 

Thus 

-21nL(t|m) = (Y _1 t) T • C _1 • (Y~H) + In det C + const, (C25) 

and 

Cl - 2L, 2Z + 1 • [L b} 

m=—l 

What is the nature of Y _1 ? In the limit of full sky coverage the spherical harmonics 
form an orthonormal basis, 

or Y T • Y I. where I is the identity matrix. In this limit Y is orthogonal, Y _1 = Y T , and 
Y _1 t = Y T t is just the spherical harmonic transform of t. 

More generally, in the presence of a sky cut, we can relate the pixel space formalism to 
the pseudo-a; m formalism as follows. The pseudo-a im coefficients may be written as 

a = Y T ■ t = Y T Y • a = M • a, (C28) 

where M = Y T Y is the coupling matrix defined in equation fl016[) . Note that we are now 
implicitly assuming that the sky mask, w, is sharp in the sense that it consists of only l's 
or O's. The case of an apodized mask needs to be studied further. It follows that 

a = 1VT 1 • a = (Y T Y) _1 Y T • t = Y" 1 (Y T )" 1 Y T • t = Y" 1 • t, (C29) 
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thus Y -1 = M _1 Y T (this reduces to Y _1 = Y T in the limit of full sky coverage). The 
relations in equation (IG29P establish the equivalence of equation flQ17j) and equation (IC26I) 
in the case of a sharp mask. Note also that equation (IU26P is the same expression one would 
get by least squares fitting the a lm coefficients on the cut sky (with uniform weight outside 
the cut) and summing them to get the power spectrum. 

Since Y -1 = M _1 Y T , the essential numerical features of the maximum likelihood esti- 
mates are encoded in the coupling matrix M. Specifically, one must ensure that M _1 has 
sufficient range to adequately capture the harmonic content in the pseudo-a; m data. 
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Table 1. Differencing Assembly (DA) Properties 
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(mm) 


(GHz) 




(°) 


(mK) 


(mK) 


(GHz) 


(GHz) 


(GHz) 


Kl 


13.17 


22.77 


1.0135 


0.807 


1.439 


1.455 


22.47 


22.52 


22.78 


Kal 


9.079 


33.02 


1.0285 


0.624 


1.464 


1.483 


32.71 


32.76 


33.02 


Ql 


7.342 


40.83 


1.0440 


0.480 


2.245 


2.269 


40.47 


40.53 


40.85 


Q2 


7.382 


40.61 


1.0435 


0.475 


2.135 


2.156 


40.27 


40.32 


40.62 


VI 


4.974 


60.27 


1.0980 


0.324 


3.304 


3.330 


59.65 


59.74 


60.29 


V2 


4.895 


61.24 


1.1010 


0.328 


2.946 


2.970 


60.60 


60.70 


61.27 


Wl 


3.207 


93.49 


1.2480 


0.213 


5.883 


5.918 


92.68 


92.82 


93.59 


W2 


3.191 


93.96 


1.2505 


0.196 


6.532 


6.571 


93.34 


93.44 


94.03 


W3 


3.226 


92.92 


1.2445 


0.196 


6.885 


6.925 


92.34 


92.44 


92.98 


W4 


3.197 


93.76 


1.2495 


0.210 


6.744 


6.780 


93.04 


93.17 


93.84 



a Effective wavelength and frequency for a thermodynamic spectrum. 

b Conversion from antenna temperature to CMB thermodynamic t emperature; 

AT = g(v)AT A ; g{v) = [{e x - l) 2 /xV]; x = hu/kT ; T = 2.725 K jMather et alill999h . 

c Full-width-at-half-maximum from radial profile of A- and B-side average beams. Note: 
beams are not Gaussian. 

d Noise per observation for resolution 9 and 10 /, Q, & U maps, to ~0.1% uncertainty. 

e Effective frequency for synchrotron (s), free-free (ff), and dust (d) emission, assuming 
spectral indices of (3 = —2.9, —2.1, +2.0, respectively, in antenna temperature units. 
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Table 2. Data Flagging Summary 



Category 


K-band 


Ka-band 


Q-band 


V-band 


W-band 


Lost or rejected data 












Lost a (%) 


0.43 


0.43 


0.43 


0.43 


0.43 


Thermal disturbance 1 ^ (%) 


0.51 


0.51 


0.51 


0.51 


0.51 


Gain/baseline step (%) 


0.02 


0.04 


0.05 


0.00 


0.06 



Total lost or rejected (%) 0.96 0.98 0.99 0.94 1.00 

Data not used in maps 
Planet in beam (%) 0.11 0.11 0.11 0.11 0.11 

a Primarily due to one solar storm induced safehold. 
b Primarily due to station-keeping maneuvers at L 2 . 



Table 3. Change in low-/ Power 



Band 


/ = l a 


I = 2 b 


/ = 3 b 




(jjlK) 




GuK 2 ) 


K 


10.1 


38.8 


9.1 


Ka 


7.3 


2.7 


3.4 


Q 


6.1 


7.1 


12.3 


V 


5.1 


7.1 


2.2 


w 


7.0 


5.8 


1.5 



H = 1 - Amplitude in the differ- 
ence map, in /iK. 

h l > 1 - Power in the difference 
map, 1(1 + 1) C//27T, in [iK 2 . 
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Table 4. Galactic Template Fit Coefficients 



DA 






b 3 b 




h s 


Ql 


0.243 


1.020 


0.195 


-3.15 


6.28 


Q2 


0.242 


1.046 


0.193 


-3.19 


6.28 


VI 


0.054 


0.700 


0.448 


-3.50 


6.57 


V2 


0.052 


0.678 


0.463 


-3.48 


6.57 


Wl 


0.000 


0.405 


1.226 




6.72 


W2 


0.000 


0.400 


1.240 




6.72 


W3 


0.003 


0.397 


1.207 




6.71 


W4 


0.000 


0.402 


1.234 




6.72 



a Units of fiK R _1 . Coefficients con- 
strained to follow a free-free spectrum. See 
text. 

b Dimensionless amplitude relative to FDS 
Model 8 evaluated at 94 GHz. Coefficients 
forced to follow (3d = +2.0. See text. 

c Synchrotron spectral index computed 
from equation (I24"I) . The W-band fits imply 
Pa{v K , vw) < -4- 
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Table 5. ILC Weights by Region 



Region a 


K-band 


Ka-band 


Q-band 


V-band 


W-band 





0.1559 


-0.8880 


0.0297 


2.0446 


-0.3423 


1 


-0.0862 


-0.4737 


0.7809 


0.7631 


0.0159 


2 


0.0358 


-0.4543 


-0.1173 


1.7245 


-0.1887 


3 


-0.0807 


0.0230 


-0.3483 


1.3943 


0.0118 


4 


-0.0781 


0.0816 


-0.3991 


0.9667 


0.4289 


5 


0.1839 


-0.7466 


-0.3923 


2.4184 


-0.4635 


6 


-0.0910 


0.1644 


-0.4983 


0.9821 


0.4428 


7 


0.0718 


-0.4792 


-0.2503 


1.9406 


-0.2829 


8 


0.1829 


-0.5618 


-0.8002 


2.8464 


-0.6674 


9 


-0.0250 


-0.3195 


-0.0728 


1.4570 


-0.0397 


10 


0.1740 


-0.9532 


0.0073 


2.7037 


-0.9318 


11 


0.2412 


-1.0328 


-0.2142 


2.5579 


-0.5521 



a See Figure [8] for the region definitions. See §5.21 for notes on 
the handling of regions and 1. 
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Table 6. Low-/ Multipole Moments (/xK) 





Z = l b 


Z = 2 C 


/ = 3 d 


a>i+3 






-11.24 ± 1.48 






-14.41 ±3.13 


22.03 ±0.37 




-239.3 ± 1.2 


-0.05 ±2.62 


-13.05 ±1.55 




2505.0 ± 11.1 


11.48 ±3.59 


-5.99 ±0.46 


ai-i 


-2223.6 ± 12.5 


4.86 ±2.62 


2.45 ±0.83 






-18.80 ±2.88 


0.70 ±0.26 








33.46 ±1.35 



a Phase convention: the coefficients of the real valued har- 
monics reported here, di m , are related to the complex coeffi- 
cients by di m = v2 Imfl (m for Z < 0, di m = aj m for I = 0, and 
a im = v^Reaim for Z > 0. Also J2 m |a Zm | 2 = £) m |a /m | 2 . 

b The dipole components in Galactic rectilinear coordi- 
nates are (x,y,z) = (5i + i,ai_i,aio); in Galactic polar coor- 
dinates the components are (d,l,b) = (3.358 ± 0.017 mK, 
263.86 ± 0.04°, 48.24 ± 0.10°). The uncertainty in the 
dipole magnitude is dominated by the 0.5% calibration un- 
certainty. 

c Uncertainties fro m a quadrature sum of calibration error 
(IJarosik et al.ll2006l ) and Galaxy subtraction error ( §5.21) . 

Uncertainties from Galaxy subtraction ( §5.2p . 



-83- 



Table 7. Maximum Likelihood Spectrum (1 = 2 — 10) 



I 


AT; 2 


-95% 


-68% 


+68% 


+95% 


2 


236 


54 


99 


796 


3827 


3 


1053 


347 


552 


2398 


6036 


4 


752 


285 


434 


1472 


3131 


5 


1582 


664 


977 


2811 


5338 


6 


601 


272 


388 


1000 


1747 


7 


1345 


639 


904 


2126 


3591 


8 


655 


330 


453 


998 


1589 


9 


621 


324 


438 


922 


1407 


10 


755 


409 


545 


1088 


1626 



Table 8. WMAP Power Spectrum Peak and Trough Data 



Quantity 


I 


AT 2 








First peak 


220.8 + 0.7 


5624 ± 30 


First trough 


412.4+1.9 


1716 + 28 


Second peak 


530.9 + 3.8 


2485 ± 44 


Second trough 


675.2 ± 11.1 


1688 + 81 
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Table 9. WMAP Source Catalog 



RA [hms] 


Dec [dm] 


ID 


K [Jy] 


Ka [Jy] 


Q [Jy] 


V [Jy] 


W [Jy] 


a 


5 GHz ID 


00 


06 


05 


—06 22 


060 


94 + nns 


2.2 ± 0.1 


2.3 ± 0.2 


2 


3 ±0.3 




—0.0 ± 0.2 


PMN T0006-0623 


00 


12 


56 


—39 53 

*JZj UO 


202 


i 4 4. o 06 


1 3 + 1 

± .O —1— U. X 


10 + 01 

X .U —1— U . J. 


1 


5 ± 0.3 




_0 2 + 03 

yj.Li _i_ u.o 


PMN 10013-3954 


00 


19 


32 


26 02 




i i 4. o 08 

J. . J. _l_ U. UO 


1.1 ± 0.2 


q 4- o 2 





8 ± 0.2 




_0 Q 4- o 5 

U.O —1— U.O 


CR6 10019+2602 


00 


19 


33 


20 19 




_L . X _l_ U.Ul 


1.1 ± 0.1 


9 + 01 

U . Zj —1— U . X 


1 


9 ± 0.4 




0.2 ± 0.4 


r!R6 T001 9+2021 

v i 1 JU O UU X V\^Li\JLi X 


00 


25 


26 


— 9fi 04 

Li\J U4: 




x.u zn u.uo 


08 + 01 
u.o zn u. i 










—0 A + 1 

U .4; HZ -L 


PMN 1009^ 9fi09 

JT 1VX1N <J\J\J DO Li\j\J Li 


00 


2G 


03 


OU X X 




i n _l n no 

X .U —1— u.uo 


10 + 01 

± .U —1— U. X 


1.2 ± 0.1 


1 


2 ±0.3 




0^4-04 

U.O —1— U . 4: 


PMN 10026-351 2 

X 1VX1 N tJUUz_,UOUXz_, 


00 


29 


35 


05 55 

UU UU 




J. .U —1— U.Ul 


1.4 ± 0.1 


1.2 ± 0.1 








4 4-04 


CP>6 I0029+0554R 

vjtxju <j uu^ji? r uoo4:Xj 


00 


42 


57 


52 08 

OLi uO 




U . iJ —1— U.UU 


fi + 08 

U .U _l— u .uo 










_0 q 4- o 8 

U . Z3 —1— U.O 


CR6 T0043+5203 


00 


47 


20 


—25 1 5 

LiO XU 


062 


X . X _l_ U.Ul 




1.1 ± 0.1 








_0 o + 05 

U .U —1— U.O 


PMN 10047-2517 

X 1VX1>I JUUI 1 Li 'J X 1 


00 


49 


50 


— ^7 40 

U I 4:U 


1 79 


i q _)_ 07 

J. .O —1— U.Ul 


1.2 ± 0.1 


i o + 09 

J. .U —1— U. Ut/ 


1 


2 ±0.2 


q 4- o 4 

U. iJ — 1— U.T: 


_0 q 4- o 3 

U.O —1— U.O 


PMN 10050-5738 

X 1VX1N OXJXJOXJ O ( OO 


00 


50 45 


—06 47 




i n + o 07 

J. .U —1— U.Ul 


1.2 ± 0.1 


8 + 01 

U.O —1— U . J- 


1 


1±0.2 


i q 4. fi 

X . ij _l_ u. u 


1+03 

U. X _l_ U.O 


PMN 10051-0650 

X 1VX1 N OXJXJ'J X UUUU 


00 


50 


47 


—42 23 




i o + 05 


1.2 ± 0.07 


q 4- o 09 

U.t/ —1— U.Ut/ 





8 ±0.2 




—0.3 ± 0.3 


PMN 10051-4226 


00 


51 


06 


—09 27 

XJiJ Li 1 


077 

U i I 


i n + o 08 

X -U —1— u.uo 


9 + 1 

U. tJ —1— U. X 


8 + 01 

U.O —1— U. X 








_0 4 + 06 

VJ.'-t _l_ u. u 


PMN 10050-0928 

X 1VX1N O UUUU \JiJ Liij 


01 


0G 


42 


—40 35 

rtKJ OU 


1 71 

X 1 X 


X .O _l— U.UU 


o o 4- 09 


18 + 01 

X .O —1— U. X 


1 


4 ±0.2 


X . O _l_ U.O 


—0.1 ± 0.2 


PMN T01 06-4034 

X 1VX1N (JUXUU tiUOt: 


01 


08 


26 


1 3 90 

X O Li\J 


07Q 
u i y 


1 .4: ZC U.Ul 


11+09 

X . _L ZIZ U. Li 


1 + 09 

X .U ZIZ U . 








_n 7 + 07 
u . ( zn u. / 


HRfi T010S-H^1Q 
vjtxju u u iuotioi y 


01 


08 


40 


01 35 


081 


2 5 + 08 

Li.tJ _1_ U.UO 


2.5 ± 0.1 


2.3 ± 0.2 


2 


4 ±0.3 




—0.1 ± 0.2 


CR6 101 08+01 35 

V. I 1 )\> o U X UO | U X o o 


01 


15 


11 


—01 2fi 




1 A _L A Q7 
X .U —1— U.Ul 


1.2 ± 0.1 


s 4- o 2 





5 ±0.2 




1+06 

U. X _l_ u .u 


PMN T01 1 ^-01 27 

X 1VX1>I tJUXXU \J i- Li 1 


01 


1G 


21 


— 11 37 

J- J- Ol 




i qj-0 117 

X .O _l_ U.Ul 


1.1 ±0.1 


1.2 ± 0.2 


1 


4 ±0.3 




_0 1+04 

U. X _l_ U.I 


PMN 101 1 6-1 1 36 

X 1VX1N (JUXXU XXOU 


01 


21 


44 


1 1 50 

X X U U 




i ij-0 07 

_L . Li _l_ U.Ul 


1.1 ±0.1 


1.2 ± 0.2 








_0 1+05 

U. X _l_ U.O 


HR6 101 21 +1 1 49 

v I 1 )\) JUlZil \^ 114:^ 


01 


25 


20 


—00 10 

UU XU 


086 
uou 


i 1 + 007 

_L . Li _l_ U.Ul 


13 + 01 

J. .O —1— U. X 


13 + 01 

X .O _l_ U. X 


1 


2 ±0.2 




1+03 

U. X _l_ U.O 


PMN 10125-0005 

X 1VX1>I O U X Li' j \J\j\J'j 


01 


32 


35 


— 16 53 


097 


1.7 ± 0.07 


1.8 ± 0.1 


1.9 ± 0.1 








0.2 ± 0.3 


PMN 10132-1654 


01 


37 


01 


47 52 

j: 1 U Li 


080 

UOU 


4. o + o f)7 

^±.Li _l_ U.Ul 


4.4 ±0.1 


4.2 ± 0.1 


3.7±0.2 


2.2 ± 0.4 


_0 1+01 

U. X _l_ U. X 


HR6 101 36+47^)1 

v. I \ >\) JUIUUTt' ^1 


01 


37 


35 


— 94 98 

Zl*± ZjO 




i 04-007 

X .U ZIZ U.Ul 


11+01 

j. . J. zn u. -L 


17 + 01 

j. . i zn vj. l 


1 


4 ±0.2 


i o + A 


rc _l o q 

U.O Z1Z U.O 


PMN T01 XI 94^0 


01 


52 


34 


99 08 
ZZ Uo 




19-1-01 
1 .Z zt U. 1 


1.0 It u.z 


19-1-09 
1 .Z it U.Z 


1 


4 ±0.3 




9 -I- A 
U.Z zt U.4 


PRfi T01 ^9-J_990^ 
udD JU10Z-|-ZZUO 


02 


04 


51 


15 13 


092 


1.4 ±0.09 


1.6 ±0.2 


1.5 ±0.1 


1 


7 ±0.4 




0.1 + 0.4 


GB6 J0204+1514 


02 


04 


58 


32 13 


085 


1.5 ±0.09 


1.4 ±0.1 


1.2 ±0.2 





7 ±0.3 




-0.4 ±0.4 


GB6 J0205+3212 


02 


10 


51 


-51 00 


158 


2.6 ±0.06 


2.8 ±0.1 


2.8 ± 0.1 


2 


6 ±0.2 


3.0 ± 1 


0.1 ±0.1 


PMN J0210-5101 


02 


20 


54 


35 58 




1.1 ±0.08 


1.1 ±0.1 


1.0±0.2 








-0.1 ±0.6 


GB6 J0221+3556 


02 


22 


45 


-34 41 


137 


0.9 ±0.04 


1.0 ±0.06 







8 ±0.2 




0.0 ±0.4 


PMN J0222-3441 


02 


23 


16 


43 04 


084 


1.8 ±0.07 


1.4 ±0.1 


1.4±0.2 


1 


2 ±0.3 


1.7 ±0.6 


-0.4 ±0.3 


GB6 J0223+4259 


02 


31 


40 


13 18 




1.4 ±0.08 


1.2 ±0.1 


1.3 ±0.1 








-0.2 ±0.4 


GB6 J0231+1323 


02 


37 


58 


28 48 


093 


3.5 ±0.08 


3.1 ±0.1 


3.5 ±0.2 


2 


9 ±0.3 


1.8 ±0.6 


-0.1 ±0.2 


GB6 J0237+2848 


02 


38 


53 


16 36 




1.3 ±0.1 


1.4 ±0.2 


1.4±0.2 


1 


7 ±0.3 




0.2 ±0.4 


GB6 J0238+1637 


02 


53 


31 


-54 41 


155 


2.5 ±0.06 


2.8 ±0.09 


2.8 ±0.1 


2 


3 ±0.2 


2.5 ±0.3 


0.1 ±0.1 


PMN J0253-5441 


03 


03 


31 


-62 12 


162 


1.4 ±0.07 


1.4 ±0.1 


1.4 ±0.1 


1 


4 ±0.2 




0.0 ±0.2 


PMN J0303-6211 


03 


08 


25 


04 04 


102 


1.3 ±0.08 


1.3 ±0.1 


1.4 ±0.1 


0.9 ±0.2 




-0.0 ±0.4 


GB6 J0308+0406 
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Table 9 — Continued 



RA [hms] 


Dec [dm] 


ID 


K [Jy] 


Ka [Jy] 


Q [Jy] 


V[Jy] 


W [Jy] 


a 


5 GHz ID 




03 


09 


47 


-61 03 


160 


1.0 ±0.07 


1.3 ±0.1 


0.9 ±0.08 


1.0 ±0.2 




-0.0 ±0.4 


PMN J0309-6058 




03 


12 


08 


-76 47 


174 


1.1 ±0.06 


1.2 ±0.08 


1.1 ±0.09 


1.1 ±0.2 


1.2 ±0.4 


0.0 ±0.3 


PMN J0311-7651 




03 


19 


45 


41 31 


094 


10.8 ±0.07 


8.4 ±0.1 


7.0 ±0.1 


4.8 ±0.3 


3.0 ±0.5 


-0.7 ±0.06 


GB6 J0319+4130 




03 


22 


25 


-37 11 


138 


18.5 ±3.6 


12.6 ± 1.9 


10.5 ±2.5 


8.3 ±2.1 




-0.8 ±0.2 


Uy 0320-37 


b 


03 


25 


11 


22 24 




1.0 ±0.1 


1.2 ±0.2 


1.1 ±0.2 






0.3 ±0.7 


GB6 J0325+2223 




03 


29 


46 


-23 53 


123 


1.1 ±0.05 


1.3 ±0.09 


1.2 ±0.1 


1.0 ±0.2 




0.1 ±0.3 


PMN J0329-2357 




03 


34 


18 


-40 07 


146 


1.6 ±0.06 


1.7 ±0.09 


1.6 ±0.1 


1.5 ±0.2 




0.0 ±0.2 


PMN J0334-4008 




03 


3G 


49 


-12 58 




1.0 ±0.07 




1.2 ±0.1 






0.2 ±0.5 


PMN J0336-1302 




03 


39 


24 


-01 44 


106 


2.6 ±0.09 


2.6 ±0.1 


2.7 ±0.2 


1.5 ±0.3 


3.4±0.5 


0.0 ±0.2 


PMN J0339-0146 




03 


40 


27 


-21 21 




1.1 ±0.07 


1.1 ±0.1 


1.2 ±0.1 


1.4 ±0.2 




0.1 ±0.3 


PMN J0340-2119 




03 


48 


51 


-27 47 


129 


1.1 ±0.04 


1.0 ±0.07 


0.8 ±0.09 


1.4 ±0.2 




-0.2 ±0.3 


PMN J0348-2749 




03 


58 


43 


10 29 




1.2 ±0.1 


0.8 ±0.2 








-1.0± 1 


GB6 J0358+1026 




04 


03 


57 


-36 04 


136 


4.0 ±0.06 


4.6 ±0.1 


4.5 ±0.1 


4.3 ±0.2 


3.9 ±0.4 


0.1 ±0.08 


PMN J0403-3605 




04 


05 


37 


-13 05 


114 


2.0 ±0.07 


2.1 ±0.1 


1.8 ±0.1 


1.6 ±0.3 




-0.2 ±0.2 


PMN J0405-1308 




04 


07 


01 


-38 24 


141 


1.1 ±0.06 


1.1 ±0.1 


1.0 ±0.1 






-0.2 ±0.5 


PMN J0406-3826 




04 


08 


49 


-75 03 




0.8 ±0.05 


0.4 ±0.07 


0.6 ±0.2 






-1.2 ±0.8 


PMN J0408-7507 




04 


11 


24 


76 55 


082 


1.0 ±0.06 


0.7±0.1 


0.8 ±0.1 






-0.4 ±0.6 


Uy 0403+76 




04 


23 


17 


-01 20 


110 


10.0 ±0.08 


10.4 ±0.1 


10.1 ±0.2 


9.4 ±0.3 


6.4±0.5 


-0.0 ±0.05 


PMN J0423-0120 




04 


23 


58 


02 19 




1.2 ±0.07 


1.0±0.1 








-0.4 ±0.8 


GB6 J0424+0226 




04 


24 


48 


-37 57 


140 


1.6 ±0.06 


1.2 ±0.2 


1.3 ±0.1 


1.4 ±0.2 


0.8 ±0.3 


-0.3 ±0.3 


PMN J0424-3756 




04 


24 


51 


00 36 


109 


1.6 ±0.09 


1.8 ±0.1 


1.8 ±0.2 


1.4 ±0.2 




0.0 ±0.3 


GB6 J0424+0036 




04 


28 


32 


-37 57 




2.2 ±0.2 


1.5 ±0.1 


1.3 ±0.1 


1.3 ±0.2 


1.7±0.4 


-0.5 ±0.3 


PMN J0428-3756 


a 


04 


33 


12 


05 21 


108 


2.5 ±0.08 


2.7±0.1 


2.6 ±0.2 


2.5 ±0.3 


3.4 ±0.6 


0.1 ±0.2 


GB6 J0433+0521 




04 


40 


18 


-43 33 


147 


3.0 ±0.07 


2.9 ±0.1 


2.8 ±0.1 


2.2 ±0.2 




-0.2 ±0.1 


PMN J0440-4332 




04 


42 


49 


-00 18 




1.0 ±0.08 


1.0 ±0.1 


1.2 ±0.2 






0.3 ±0.6 


PMN J0442-0017 




04 


49 


21 


-81 00 


175 


1.8 ±0.06 


1.8 ±0.09 


1.7 ±0.09 


1.7 ±0.2 


1.2 ±0.3 


-0.1 ±0.2 


PMN J0450-8100 




04 


53 


27 


-28 06 


131 


1.5 ±0.06 


1.6 ±0.09 




1.3 ±0.2 


1.0 ±0.3 


-0.1 ±0.3 


PMN J0453-2807 




04 


55 


53 


-46 17 


151 


3.7 ±0.06 


4.1 ±0.1 


4.0 ±0.1 


3.8 ±0.2 


1.8 ±0.4 


0.1 ±0.1 


PMN J0455-4616 




04 


57 


00 


-23 22 


128 


2.7 ±0.05 


2.7 ±0.09 


2.6 ±0.1 


2.2 ±0.2 


1.6 ±0.4 


-0.1 ±0.1 


PMN J0457-2324 




05 


01 


16 


-01 58 




1.2 ±0.08 


1.2 ±0.1 


1.3 ±0.2 


1.1 ±0.2 


1.5 ±0.5 


0.1 ±0.3 


PMN J0501-0159 




05 


0G 


57 


-61 08 


154 


2.2 ±0.05 


1.9 ±0.08 


1.7 ±0.09 


1.2 ±0.2 


0.8 ±0.2 


-0.5 ±0.2 


PMN J0506-6109 


a 


05 


13 


44 


-20 15 




0.8 ±0.05 


1.1 ±0.2 


0.8 ±0.1 


0.5 ±0.2 




-0.1 ±0.5 


PMN J0513-2016 




05 


13 


58 


-21 55 


127 


1.0 ±0.05 


1.0 ±0.07 


1.0 ±0.09 


0.7 ±0.3 


0.9 ±0.3 


0.0 ±0.3 


PMN J0513-2159 




05 


15 


05 


-45 58 






0.4 ±0.1 


0.9 ±0.1 


1.2 ±0.3 




1.4 ± 1 


PMN J0515-4556 


a 


05 


19 


14 


-05 40 


116 


2.3 ±0.07 


2.3 ±0.3 


1.8 ±0.2 


1.1 ±0.3 




-0.4 ±0.4 




c 


05 


19 


42 


-45 46 


150 


6.5 ±0.06 


5.2 ±0.1 


4.2 ±0.1 


2.9 ±0.2 


2.0 ±0.4 


-0.7 ±0.08 


PMN J0519-4546 


a 


05 


23 


02 


-36 27 


139 


3.9 ±0.06 


3.7±0.1 


3.4 ±0.1 


3.4 ±0.2 


2.6 ±0.3 


-0.2 ±0.1 


PMN J0522-3628 
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RA [hms] 


Dec [dm] 


ID 




K [Jy] 


Ka [Jy] 


Q [Jy] 


V[Jy] 


W [Jy] 


a 


5 GHz ID 


05 25 


37 


-48 


26 




1 


.0±0.06 


1.4 ±0.09 


1.5 ±0.1 


1 


.4 ±0.3 




0.7±0.3 


PMN J0526-4830 


05 27 


32 


-12 


40 


122 


1 


.4 ±0.07 


1.7±0.1 


1.4±0.2 


1 


.2 ± 0.2 




0.0 ±0.3 


PMN J0527-1241 


05 38 


52 


-44 


05 


148 


5 


.4 ±0.06 


5.7±0.1 


5.7±0.1 


5 


.0± 0.2 


4.6 ±0.4 


0.0 ±0.06 


PMN J0538-4405 


05 40 


46 


-54 


16 


152 


1 


.4 ±0.06 


1.5 ±0.09 


1.8 ±0.1 


1 


.4 ±0.2 




0.1 ±0.2 


PMN J0540-5418 


05 42 


24 


49 


51 


095 


1 


.5 ±0.08 


1.1 ±0.1 


1.1 ±0.2 





.6 ±0.3 




-0.7 ±0.5 


GB6 J0542+4951 


05 50 


23 


-57 


32 


153 


1 


.3 ±0.05 




1.2 ±0.1 





.8 ±0.2 


0.8 ±0.4 


-0.3 ±0.3 


PMN J0550-5732 


05 55 


58 


39 


44 


100 


2 


.9 ±0.07 


2.1 ±0.1 


1.1 ±0.2 








-1.0 ±0.4 


GB6 J0555±3948 


06 00 


11 


-45 


28 







.6 ±0.07 


1.1 ±0.09 


0.9 ±0.09 





.9 ±0.2 


1.2 ±0.4 


0.4 ±0.4 


PMN J0559-4529 


06 07 


08 


67 


22 


091 


1 


,1±0.04 




0.7 ±0.1 








-0.8 ±0.7 


GB6 J0607±6720 


06 08 


44 


-22 


19 




1 


.0±0.04 


1.0 ±0.07 


1.0 ±0.1 








0.0 ±0.4 


PMN J0608-2220 


06 09 


38 


-15 


41 


126 


3 


.7 ±0.07 


3.5 ±0.1 


3.4 ±0.1 


2 


.5 ±0.3 




-0.2 ±0.1 


PMN J0609-1542 


06 23 


11 


-64 


36 







.8 ±0.04 


0.8 ±0.05 


0.9 ±0.06 





.9 ±0.1 


0.9 ±0.2 


0.1 ±0.2 


PMN J0623-6436 


06 29 


26 


-19 


57 


130 


1 


.5 ±0.05 


1.4 ±0.08 


1.6 ±0.1 


1 


.5 ±0.3 




0.0 ±0.3 


PMN J0629-1959 


06 33 


47 


-22 


17 


135 





.6 ±0.07 


0.8 ±0.07 


0.9 ±0.1 


1 


.1 ± 0.2 




0.6 ±0.5 


PMN J0633-2223 


06 35 


50 


-75 


17 


167 


4 


.7 ±0.05 


4.4 ±0.08 


4.0 ±0.09 


3 


,0± 0.2 


2.6 ±0.5 


-0.3 ±0.07 


PMN J0635-7516 


06 36 


29 


-20 


32 


134 


1 


,1±0.06 


1.1 ±0.07 


0.9 ±0.1 








-0.2 ±0.4 


PMN J0636-2041 


06 39 


46 


73 


26 


087 





.8 ±0.05 


0.8 ±0.1 


1.1 ±0.1 






1.2 ±0.3 


0.4 ±0.4 


GB6 J0639+7324 


06 46 


30 


44 


49 


099 


3 


,1±0.08 


2.5 ±0.1 


2.2 ±0.1 


1 


.8 ±0.3 




-0.6 ±0.2 


GB6 J0646+4451 


07 21 


00 


04 


01 







.9 ±0.08 




0.7 ±0.1 








-0.5 ±0.7 


GB6 J0721+0406 


07 21 


54 


71 


22 




1 


.8 ±0.06 


2.0 ±0.09 


2.2 ±0.1 


1 


.9 ±0.2 


1.7 ±0.2 


0.2 ±0.1 


GB6 J0721+7120 


07 25 


49 


-00 


49 




0.9 ±0.1 


1.0 ±0.2 


1.2 ±0.1 





.8 ±0.3 




0.3 ±0.6 


PMN J0725-0054 


07 27 


05 


67 


48 







.6 ±0.05 


0.5 ±0.08 


0.7 ±0.1 





.7 ±0.3 




0.0 ±0.6 


GB6 J0728±6748 


07 34 


15 


50 


21 




1 


,1±0.07 


1.3 ±0.1 


1.2 ±0.1 


1 


.1 ± 0.2 




0.1 ±0.3 


GB6 J0733±5022 


07 38 


15 


17 


44 


113 


1 


.4 ±0.08 


1.4 ±0.1 


1.3 ±0.1 


1 


.3 ±0.3 




-0.1 ±0.3 


GB6 J0738+1742 


07 39 


14 


01 


36 


124 


2 


.0±0.07 


2.1 ±0.1 


2.5 ±0.1 


2 


.7 ±0.3 


3.1 ±0.7 


0.3 ±0.2 


GB6 J0739+0136 


07 41 


13 


31 


11 


107 


1 


.3 ±0.07 




0.9 ±0.2 


1 


.3 ±0.4 




-0.3 ±0.5 


GB6 J0741+3112 


07 43 


38 


-67 


27 


161 


1 


.3 ±0.05 


0.8 ±0.08 


0.7 ±0.09 


1 


,1±0.3 




-0.8 ±0.4 


PMN J0743-6726 


07 45 


27 


10 


12 


118 


1 


.1±0.09 


1.0 ±0.1 


0.8 ±0.2 








-0.5 ±0.6 


GB6 J0745+1011 


07 46 


02 


-00 


45 




1 


.2 ±0.08 


1.1 ±0.1 


0.8 ±0.1 








-0.6 ±0.6 


PMN J0745-0044 


07 50 


53 


12 


30 


117 


2 


.6 ±0.08 


2.5 ±0.1 


2.6 ±0.2 


1 


.9 ±0.2 




-0.1 ±0.2 


GB6 J0750+1231 


07 53 


31 


53 


54 




1 


,0±0.08 


0.8 ±0.1 


1.0 ±0.3 


1 


,0±0.2 




-0.1 ±0.5 


GB6 J0753±5353 


07 57 


00 


09 


57 


120 


1 


.5 ±0.09 


1.6 ±0.1 


1.6 ±0.2 


1 


.3 ±0.4 


1.2 ±0.6 


0.1 ±0.3 


GB6 J0757±0956 


08 08 


18 


-07 


50 


133 


1 


.4 ±0.07 


1.5 ±0.1 


1.4 ±0.1 


1 


.9 ±0.2 


1.2 ±0.4 


0.1 ±0.2 


PMN J0808-0751 


08 13 


30 


48 


17 




1 


,0±0.08 


1.0 ±0.1 




1 


,1±0.3 




0.1 ±0.6 


GB6 J0813±4813 


08 16 


23 


-24 


25 


145 


1 


.0±0.06 


1.1 ±0.08 


1.2 ±0.1 





.9 ±0.2 




0.2 ±0.3 


PMN J0816-2421 


08 24 


56 


39 


14 




1 


,1±0.09 




1.0 ±0.2 


1 


,1±0.4 




-0.1 ±0.5 


GB6 J0824+3916 


08 25 


48 


03 


11 


125 


1 


.4 ±0.08 


1.6 ±0.1 


1.5 ±0.2 


1 


.6 ±0.2 




0.2 ±0.3 


GB6 J0825±0309 
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RA [hms] 


Dec [dm] 


ID 


K [Jy] 


Ka [Jy] 


Q [Jy] 


V[Jy] 


W [Jy] 


a 


5 GHz ID 


08 


30 


58 


24 10 


112 


1.5 ±0.09 


1.5 ±0.1 


1.6 ±0.2 


2.0 ±0.3 




0.2 ±0.3 


GB6 J0830+2410 


08 


36 


46 


-20 15 


144 


2.8 ±0.07 


2.5 ±0.1 


2.3 ±0.1 


2.1 ±0.2 


1.1 ±0.5 


-0.3 ±0.2 


PMN J0836-2017 


08 


38 


12 


58 20 






1.1 ±0.09 


1.0 ±0.1 






-0.6 ±2 


GB6 J0837+5825 


08 


40 


38 


13 12 


121 


1.8 ±0.08 


2.2 ±0.1 


1.8 ±0.2 


1.3 ±0.3 




-0.0 ±0.3 


GB6 J0840+1312 


08 


41 


23 


70 53 


089 


1.7 ±0.06 


1.8 ±0.1 


1.8 ±0.1 


1.5 ±0.2 




0.0 ±0.2 


GB6 J0841+7053 


08 


•54 


46 


20 05 


115 


3.8 ±0.09 


4.5 ±0.1 


4.2 ±0.2 


4.2 ±0.3 


3.6 ±0.5 


0.1 ±0.1 


GB6 J0854±2006 


09 


02 


16 


-14 14 




1.1 ±0.07 


1.2 ±0.1 


1.3 ±0.1 




1.8 ±0.4 


0.3 ±0.3 


PMN J0902-1415 


09 


07 


54 


-20 17 




1.2 ±0.06 


1.0 ±0.1 








-0.5 ±0.7 


PMN J0907-2026 


09 


09 


17 


01 18 


132 


1.9 ±0.08 


1.9 ±0.1 


1.8 ±0.2 


1.6 ±0.3 




-0.1 ±0.3 


GB6 J0909+0121 


09 


09 


56 


42 52 




0.9 ±0.1 


1.1 ±0.2 


1.2 ±0.2 


1.0 ±0.2 




0.2 ±0.6 


GB6 J0909+4253 


09 


14 


42 


02 49 




1.3 ±0.07 


1.6 ±0.1 


1.4 ±0.2 




1.6 ±0.6 


0.2 ±0.4 


GB6 J0914+0245 


09 


18 


11 


-12 03 


143 


2.0 ±0.08 


1.1 ±0.2 


1.0 ±0.2 






-1.4 ±0.7 


PMN J0918-1205 


09 


20 


40 


44 40 




1.2 ±0.09 


1.1 ±0.2 


1.3 ±0.1 


1.4 ±0.4 




0.2 ±0.4 


GB6 J0920+4441 


09 


21 


07 


62 16 




0.9 ±0.06 




0.9 ±0.2 




1.3 ±0.6 


0.2 ±0.6 


GB6 J0921+6215 


09 


21 


40 


-26 19 




1.2 ±0.08 


1.3 ±0.1 


1.2 ±0.1 


1.2 ±0.2 


1.8 ±0.7 


0.0 ±0.3 


PMN J0921-2618 


09 


27 


05 


39 01 


105 


6.5 ±0.09 


5.5 ±0.1 


5.2 ±0.1 


4.4 ±0.3 


2.4 ±0.4 


-0.4 ±0.09 


GB6 J0927+3902 


09 


48 


58 


40 37 


104 


1.4 ±0.07 


1.7±0.1 


1.4 ±0.1 


1.1 ±0.2 




0.0 ±0.3 


GB6 J0948±4039 


09 


•55 


53 


69 36 


088 


1.3 ±0.07 


1.1 ±0.09 


1.1 ±0.1 


1.0 ±0.2 


0.8 ±0.4 


-0.3 ±0.3 


GB6 J0955±6940 


09 


•57 


24 


55 26 




0.9 ±0.08 


0.8 ±0.1 


1.0 ±0.1 






0.0 ±0.5 


GB6 J0957+5522 


09 


•58 


07 


47 21 


098 


1.5 ±0.07 


1.5 ±0.1 


1.4 ±0.1 






-0.1 ±0.3 


GB6 J0958+4725 


10 


14 


19 


23 04 


119 


1.2 ±0.08 


1.0 ±0.1 


0.9 ±0.2 






-0.4 ±0.6 


GB6 J1014+2301 


10 


18 


50 


-31 29 




0.9 ±0.06 






1.0 ±0.4 




0.1 ±0.9 


PMN J1018-3123 


10 


32 


38 


41 18 


103 


1.1 ±0.07 


1.0 ±0.1 


0.9 ±0.2 


1.5 ±0.6 




-0.1 ±0.5 


GB6 J1033+4115 


10 


37 


20 


-29 34 




1.3 ±0.07 


1.3 ±0.1 


1.1 ±0.1 


1.5 ±0.3 


1.0 ±0.4 


-0.0 ±0.3 


PMN J1037-2934 


10 


38 


42 


05 11 


142 


1.8 ±0.08 


1.8 ±0.1 


1.5 ±0.2 


1.6 ±0.3 




-0.2 ±0.3 


GB6 J1038+0512 


10 


41 


25 


06 11 




1.2 ±0.1 


1.5 ±0.2 


1.4 ±0.2 


1.1 ±0.3 




0.1 ±0.5 


GB6 J1041+0610 


10 


41 


37 


-47 40 


163 


1.1 ±0.06 


0.8 ±0.1 


0.3 ±0.1 






-1.2 ±0.8 


PMN J1041-4740 


10 


47 


24 


71 43 


083 


1.0 ±0.08 


0.9 ±0.1 


0.8 ±0.1 


0.7 ±0.3 




-0.3 ±0.5 


GB6 J1048+7143 


10 


47 


53 


-19 11 




1.1 ±0.07 


1.0 ±0.1 


0.8 ±0.2 


1.1 ±0.3 




-0.3 ±0.5 


PMN J1048-1909 


10 


•58 


28 


01 34 


149 


4.5 ±0.07 


4.5 ±0.1 


4.5 ±0.2 


4.3 ±0.3 


3.1 ± 1 


-0.0 ±0.1 


GB6 J1058+0133 


10 


•59 


28 


-80 03 


176 


1.9 ±0.06 


2.2 ±0.09 


2.2 ±0.1 


2.5 ±0.2 




0.3 ±0.2 


PMN J1058-8003 


11 


02 


03 


-44 00 




0.6 ±0.05 


0.9 ±0.08 




0.5 ±0.2 




0.5 ±0.6 


PMN J1102-4404 


11 


07 


13 


-44 46 


166 


1.6 ±0.05 


1.5 ±0.07 


1.2 ±0.09 


1.5 ±0.2 


1.0 ±0.4 


-0.2 ±0.2 


PMN J1107-4449 


11 


17 


47 


-46 35 




1.0 ±0.04 


0.7 ±0.07 


0.6 ±0.1 






-1.0 ±0.5 


PMN J1118-4634 


11 


18 


50 


12 38 




1.0 ±0.08 




1.1 ±0.2 


1.0 ±0.3 




0.1 ±0.5 


GB6 J1118+1234 


11 


27 


03 


-18 58 


159 


1.5 ±0.07 


1.4 ±0.1 


1.3 ±0.1 






-0.2 ±0.4 


PMN J1127-1857 


11 


30 


14 


-14 52 


157 


1.5 ±0.08 


1.5 ±0.1 


1.5 ±0.2 


1.3 ±0.3 




-0.1 ±0.3 


PMN J1130-1449 
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RA [hms] 


Dec [dm] 


ID 


K [Jy] 


Ka [Jy] 


Q [Jy] 


V [Jy] 


W [Jy] 


a 


5 GHz ID 


11 


30 


45 


38 


13 


101 


1.3 ±0.07 


1.2 ±0.1 


1.1 ±0.1 


1.0 ±0.2 




-0.3 ±0.4 


GB6 J1130±3815 


11 


4G 


15 


-48 


38 




0.7 ±0.05 


0.6 ±0.1 








-0.4 ± 1 


PMN J1145-4836 


11 


4G 


40 


40 


01 




0.9 ±0.07 


0.9 ±0.1 


1.1 ±0.1 




1.5 ±0.5 


0.4 ±0.4 


GB6 J1146±3958 


11 


47 


09 


-38 


11 


169 


1.9 ±0.07 


2.1 ±0.1 


1.9 ±0.1 


1.8 ±0.2 




-0.0 ±0.2 


PMN J1147-3812 


11 


49 


45 


-79 


32 




1.2 ±0.05 




0.7±0.3 






-0.9 ±2 


PMN J1147-7936 


11 


53 


16 


49 


31 


090 


2.0 ±0.05 


2.0 ±0.08 


2.2 ±0.1 


1.9 ±0.2 


1.6 ±0.3 


0.0 ±0.2 


GB6 J1153±4931 


11 


54 


46 


81 


04 


078 


1.0 ±0.07 


0.9 ±0.1 




1.2 ±0.5 




0.0 ±0.7 


Uy 1150+81 


11 


59 


37 


29 


14 


111 


2.0 ±0.07 


2.2 ±0.1 


2.0 ±0.1 


2.0 ±0.3 


1.5 ±0.6 


0.0 ±0.2 


GB6 J1159±2914 


12 


03 


21 


48 


06 




0.8 ±0.05 


0.8 ±0.07 


0.7 ±0.1 


0.6 ±0.2 




-0.3 ±0.4 


GB6 J1203±4803 


12 


08 


59 


-24 


05 


172 


1.3 ±0.07 




0.9 ±0.1 






-0.7 ±0.7 


PMN J1209-2406 


12 


15 


54 


-17 


30 


173 


1.5 ±0.08 


1.4 ±0.1 


1.2 ±0.1 






-0.4 ±0.4 


PMN J1215-1731 


12 


18 


55 


48 


32 




0.8 ±0.04 


0.8 ±0.06 


0.8 ±0.09 


0.8 ±0.2 




-0.1 ±0.4 


GB6 J1219±4830 


12 


19 


20 


05 


48 


164 


2.6 ±0.07 


2.2 ±0.1 


2.0 ±0.2 


1.5 ±0.3 




-0.5 ±0.2 


GB6 J1219+0549A 


12 


22 


54 


-83 


07 


178 


1.0 ±0.05 


1.0 ±0.1 


0.9 ±0.09 


0.9 ±0.3 




-0.1 ±0.4 


PMN J1224-8312 


12 


27 


51 


11 


24 








0.5 ±0.2 


0.5 ±0.2 




-0.0 ±4 


GB6 J1228±1124 


12 


29 


06 


02 


03 


170 


19.1 ±0.08 


17.3 ±0.1 


15.7±0.1 


13.1 ±0.3 


8.5 ±0.4 


-0.4 ±0.03 


GB6 J1229±0202 


12 


30 


50 


12 


23 


165 


19.1 ±0.07 


15.5 ±0.1 


13.1 ±0.1 


9.3 ±0.3 


5.9 ±0.5 


-0.6 ±0.03 


GB6 J1230±1223 


12 


31 


09 


13 


51 




0.2 ±0.1 






0.7 ±0.3 




1.1 ±1 


GB6 J1231±1344 


12 


4G 


52 


-25 


47 


177 


1.3 ±0.08 


1.5 ±0.2 


1.8 ±0.2 


1.1 ± 0.2 


1.3 ±0.5 


0.2 ±0.3 


PMN J1246-2547 


12 


56 


11 


-05 


47 


181 


17.1 ±0.08 


17.9 ±0.1 


18.0 ±0.1 


16.6 ±0.3 


13.1 ±0.5 


0.0 ±0.03 


PMN J1256-0547 


12 


58 


06 


-31 


58 


180 


1.5 ±0.07 


1.2 ±0.1 


1.3 ±0.2 






-0.4 ±0.4 


PMN J1257-3154 


12 


58 


55 


-22 


23 




1.0 ±0.08 


0.9 ±0.1 


1.0 ±0.2 






-0.1 ±0.7 


PMN J1258-2219 


13 


02 


49 


48 


56 








0.9 ±0.07 




1.1 ±0.4 


0.2 ±0.9 




13 


05 


54 


-49 


28 




1.1 ±0.06 


1.1 ±0.1 




0.9 ±0.2 




-0.1 ±0.5 


PMN J1305-4928 


13 


09 


25 


11 


55 




1.0 ±0.07 


0.9 ±0.1 


1.0 ±0.1 


1.2 ±0.3 




0.1 ±0.4 


GB6 J1309±1154 


13 


10 


35 


32 


21 


052 


2.9 ±0.07 


3.0 ±0.1 


2.7 ±0.1 


2.4 ±0.2 


2.1 ±0.6 


-0.1 ±0.1 


GB6 J1310±3220 


13 


1G 


05 


-33 


37 


182 


1.5 ±0.07 


1.6 ±0.1 


1.7±0.1 


1.7 ±0.2 


1.7 ±0.4 


0.1 ±0.2 


PMN J1316-3339 


13 


24 


29 


-10 


46 




0.9 ±0.09 


0.9 ±0.1 


1.0 ±0.1 


1.4 ±0.3 


1.5 ±0.6 


0.4 ±0.4 


PMN J1324-1049 


13 


27 


25 


22 


12 




1.0 ±0.07 


1.0 ±0.09 


0.8 ±0.1 


1.4 ±0.5 


1.1 ±0.4 


0.1 ±0.4 


GB6 J1327±2210 


13 


29 


15 


32 


01 


040 


0.8 ±0.06 




0.8 ±0.3 






-0.0 ± 1 


GB6 J1329±3154 


13 


30 


55 


25 


04 




1.1 ±0.07 


1.1 ±0.1 


1.1 ±0.1 






-0.0 ±0.5 


GB6 J1330±2509 


13 


31 


19 


30 


31 


026 


2.1 ±0.07 


1.8 ±0.1 


1.4 ±0.1 






-0.6 ±0.3 


GB6 J1331±3030 


13 


32 


54 


02 


00 




1.3 ±0.07 


1.3 ±0.1 


1.1 ±0.1 


1.2 ±0.3 


0.9 ±0.4 


-0.2 ±0.3 


GB6 J1332+0200 


13 


33 


33 


27 


23 




0.8 ±0.07 


0.9 ±0.1 


0.9 ±0.1 


0.7±0.2 




0.1 ±0.5 


GB6 J1333±2725 


13 


3G 


49 


-33 


58 


185 


1.9 ±0.07 


1.4±0.1 


1.1 ±0.1 


1.7 ±0.4 




-0.6 ±0.3 


PMN J1336-3358 


13 


37 


39 


-12 


57 


188 


5.9 ±0.08 


6.3 ±0.1 


6.4 ±0.1 


6.1 ±0.2 


4.1 ±0.4 


0.0 ±0.07 


PMN J1337-1257 


13 


43 


57 


66 


01 




0.8 ±0.08 


0.4 ±0.1 


0.6 ±0.1 




1.4 ±0.3 


0.2 ±0.4 


GB6 J1344+6606 
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RA [hms] 


Dec [dm] 


ID 


K [Jy] 


Ka [Jy] 


Q [Jy] 


V [Jy] 


W [Jy] 


a 


5 GHz ID 


13 


54 


51 


-10 41 


197 


1.7 ±0.08 


1.4 ±0.1 


1.6 ±0.2 


1.5 ±0.3 




-0.1 ±0.3 


PMN J1354-1041 


13 


56 


56 


19 19 


004 


1.3 ±0.07 


1.4 ±0.1 


1.2 ±0.1 


1.5 ±0.3 


0.9 ±0.4 


-0.0 ±0.3 


GB6 J1357+1919 


13 


57 


04 


76 44 




0.7 ±0.07 


0.9 ±0.2 


0.8 ±0.1 


0.8 ±0.2 


1.3 ±0.4 


0.3 ±0.4 




13 


57 


06 


-15 29 




0.9 ±0.1 




1.0 ±0.2 


1.1 ±0.3 




0.1 ±0.6 


PMN J1357-1527 


14 


08 


53 


-07 49 


203 


1.2 ±0.09 


1.0 ±0.1 


1.1 ±0.2 




0.9 ±0.4 


-0.2 ±0.5 


Uy 1406-076 


14 


11 


16 


52 16 




0.8 ±0.05 


0.5 ±0.1 








-1.2 ± 1 


GB6 J1411+5212 


14 


19 


34 


54 25 




0.9 ±0.07 


0.8 ±0.1 


0.9 ±0.1 


1.5 ±0.3 


1.6 ±0.4 


0.4 ±0.3 


GB6 J1419+5423 


14 


19 


38 


38 23 


042 


0.9 ±0.05 


1.1 ±0.08 


1.0 ±0.08 


1.1 ± 0.1 




0.1 ±0.3 


GB6 J1419+3822 


14 


20 


00 


27 02 




0.9 ±0.07 


1.0 ±0.09 


1.1 ±0.1 


1.1 ±0.3 




0.3 ±0.4 


GB6 J1419+2706 


14 


27 


31 


-33 03 


193 


1.0 ±0.07 


1.5 ±0.1 


1.7 ±0.2 


1.6 ±0.2 




0.6 ±0.3 


PMN J1427-3306 


14 


27 


50 


-42 06 


191 


3.2 ±0.07 


3.0 ±0.1 


2.9 ±0.1 


2.4 ±0.3 




-0.2 ±0.2 


PMN J1427-4206 


14 


40 


05 


49 58 










0.8 ±0.2 


1.8 ±0.8 


2.0 ±3 


GB6 J1439+4958 


14 


42 


56 


51 56 




0.9 ±0.06 


1.0 ±0.08 


0.8 ±0.09 


1.1 ±0.2 




0.1 ±0.4 


GB6 J1443+5201 


14 


58 


39 


71 40 


071 


1.4 ±0.07 


1.4 ±0.1 


1.1 ±0.1 


0.9 ±0.3 


1.0 ±0.3 


-0.4 ±0.3 


GB6 J1459+7140 


15 


04 


34 


10 29 


006 


1.8 ±0.07 


1.8 ±0.1 


1.6 ±0.1 






-0.2 ±0.3 


GB6 J1504+1029 


15 


06 


52 


-16 44 




1.4 ±0.08 


1.5 ±0.2 


1.0 ±0.1 






-0.3 ±0.5 


PMN J1507-1652 


15 


10 


32 


-05 48 




1.1 ±0.08 




1.4 ±0.2 






0.4 ±0.6 


PMN J1510-0543 


15 


12 


46 


-09 05 


207 


1.7 ±0.08 


1.8 ±0.1 


1.8 ±0.2 


1.9 ±0.3 


1.3 ±0.4 


0.1 ±0.3 


Uy 1510-08 


15 


14 


00 


-10 13 




0.9 ±0.1 


0.9 ±0.2 


1.1 ±0.2 


0.8 ±0.3 




0.1 ±0.6 


PMN J1513-1012 


15 


16 


38 


00 13 


002 


1.6 ±0.07 


1.9 ±0.1 


1.7 ± 0.1 


1.5 ±0.2 




0.1 ±0.2 


GB6 J1516+0015 


15 


17 


44 


-24 21 


205 


1.8 ±0.08 


1.8 ±0.1 


1.9 ±0.2 


1.8 ±0.3 


2.3 ± 1 


0.0 ±0.3 


PMN J1517-2422 


15 


40 


51 


14 47 




1.1 ±0.07 




0.8 ±0.2 


1.0 ±0.2 




-0.2 ±0.4 


GB6 J1540+1447 


15 


49 


27 


50 37 




1.0 ±0.07 


0.8 ±0.1 


0.9 ±0.1 


0.8 ±0.2 


1.1 ±0.3 


-0.1 ±0.3 


GB6 J1549+5038 


15 


49 


31 


02 36 


005 


2.2 ±0.08 


2.4 ±0.1 


2.1 ±0.1 


1.9 ±0.3 


1.9 ±0.5 


-0.1 ±0.2 


GB6 J1549+0237 


15 


50 


38 


05 26 


007 


2.4 ±0.07 


2.4 ±0.1 


2.0 ±0.1 


1.9 ±0.2 


1.9 ±0.4 


-0.2 ±0.2 


GB6 J1550+0527 


15 


56 


13 


-79 12 




0.8 ±0.07 


0.4 ±0.1 


0.5 ±0.1 






-1.1 ± 1 


PMN J1556-7914 


16 


01 


56 


33 29 




0.9 ±0.06 


0.8 ±0.08 


0.9 ±0.08 


0.7±0.2 




-0.1 ±0.4 


GB6 J1602+3326 


16 


08 


50 


10 28 


009 


2.3 ±0.07 


2.4 ±0.1 


2.2 ±0.1 


1.9 ±0.2 


2.1 ±0.7 


-0.1 ±0.2 


GB6 J1608+1029 


16 


13 


41 


34 12 


023 


3.9 ±0.06 


3.4 ±0.1 


3.3 ±0.1 


2.8 ±0.2 


1.8 ±0.4 


-0.3 ±0.1 


GB6 J1613+3412 


16 


17 


55 


-77 16 


183 


2.4 ±0.06 


2.2 ±0.09 


2.0 ±0.1 


1.7 ±0.2 




-0.3 ±0.2 


PMN J1617-7717 


16 


33 


15 


82 27 


076 


1.3 ±0.05 


1.5 ±0.08 


1.4 ±0.1 


1.4 ±0.3 




0.3 ±0.3 




16 


35 


16 


38 07 


033 


4.2 ±0.06 


4.9 ±0.1 


4.9 ±0.1 


4.5 ±0.2 


3.9 ±0.4 


0.1 ±0.07 


GB6 J1635+3808 


16 


37 


27 


47 14 




0.8 ±0.06 




0.9 ±0.1 


0.7 ± 0.2 




0.1 ±0.4 


GB6 J1637+4717 


16 


38 


12 


57 22 


056 


1.3 ±0.05 


1.4 ±0.1 


1.5 ±0.1 


1.6 ±0.2 


1.2 ±0.3 


0.1 ±0.2 


GB6 J1638+5720 


16 


42 


39 


68 53 


069 


1.0 ±0.06 


1.0 ±0.1 


0.9 ±0.1 


1.0 ±0.2 


0.8 ±0.3 


-0.2 ±0.3 


GB6 J1642+6856 


16 


42 


56 


39 48 


035 


6.9 ±0.06 


6.5 ±0.1 


6.0 ±0.1 


5.3 ±0.2 


4.0 ±0.4 


-0.3 ±0.06 


GB6 J1642+3948 


16 


43 


19 


-77 12 




0.8 ±0.08 


0.9 ±0.1 


0.8 ±0.1 


0.6 ±0.2 




0.0 ±0.5 


PMN J1644-7715 
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RA [hms] 


Dec [dm] 


ID 




K [Jy] 


Ka [Jy] 


Q [Jy] 


V[Jy] 




W [Jy] 


a 


5 GHz ID 


16 48 


26 


41 14 








0.8±0.1 


1.2 ±0.2 










1.6 ±2 


GB6 J1648+4104 


16 51 


06 


04 57 


010 


1 


6 ±0.09 


1.0 ±0.1 


1.1 ±0.2 


1 


0± 


3 




-0.7 ±0.4 


GB6 J1651±0459 


16 54 


12 


39 39 


036 


1 


1±0.05 


1.2 ±0.09 










1.0 ±0.4 


0.0 ±0.4 


GB6 J1653±3945 


16 56 


57 


57 06 







6 ±0.07 


0.8 ±0.1 


0.6 ±0.1 





9± 


2 


1.0 ±0.4 


0.3 ±0.4 


GB6 J1657+5705 


16 57 


54 


47 49 




1 


0±0.05 




0.6 ±0.09 










-1.0 ±0.6 




16 58 


08 


07 42 


013 


1 


1±0.07 


1.2 ±0.09 


1.0 ±0.1 


1 


5±0 


5 


1.2 ±0.6 


-0.0 ±0.4 


GB6 J1658+0741 


16 59 


40 


68 29 










0.7 ±0.09 


1 


0±0 


2 


0.8 ±0.3 


0.4 ±0.9 


GB6 J1700+6830 


17 03 


40 


-62 13 


198 


1 


8 ±0.06 


1.9 ±0.1 


2.0 ± 0.1 


1 


7±0 


2 




0.0 ±0.2 


PMN J1703-6212 


17 07 


39 


01 47 







9 ±0.08 


0.9 ±0.1 


0.9 ±0.1 





7±0 


2 




-0.0 ±0.5 


GB6 J1707+0148 


17 16 


02 


68 40 







6 ±0.05 


0.6 ±0.07 


0.7 ±0.07 





6±0 


1 




0.1 ±0.4 


GB6 J1716+6836 


17 24 


00 


-65 00 


196 


2 


3 ±0.06 


2.0 ±0.09 


1.7±0.1 


1 


2±0 


2 


1.2 ±0.4 


-0.5 ±0.2 


PMN J1723-6500 


17 27 


17 


45 30 


043 


1 


2 ±0.06 


1.3 ±0.09 


1.1 ±0.1 


1 


4±0 


2 


1.5 ±0.4 


0.1 ±0.2 


GB6 J1727+4530 


17 34 


13 


38 56 


038 


1 


1±0.06 


1.3 ±0.1 


1.2 ±0.1 


1 


3±0 


2 


1.1 ±0.4 


0.2 ±0.3 


GB6 J1734+3857 


17 37 


11 


-79 34 


186 





9 ±0.05 


1.1 ±0.08 


1.1 ±0.09 










0.3 ±0.4 


PMN J1733-7935 


17 40 


11 


47 39 







8 ±0.06 


0.9 ±0.08 


1.1 ±0.1 





8±0 


2 




0.2 ±0.4 


GB6 J1739±4738 


17 40 


33 


52 12 


048 


1 


3 ±0.06 


1.2 ±0.08 


1.3 ±0.1 


1 


5±0 


2 




0.1 ±0.2 


GB6 J1740+5211 


17 48 


57 


70 06 


068 





6 ±0.05 


0.7 ±0.06 


0.9 ±0.06 


1 


2±0 


1 




0.7 ±0.3 


GB6 J1748+7005 


17 53 


18 


44 05 







8 ±0.07 


0.6 ±0.1 


0.8 ±0.1 










-0.2 ±0.6 


GB6 J1753+4410 


17 53 


34 


28 48 


022 


2 


2 ±0.06 


2.3 ±0.09 


2.4 ±0.1 


2 


6±0 


2 


1.8 ±0.4 


0.1 ±0.1 


GB6 J1753+2847 


17 58 


52 


66 32 


064 





6 ±0.02 


0.6 ±0.03 


0.6 ±0.08 


0.4 ±0 


1 




0.0 ±0.3 


GB6 J1758+6638 


17 59 


53 


38 52 







9 ±0.06 


0.9 ±0.1 


0.8 ±0.1 










-0.2 ±0.6 


GB6 J1800+3848 


18 00 


24 


78 28 


072 


1 


8 ±0.06 


1.7 ±0.08 


1.7 ±0.1 


1 


6±0 


2 


1.0 ±0.2 


-0.2 ±0.2 


Uy 1803±78 


18 01 


33 


44 04 




1 


2 ±0.06 


1.3 ±0.08 


1.4 ±0.1 


1 


3±0 


2 


1.1 ±0.3 


0.2 ±0.2 


GB6 J1801±4404 


18 03 


02 


-65 07 


199 


1 


3 ±0.06 


1.3 ±0.09 


1.5 ±0.1 


1 


2±0 


3 


0.9 ±0.4 


0.1 ±0.3 


PMN J1803-6507 


18 06 


47 


69 49 


067 


1 


4 ±0.04 


1.4 ±0.07 


1.2 ±0.08 


1 


5±0 


1 


1.3 ±0.4 


-0.0 ±0.2 


GB6 J1806+6949 


18 20 


04 


-55 21 







9 ±0.08 


0.6 ±0.1 












-0.9 ±2 


PMN J1819-5521 


18 20 


05 


-63 43 


200 


1 


5 ±0.06 


1.4 ±0.1 


1.2 ±0.1 


1 


2±0 


2 


1.5 ±0.4 


-0.2 ±0.3 


PMN J1819-6345 


18 24 


07 


56 50 


053 


1 


5 ±0.05 


1.4 ±0.09 


1.4 ±0.09 


1 


5±0 


2 


1.3 ±0.4 


-0.1 ±0.2 


GB6 J1824+5650 


18 29 


41 


48 44 


046 


2 


6 ±0.06 


2.5 ±0.09 


2.5 ±0.1 


1 


8±0 


2 


1.6 ±0.3 


-0.2 ±0.1 


GB6 J1829±4844 


18 35 


00 


32 46 







8 ±0.07 


0.9 ±0.08 


0.7 ±0.09 










-0.0 ±0.6 


GB6 J1835+3241 


18 37 


29 


-71 06 


192 


1 


9 ±0.05 


1.8 ±0.07 


1.5 ±0.08 


1 


2±0 


1 




-0.4 ±0.2 


PMN J1837-7108 


18 40 


29 


79 46 


073 


1 


2 ±0.05 


0.9 ±0.09 


1.0 ±0.2 





5±0 


2 




-0.6 ±0.4 


Uy 1845±79 


18 42 


47 


68 07 


066 


1 


0±0.04 


1.2 ±0.06 


1.2 ±0.06 





9±0 


1 




0.1 ±0.2 


GB6 J1842+6809 


18 48 


30 


32 21 







7 ±0.06 


0.9 ±0.1 


0.8 ±0.1 


1 


3±0 


3 




0.5 ±0.5 


GB6 J1848+3219 


18 49 


34 


67 05 


065 


1 


3 ±0.05 


1.4 ±0.08 


1.4 ±0.06 


1 


3±0 


1 


1.5 ±0.2 


0.1 ±0.2 


GB6 J1849±6705 


18 50 


40 


28 23 


028 


1 


4 ±0.05 


1.2 ±0.08 


0.9 ±0.08 





7±0 


2 


1.0 ±0.5 


-0.6 ±0.3 


GB6 J1850±2825 


19 02 


56 


31 53 


034 


1 


3 ±0.05 


1.1 ±0.08 


0.8 ±0.08 





8±0 


2 




-0.6 ±0.3 


GB6 J1902+3159 
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RA [hms] 


Dec [dm] 


ID 


K [Jy] 


Ka [Jy] 


Q [Jy] 


V[Jy] 


W [Jy] 


a 


5 GHz ID 


19 


23 


30 


-21 06 


008 


2.1 ±0.08 


2.2 ±0.1 


2.1 ±0.1 


2.5 ±0.3 


2.0 ±0.8 


0.1 ±0.2 


PMN J1923-2104 


19 


27 


27 


61 18 


059 


1.1 ±0.05 


1.0 ±0.09 


1.1 ±0.1 


1.0 ±0.2 




-0.1 ±0.3 


GB6 J1927+6117 


19 


27 


37 


73 57 


070 


3.5 ±0.05 


3.4 ±0.08 


2.9 ±0.08 


2.8 ±0.2 


2.0 ±0.4 


-0.3 ±0.09 


GB6 J1927+7357 


19 


37 


05 


-39 56 




0.7 ±0.08 


0.9 ±0.2 


1.1 ±0.2 


1.3 ±0.3 




0.6 ±0.5 


PMN J1937-3957 


19 


38 


11 


-63 45 




0.8 ±0.06 


0.8 ±0.09 








-0.1 ±0.8 


PMN J1939-6342 


19 


•52 


17 


02 33 




0.9 ±0.09 


0.6 ±0.1 


0.6 ±0.1 






-0.8 ±0.8 


GB6 J1952+0230 


19 


55 


35 


51 32 


051 




1.1 ±0.1 


0.8 ±0.1 


1.1 ±0.3 




-0.3 ± 1 


GB6 J1955+5131 


19 


58 


00 


-38 45 


003 


3.1 ±0.08 


3.3 ±0.1 


3.0 ±0.1 


2.8 ±0.2 


2.7 ± 1 


-0.1 ±0.1 


PMN J1957-3845 


20 


00 


56 


-17 49 


Oil 


1.6 ±0.08 


1.5 ±0.1 


1.6 ±0.1 


1.6 ±0.2 


2.5 ± 1 


-0.0 ±0.3 


PMN J2000-1748 


20 


05 


57 


77 55 




0.7 ±0.06 






0.9 ±0.3 


1.0 ±0.3 


0.3 ±0.4 


Uy 2007+77 


20 


10 


23 


72 30 




0.7 ±0.08 


0.9 ±0.1 


0.8 ±0.09 


1.1 ±0.2 




0.4 ±0.5 


GB6 J2009+7229 


20 


11 


17 


-15 47 


014 


1.5 ±0.07 


1.6 ±0.1 


1.5 ±0.2 


1.1 ±0.3 




-0.1 ±0.4 


PMN J2011-1546 


20 


22 


30 


61 36 


063 


1.5 ±0.07 


1.4 ±0.09 


1.3 ±0.1 


0.9 ±0.2 




-0.3 ±0.3 


GB6 J2022±6137 


20 


35 


03 


-68 45 


194 


0.8 ±0.06 


0.9 ±0.1 


0.8 ±0.1 


0.7±0.2 




-0.0 ±0.5 


PMN J2035-6846 


20 


35 


22 


10 58 




0.7 ±0.08 


1.2 ±0.2 




1.0 ±0.2 


0.9 ±0.4 


0.3 ±0.4 


GB6 J2035+1055 


20 


5G 


07 


-47 17 


208 


1.7 ±0.06 


1.9 ±0.09 


1.9 ±0.1 


1.9 ±0.2 


2.4 ±0.9 


0.2 ±0.2 


PMN J2056-4714 


21 


09 


26 


-41 14 


001 


1.5 ±0.06 


1.7±0.1 


1.2 ±0.1 


0.9 ±0.2 


1.2 ±0.5 


-0.2 ±0.3 


PMN J2109-4110 


21 


09 


37 


35 37 


049 


1.0 ±0.07 


0.8 ±0.1 


0.7 ±0.1 




1.1 ±0.3 


-0.2 ±0.4 


GB6 J2109+3532 


21 


23 


43 


05 36 


027 


2.2 ±0.07 


1.8 ±0.1 


1.7 ± 0.1 


1.5 ±0.3 




-0.4 ±0.3 


GB6 J2123+0535 


21 


31 


33 


-12 06 


017 


2.7 ±0.08 


2.4 ±0.1 


2.3 ±0.1 


1.9 ±0.2 




-0.3 ±0.2 


PMN J2131-1207 


21 


34 


05 


-01 54 


020 


1.9 ±0.08 


1.8 ±0.1 


1.7 ±0.1 


1.5 ±0.3 


1.1 ±0.4 


-0.2 ±0.3 


PMN J2134-0153 


21 


3G 


37 


00 42 


025 


4.3 ±0.08 


3.5 ±0.1 


3.0 ±0.2 


1.4 ±0.2 


1.3 ±0.4 


-0.7 ±0.2 


GB6 J2136±0041 


21 


39 


17 


14 24 


041 


2.3 ±0.07 


2.2 ±0.1 


1.9 ±0.1 


1.6 ±0.2 




-0.3 ±0.2 


GB6 J2139+1423 


21 


43 


25 


17 40 


044 


1.2 ±0.06 


1.4 ±0.09 








0.5 ±0.5 


GB6 J2143+1743 


21 


48 


05 


06 57 


037 


7.8 ±0.07 


7.5 ±0.1 


7.2 ±0.1 


6.6 ±0.2 


5.1 ±0.5 


-0.2 ±0.06 


GB6 J2148±0657 


21 


48 


49 


-77 57 


184 


1.4 ±0.05 


1.3 ±0.09 








-0.3 ±0.5 


PMN J2146-7755 


21 


51 


45 


-30 26 




1.0 ±0.07 


1.2 ±0.1 


1.2 ±0.1 


1.5 ±0.2 


1.0 ±0.5 


0.3 ±0.3 


PMN J2151-3028 


21 


57 


05 


-69 42 


190 


3.6 ±0.06 


2.9 ±0.1 


2.6 ±0.1 


2.1 ±0.2 


1.5 ±0.7 


-0.6 ±0.1 


PMN J2157-6941 


21 


58 


07 


-15 02 


018 


2.0 ±0.08 


2.0 ±0.1 


1.7 ±0.1 


1.5 ±0.3 




-0.2 ±0.3 


PMN J2158-1501 


22 


02 


49 


42 17 


058 


2.9 ±0.06 


3.0 ±0.09 


3.1 ±0.09 


2.7±0.2 




0.1 ±0.1 


GB6 J2202+4216 


22 


03 


20 


31 46 


054 


2.8 ±0.06 


2.7±0.1 


2.2 ±0.1 


1.8 ±0.3 


1.8 ±0.4 


-0.4 ±0.2 


GB6 J2203+3145 


22 


03 


26 


17 24 


045 


1.5 ±0.08 


1.7±0.1 


1.7 ± 0.1 


1.7 ±0.2 




0.2 ±0.3 


GB6 J2203+1725 


22 


0G 


09 


-18 39 


016 


1.7±0.07 


1.7±0.1 


1.3 ±0.1 


1.5 ±0.4 




-0.3 ±0.3 


PMN J2206-1835 


22 


07 


08 


-53 49 




1.0 ±0.06 


0.9 ±0.09 


0.8 ±0.1 






-0.3 ±0.6 


PMN J2207-5346 


22 


11 


33 


23 51 


050 


1.0 ±0.07 




1.4 ±0.1 


1.0 ±0.3 




0.5 ±0.4 


GB6 J2212±2355 


22 


18 


52 


-03 35 


030 


2.4 ±0.08 


2.0 ±0.1 


2.0 ±0.2 


1.7 ±0.2 




-0.4 ±0.2 


PMN J2218-0335 


22 


25 


46 


-04 56 


029 


4.9 ±0.07 


4.6 ±0.1 


4.1 ±0.1 


3.8 ±0.3 


5.4 ±2 


-0.2 ±0.1 


PMN J2225-0457 
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Table 9 — Continued 



RA [hmsj 


Dec [dm] 


ID 


I - r T -i 

K [Jy] 


TS~ r T 1 

Ka [Jy] 


/-\ r t 1 

Q [Jy] 


V [Jy] 


tit r t 1 

W [Jy] 


a 


5 GHz ID 


oo oq 41 

Z Z Z C? t: J. 


—08 32 

UO OZ/ 


024 


o + n f)8 

Z .O —1— VJ ,UO 


2 7 + 1 

Z . I _l_ 0. J- 


2.6 ± 0.2 


q 4- 3 

O . Z _l_ O.O 


2.2 ± 0.4 


0.2 ± 0.2 


PMN 12229-0832 

J. 1V_L1>I O Z Z LitJ ' OOO Z 


22 2Q 4fi 

Z Z Zj , J 1 U 


—20 51 




9 + 08 

\J . O _l_ o.oo 


0.8 ± 0.1 


0.8 ± 0.1 


1.1 ± 0.4 




—0.1 ± 0.6 


PMN T2229-2049 


22 32 38 

_Z OZj (JO 


11 44 


047 


Q o _|_ n no 
O.O —1— o.oo 


3.3 ± 0.1 


3.4 ± 0.1 


3.7 ± 0.3 


2.8 ± 0.3 


0.1 ± 0.1 


CR6 T2232+1 143 

V..I 1 J\) O ZjZjOZj 1 J LtU 


22 35 12 

i Z OO -L Z 


—48 34 


206 


1.7 ± 0.06 


1.9 ± 0.1 


1.8 ± 0.1 


1.9 ± 0.2 


2.1 ± 0.5 


0.1 ± 0.2 


PMN T2235-4835 

J. XVJ.1 > (J ilj ilj O O tiOOO 


22 3fi 18 


28 25 


057 


i o + n 08 

J. • Z i o.oo 


1.3 ± 0.1 


1.1 ± 0.1 






—0.0 ± 0.4 


CR6 T2236+2828 

VJJJU O ZjZjOO I zozo 


22 3Q 39 

Z;Z; Ul/ O O 


—57 00 


201 


1.1 ± 0.05 


1.3 ± 0.07 








0.5 ± 0.4 


PMN 12239-5701 

_L 1V±1> O Z ZO iJ OIOJ- 


22 4fi 13 

ZZ iO J_ O 


— 12 10 


021 


1.2 ± 0.07 


1.2 ± 0.2 


1.1 ± 0.2 


1.0 ± 0.4 




—0.2 ± 0.6 


PMN 12246-1206 

i. XVJ.1 > (JZZtiO J-ZOO 


22 54 00 


16 08 


055 


7.0 ± 0.07 


6.7 ± 0.1 


6.8 ± 0.1 


6.3 ± 0.2 


5.6 ± 0.5 


—0 1 ± 06 

\f . -L _1_ O.OO 


GR6 12253+1608 

VI 1 J \ J tj ZjZjOO 1 -LOOO 


22 55 3fi 

Z;Z; OO 


42 02 




1.1 ± 0.05 


8 + 07 

O . O —1— U.UI 


0.9 ± 0.2 


0.5 ± 0.2 




—0.6 ± 0.5 


GR6 T2255+4202 

V I 1 J \) t j ZjZjOO 1 TiZOZ 


22 56 18 

ZiZi uu io 


—20 11 


01 9 


n q + n 08 

O . Zj —1— O.OO 


7 + 1 

O. 1 —1— O. -L 


8 + 01 

O.O —1— o . ± 


1.1 ± 0.4 




_0 2 + 6 

O .Zi _l_ O. U 


PMN T2256-201 1 

J. 1VJ.1>I U ZZOO ZO-L J- 


22 58 Ofi 

Z;Z; OO OVJ 


—27 57 


012 


6.7 ± 0.07 


6.7 ± 0.1 


6.6 ± 0.1 


5.9 ± 0.3 


5.2 ± 0.5 


_0 1 + 06 

O.J- —1— o.oo 


PMN T2258-2758 

±. XVJ.1 > O Z ZOO ZIOO 


23 15 57 


—50 17 


204 


i o + o 05 

J. • Z i o.oo 


i o + o 08 

J. • Z i o.oo 


q + o 09 

\l . <J _1_ O.Ot/ 






—0.2 ± 0.3 


PMN T2315-5018 

_L XV ±1 > (J ilj O O OOJ.O 


23 22 54 

ZjO ZjZj Ot: 


51 07 




8 + 06 

O.O —1— o.oo 


7 + o 09 

O.I —1— O . \J-J 


0.8 ± 0.1 






—0.3 ± 0.6 


GR6 T2322+5057 a 


23 27 42 

— O Z I t:Z 


09 39 




8 + 08 

O.O —1— o.oo 


1.3 ± 0.2 


1.2 ± 0.2 


0.9 ± 0.3 


1.1 ± 0.5 


0.5 ± 0.5 


GR6 T2327+0940 

V. I 1 J\ > O Z;OZ; 1 \ UJt:U 


23 2Q 10 

ZjO Zj ij J- O 


—47 32 




1.4 ± 0.05 


1.2 ± 0.1 


1.4 ± 0.1 


1.3 ± 0.2 


1.2 ± 0.3 


—0.1 ± 0.2 


PMN T2329-4730 

J. XVJ.1 > O ZO Z O t:IOO 


23 30 52 

Z;0 OO oz 


10 56 




1.0 ± 0.07 


1.2 ± 0.1 


1.0 ± 0.1 


1.3 ± 0.5 




0.1 ± 0.4 


GR6 T2330+1 100 

V. I 1 )\> o Z O O VJ \ J- J- oo 


23 31 19 


-15 59 


032 


1.3 ±0.08 


0.9 ±0.1 


0.8 ±0.2 


1.1 ±0.5 




-0.5 ±0.6 


PMN J2331-1556 


23 33 43 


-23 40 




0.9 ±0.06 


0.9 ±0.1 


1.2 ±0.2 


1.2 ±0.2 




0.3 ±0.4 


PMN J2333-2343 a 


23 34 11 


07 34 




1.1 ±0.08 


1.0 ±0.1 


1.1 ±0.1 


0.8 ±0.2 




-0.2 ±0.4 


GR6 J2334+0736 


23 35 05 


-01 29 




0.7 ±0.08 


1.2 ±0.1 


1.3 ±0.1 


0.7±0.2 




0.7±0.6 


PMN J2335-0131 


23 35 30 


-52 42 


195 


1.0 ±0.04 


0.7 ±0.07 


0.5 ±0.1 


0.9 ±0.3 




-0.8 ±0.5 


PMN J2336-5236 


23 46 42 


09 27 




1.1 ±0.07 


1.0 ±0.1 


0.7±0.1 


1.0 ±0.4 


1.1 ±0.5 


-0.3 ±0.4 


GR6 J2346+0930 a 


23 48 16 


-16 31 


039 


1.8 ±0.08 


1.8 ±0.1 


2.0 ±0.1 


1.6 ±0.2 




0.1 ±0.2 


PMN J2348-1631 


23 54 25 


45 50 


074 


1.6 ±0.06 


1.2 ±0.1 


1.3 ±0.1 


1.2 ±0.2 




-0.4 ±0.3 


GR6 J2354+4553 


23 56 21 


49 52 


075 


0.9 ±0.04 


0.8 ±0.07 








-0.3 ±0.6 


g 


23 57 54 


-53 13 


189 


1.3 ±0.05 


1.1 ±0.09 


1.1 ±0.1 




0.8 ±0.3 


-0.4 ±0.3 


PMN J2357-5311 


23 58 08 


-10 13 




1.0 ±0.07 


1.2 ±0.1 




1.1 ±0.2 




0.2 ±0.4 


PMN J2358-1020 


23 58 52 


-60 50 


187 


1.8 ±0.06 


1.4 ±0.09 


1.2 ±0.09 






-0.7 ±0.3 


PMN J2358-6054 



indicates the source has multiple possible identifications. 

b Source J0322-3711 (Fornax A) is extended, and the fluxes listed were obtained by aperture photometry. 
c Source J0519-0540 is a blend of the Lynds Rright Nebulae LRN 207.65-23.11 and LRN 207.29-22.66. 

d Source J1357+7644is outside of the declination range of the GR6 and PMN catalogs. Identified as QSO NVSSJ135755+764320 
by Trushkin (2006, private communication). 

e Source J1633+8227 is outside of the declination range of the GR6 and PMN catalogs. It was identified as NGC 6251 by 
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Trushkinl (|2003h . 

f Source J1657+4749 is identified as QSO GB6J1658+4737 by Trushkin (2006, private communication). Offset from the WMAP 
position is 11.6 arcminutes. 

g Source J2356+4952 is identified as GB6 J2355+4950 by Trushkin (2003). Offset from the WMAP position is 11.9 arcminutes. 




Fig. 1. — Full-sky maps in Galactic coordinates smoothed with a 0?2 Gaussian beam, shown 
in Mollweide projection, top left: K-band (23 GHz), middle left: Ka-band (33 GHz), bottom 
left: Q-band (41 GHz), top right: V-band (61 GHz), bottom right: W-band (94 GHz). 
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Fig. 2. — The number of independent W-band observations per pixel in Galactic coordinates, 
top: year-1, middle: year-2, and bottom: year-3. The number of observations is greatest near 
the ecliptic poles and in rings approximately 141° from each pole (determined by the angular 
separation between the two bore-sight directions). The number of observations is least in 
the ecliptic plane. The small circular cuts in the ecliptic are where Mars, Saturn, Jupiter, 
Uranus, and Neptune are masked so as not to contaminate the CMB signal. The coverage 
is quite consistent from year to year, with the planet cuts being responsible for the largest 
fractional variation. 
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Fig. 3. — Comparison of the three-year maps with the previously released one-year maps. 
The data are smoothed to 1° resolution and are shown in Galactic coordinates. The frequency 
bands K through W are shown top to bottom. The first-year maps (left) and the three-year 
maps (middle) are shown scaled to ±200 /xK. The difference maps (right) are degraded to 
pixel resolution 4 and scaled to ±20 fiK. The small difference in low-Z power is mostly due 
to im provements in the gain model of the instrument as a function of time (IJarosik et al. 
2006h . See g3]and Table EJ 
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Fig. 4. — Results from the MEM foreground degeneracy analysis, (a) The spectrum of 
total foreground emission (diamonds) compared to the sum of the MEM components (lines 
connecting the diamonds), averaged over the full sky. The three component spectra of 
the model are also shown, as indicated, (b-f) The contour plots illustrate the degeneracy 
between selected components. The panels show the change in the MEM functional (in units 
of Ax 2 ) as a function of two global model parameters. For example, panel (b) shows the 
change in the MEM functional that results from adding a constant value to the synchrotron 
and/or free-free solutions, while panel (f) shows the dependence of the functional on a global 
steepening of the synchrotron spectrum, modeled as (3 a = /3 + djSs/dlogv ■ (logz/ — logi^). 
The only strong degeneracy is between the synchrotron and free-free amplitudes in (b). This 
effect is mitigated by the prior distributions assumed for each component (see text). The 
dust spectral index, (3^, was found to be essentially unconstrained, so contour plots for that 
quantity are not shown. 
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Fig. 5. — Galactic signal component maps from the Maximum Entropy Method (MEM) 
analysis ( §4.51) . top-bottom: synchrotron, free-free, and dust emission with logarithmic tem- 
perature scales, as indicated, left: Input prior maps for each component, right: Output 
maps based on three-year WMAP data for each component. See text for discussion. 
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Fig. 6. — The mean Galaxy spectrum from 408 MHz to 94 GHz, based on data from Haslam, 
Tenerife and WMAP. Each point is the mean signal in the Galactic latitude range 20° < 
|fe| < 50° with the mean at the Galactic poles subtracted. The WMAP values (23-94 GHz) 
were derived from the MEM dust model (blue), the MEM free-free model (green), and the 
MEM non-thermal component (red), and the sum of the three (black). The Haslam (408 
MHz) value is computed directly from the map (black). The dashed black curve is the 
interpolation of the total signal between 408 MHz and 23 GHz. The dashed green curve is 
the extrapolation of the WMAP MEM free-free result to 408 MHz assuming (3 S = -2.14. 
Th e Tenerife values were obtaine d by taking the template correlatio n coefficients reported 
bv lde Oliveira-Costa et all J2002h and lde Oliveira-Costa etaJ j^OoJ (fit to \b\ > 20° data), 
scaling the templates to obtain model emission maps, then evaluating the mean signal in the 
same way as with the WMAP and Haslam data. The Tenerife correlation with Ha provided 
only upper limits for the free-free signal. See §4.51 for more detail. 
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Fig. 7. — An overview of the microwave sky. The yellow, salmon, and red shaded regions 
indicate the KpO, Kp2, and Kp8 diffuse em i ssion m asks where the Galactic foreground 
signal is especially strong. See iBennett et al.l (12003a ) for a discussion of how these masks 
were constructed. The KpO and Kp2 masks are useful for cosmoloKical analysis , and the 
Kp8 mask closely follows the "processing" mask described in iJarosik et al.l (120061 ). which is 
used for reducing systematic errors in t he sky maps . The violet shading shows the "P06" 
polarization analysis mask described in iPage et al.l (120061 ). The small blue dots indicate 
the position of point sources detected by WMAP. Some well-known sources and regions are 
specifically labeled. 
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Fig. 8. — top: Full-sky map color-coded to show the 12 regions that were used to generate 
the three-year ILC map (see §5.2|) . bottom: The mean ILC residual map from 100 Monte 
Carlo simulations of CMB signal, Galactic foreground signal and instrument noise. The 
CMB signal was drawn from a ACDM power spectrum that was modified to reproduce the 
power measured in the first-year spectrum at I = 2, 3. The first-year MEM foreground model 
was used to generate the (fixed) Galactic maps. In the simulations, systematic errors arose 
in the recovered ILC maps due to a combination of effects: (i) a tendency for the minimum 
variance method to exploit (anti) alignments between the CMB and foreground signal, and 
(ii) variations in the spectra of the foreground signal, due mostly to changes in the admixture 
of synchrotron, free-free, and dust emission (see §5.21) . This bias map was used to correct 
the three-year ILC map (Figure [9]). 
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Fig. 9. — top: The first-year ILC map reproduced from iBennett et al.l (12003d ). middle: The 
three-year ILC map produced following the steps outlined in §5.21 bottom: The difference 
between the two (1-yr — 3-yr). The primary reason for the difference is the new bias 
correction ( Figure El). The low-/ change noted in 33] and shown in Figure [3] is also apparent. 
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Fig. 10. — Galactic foreground removal with spatial templates. All maps in this figure are 
three-year maps that have had the ILC estimate of the CMB signal subtracted off to highlight 
the foreground emission. The maps have been degraded to pixel resolution 5, are displayed in 
Galactic coordinates, and are scaled to ±30 fiK. The white contour indicates the perimeter of 
the Kp2 sky cut, outside of which the template fits were evaluated. The frequency bands Q 
through W are shown top to bottom, (left) Sky maps prior to the subtraction of the best-fit 
foreground model ( §5.3p . (middle) The same sky maps with the first-year template-based 
model subtracted. Note the high-latitude residuals in the vicinity of the North Polar Spur 
and around the inner Galaxy due to the use of the Haslam 408 MHz map as a synchrotron 
template, (right) The same sky maps with the three-year template-based model subtracted. 
This model substitutes K- and Ka-band data for the Haslam data which produces lower 
residuals outside the Kp2 sky cut. There are still isolated spots with residual emission of 
order 30 //K in the vicinity of the Gum Nebula and the Ophiuchus Complex (see Figure [7]). 
Note also that substantial errors (>30 fiK) remain inside the Kp2 cut due to limitations in 
the template model. 
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Fig. 11. — Measurement of the source number count distribution dN/dS. red diamonds: from 
the WMAP three-year point source catalog at 40.7 GHz (Q -band). stars: fro m the Australian 
Telescope Compact Array (ATCA) 18 GHz pilo t survey (IRicci et al.ll2004j ); triangles: from 
the 9C survey at 1 5 GHz (IWaldram et al.ll2003l ); and squares: from the Very Small Array 



(VSA) at 33 GHz ( Cleary et al. 20051 ) . The parallelogram is from the Cosmic Background 



Imager (CBI) experiment at 31 GHz (IMason et al.ll2003l ). Several models for dN/dS are 



shown for comparison: the solid curve is the 44 GHz dN/ dS model from iToffolatti et al. 



(119981 ) . the dotted curve is the Toffolatti model rescaled by 0.66 (found to b e a good fit to 
the fir st-year data), and the dashed curve is an updated 40.7 GHz model from lde Zotti et al. 
(l2005h . 
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Fig. 12. — The predicted Ci uncertainty (inverse Fisher matrix) at low I, in /zK 4 . The red 
curves show the result when pixel-pixel noise correlations are ignored. The smooth rise at low 
I reflects our approximate representation of 1/f noise in the noise bias model, equation (j42j) . 
The black curves account for the full structure of the pixel-pixel inverse covariance matrix, 
including 1/f noise and map- making covariance. For the TT spectrum (left pair of curves 
in each panel), the noise is negligible compared to the si gnal, so this stru cture can be safely 
ignored. H o wever , the TE analysis must account for it (IPage et al.ll2006l ). See Figure 16 of 
Page et al.l (120061 ) for an analogous plot of the EE and BB uncertainties. 
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Fig. 13. — (top) Q-, V-, & W-band cross-power spectra before foreground subtraction, eval- 
uated outside the Kp2 sky cut, relative to the final combined spectrum, (bottom) Same 
cross-power spectra after subtracting a template-based diffuse foreground model ( §5.3[) and 
the best-fit residual point source contamination ( §7.21) . 
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Fig. 14. — Maps of power spectrum modes I = 2 — 8 computed from full-sky fits to the ILC 
map, shown at top left. Many authors note peculiar patterns in the phase of these modes, 
and many claim that the behavior is inconsistent with Gaussian random-phase fluctuations, 
as predicted by inflation. For example, the I = 5 mode appears strikingly symmetric (a 
non-random distribution of power in m), while the I = 2 and 3 modes appear unusually 
aligned. The significance of these a posteriori observations is being actively debated. See §8] 
for a more detailed discussion. 
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Fig. 15.— The posterior likelihood of 1(1 + 1)Q/2tt given the WMAP data, for I = 2 - 10. 
The curves are computed using data from the ILC map evaluated outside the Kp2 sky cut. 
The vertical blue lines show the values inferred using the pseudo-C/ estimate. These values 
are well within the posterior likelihood distribution, but there is a tendency for them to be 
lower than the peak, especially for I = 2, 3, 7. The vertical red lines show the predicted 
power spectrum from the best- fit ACDM model (fit to WMAP data). These points are well 
within the posterior likelihood in all cases; the faint dashed lines indicate the 68% and 95% 
confidence regions of the distribution (see §7.4p . Note that the curves approach a Gaussian 
distribution as I increases. 
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Fig. 16. — The binned three-year angular power spectrum (in black) from I — 2 — 1000, 
where it provides a cosmic variance limited measurement of the first acoustic peak, a robust 
measurement of the second peak, and clear evidence for rise to the third peak. The points 
are plotted with noise errors only (see text). Note that these errors decrease linearly with 
conti nued observing time . The red curve is the best-fit ACDM model, fit to WMAP data 
only (ISpergel et al.l I2006T ) . and the band is the binned la cosmic variance error. The red 
diamonds show the model points when binned in the same way as the data. 
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Fig. 17. — x 2 vs - I f° r the full power spectrum relative to the best- fit ACDM model, fit to 
WMAP data only. The x 2 per I has been averaged in /-bands of width Al = 15. The dark to 
light grey shading indicates the 1, 2, and 3cr confidence intervals for this distribution. The 
dashed line indicates the mode. 
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Fig. 18. — The WMAP three-year power spectrum (in black) compared t o other recent mea - 
sureme nts of the CMB angula r power spectrum, in cluding; Boo merang (IJones et al. | 2005h . 



Acbar flKuo et al.ll2004f ). CBI flReadhead et al.ll2004f ). and VSA ([Dickinson et al 



200J). For 



clarity, the I < 600 data from Boomerang and VSA are omitted; as the measurements are 
consistent with WMAP, but with lower weight. These data impressively confirm the turnover 
in the 3rd acoustic peak and probe the onset of Silk damping. With improved sensitivity 
on sub-degree scales, the WMAP data are becoming an increasingly important calibration 
source for high-resolution experiments. 
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Fig. 19. — Comparison of the three-year angular power spectrum with the first-year result. 
(top) The three-year combined spectrum, in black, is shown with the first-year spectrum, as 
published, in red. The best-fit three-year ACDM model is shown in grey for comparison. For 
/ < 30, the difference is due to a change in estimation methodology: the three- year spectrum 
is based on a maximum likelihood estimate, while the first-year results is based on a pseudo- 
C[ estimate. For I > 100 the difference is due primarily to (i) an improved determination of 
the WMAP beam response, and (ii) the improved sensitivity of the three-year data, (middle) 
Ratio of the three-year spectrum to the first-year spectrum. For Z < 30 we plot the ratio 
of the two pseudo-Crbased spectra to show the consistency of the underlying data. The 
red curve is the ratio of the first-year window function to the three-year window function. 
(bottom) Same as top panel except that the first-year spectrum has been multiplied by the 
window function ratio depicted in the middle panel, and the maximum likelihood estimate 
has been substituted for / < 30. 
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Fig. 20. — The unbinned three-year angular power spectrum (in black) from 1 = 2 — 400, 
where it provides a cosmic variance limited measurement of the first acoustic peak. The 
first-year spectrum is shown in green for comparison. The red curve with the grey band is 
the best-fit ACDM model and la error band per I. The width of the band is dominated 
by cosmic variance to / = 400. One feature that was singled out in the first-year spectrum 
was the "bite" at I ~ 208. The feature is still visible in the three-year spectrum, but not 
prominently. We believe this feature was predominantly a noise fluctuation in the first-year 
data. 
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Fig. 21. — Full-sky maps in ecliptic coordinates smoothed to 1° resolution, shown in Lambert 
Azimuthal Equal Area projection, top: three-year ILC map, bottom: three-year W-band map 
with the Kp2 Galaxy mask superposed, but no other foreground removal applied. Note the 
qualitative differences in large-scale power between the two hemispheres, as has been noted 
by several authors (see text). The statistical significance of this difference continues to be 
an area of active study. 
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Fig. 22. — Angular power spectra Cf T & Cf E from the three-year WMAP data, top: The 
TT data are as shown in Figure [161 The TE data are shown in units of 1(1 + l)Cj/27r, on the 
same scale as the TT sign al for comparis o n, bott om: The TE data, in units of (I + l)Cj/27T. 
This updates Figure 12 of iBennett et al.l (j2003bl ). 
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Fig. 23. — An estimate of the multiplicative error introduced in the angular power spectrum, 
Cjf, due to the effects of beam asymmetry. See §7.1.11 and Appendix [B] for details on how 
this estimate was obtained. Note that the final power spectrum does not include data from 
Q band. 



